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ABSTRACT 


This  thesis  investigates  the  estimation  of  Lanchester 
attrition- rate  coefficients  in  a  deterministic  aggregated 
combat  model.  Numerical  values  for  the  attrition-rate  coef- 
ficients in  the  Lanchester-type  model  are  taken  to  be  the 
same  as  those  for  a  corresponding  continuous-time  Markov- 
chain  model,  and  maximum  likelihood  estimators  are  developed 
for  these  Markov-chain  coefficients.   A  discussion  of  the 
historical  background  of  Lanchester* s  equations  preceeds  the 
thaoratical  development  of  components  of  the  attrition  rate 
equations.   Using  the  functional  forms  and  procedures  devel- 
oped by  Gordon  M.  Clark  in  his  doctoral  thesis,  a  simple 
computer  simulation  program  is  developed  for  a  short  dura- 
tion battle  consisting  of  homogeneous  forces.   Recommenda- 
tions are  made  to  modify  this  program  to  model  heterogeneous 
forces  in  a  similar  battle.   A  discussion  of  the  theoretical 
background  supports  these  recommendations. 
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I-   INTHODUCTIQN 

The  purpose  of  this  thesis  is  to  investigate  methods  for 
estimating  Lanchester  attrition-rate  coefficients,  used  in  a 
deterministic  aggregated  combat  model,  from  data  generated 
by  the  detailed,  high  resolution  Monte  Carlo  combat  simula- 
tion model  called  STAR  (Simulation  of  Tactical  Alternative 
Responses) .   An  additional  purpose  is  to  construct  a  simple 
Lanchester-type  combat  model  capable  of  using  the  estimated 
attrition  rate  coefficients  to  predict  the  outcome  of  STAR 
simulations. 

STAR  is  a  model  which  has  been  developed  primarily  by 
U.S.  Army  students  attending  the  Naval  Postgraduate  School. 
It  has  recently  gained  acceptance  in  the  Army  Community  as  a 
viable  combat  model,  useful  for  studies  such  as  the  United 
States  Army  Mortar  Study  to  be  conducted  in  1982.   The  cur- 
rent version  of  STAR  is  written  in  SIMSCEIPT  and  requires 
about  2  megabytes  of  core  storage.   The  run  time  for  one 
replication  of  STAR  involving  a  company  size  force  defending 
against  a  battalion  size  attacking  force  is  approximately 
four  CPU  minutes  on  the  IBM  3033  computer. 

The  basic  conceptual  difficulties  found  in  this  thesis 
are  summarized  in  Figure  1.   The  numerical  values  for  the 


attrition  rate   coefficients   in    the   deterministic   Lanchester- 
type   lodel   are    assumed  to   be  the  same   as   those  for  a    corre- 
sponding continuous-time   Markov  chain  model.      With  this 
assumption,    maximum  likelihood   estimators    (MLE*s)    can    be 
developed   for   the   Markov   chain   model   coefficients.      Specific 
data   needed   for   the   MLE   computations   is   extracted   from  the 
high  resolution   Monte   Carlo   simulation   and  used   to  compute 
•♦estimated"   coefficient   values.      These    values   for   the 
continuous   time    Markov  chain   model   are   then    substituted  into 
the   "analogous"    deterministic   Lanchester-type   model. 

Constructing   a   Lanchester-type   attrition    model  such  as 
that   proposed   in  this    thesis   is   of   significant  value.      With 
respect    to   the    hierarchy   of    models   concept  [ Ref .     1],    it 
could    possibly    fill   the  requirement    as   a    'hook*    or   link 
between    STAR   and  low    resolution    models   that    exist   at    higher 
levels    in  the   pyramid    of    models.      However,    possibly  the  most 
important  contribution  this   model  can  make   is   as    a   filter 
model    for   STAR   to   reduce   costs    and  time    to   run  the  many  var- 
iations  of   a   scenario    in   STAR.      To   explain    further,    consider 
the   upcoming   mortar  study.      To    conduct   such    a  study    many 
runs/replications   of   a  battle   using   the    STAR    model  will  have 
to  be    made.      These   numerous   runs  are   required   to    compare 
various   types   of   mortars,    various  types    of   force    mixes. 
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Figure    1:      Relationships    of   Models 
weapon    aixes,    and   scenarios.      Since   3TAH    is    a   aonrs   Carlo 
simulation,    Tariance-reducing   replications   are  required. 
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This  raquires  significant  computer  time  at  a  significant 
cost.   By  using  a  simpler  analytical  model  which  produces 
comparable  output  to  a  STAR  run,  many  force  mixes/weapon 
mixes/scenarios  can  be  eliminated  from  further  consideration 
at  a  minimal  cost.   STAR  can  then  be  used  to  concentrate 
only  on  the  more  promising  options  for  force/weapons  mixes. 

This  thesis  records  the  progress  made  in  developing  the 
small,  simple,  Lanchester-type  combat  model  discussed  above. 
It  also  lays  the  foundation  for  future  work  on  this  model. 
Chapter  2  of  this  thesis  provides  an  introduction  to  Lan- 
chester-type combat  models.   Chapter  3  outlines  the  develop- 
ment of  the  Deterministic  Aggregated  Model  of  STAR  on  the 
Apple  computer  (DAMSTAC) ,  which  represents  the  first  srep  in 
building  an  analytical  model  comparable  to  STAR.   Chapter  4 
explores  various  Lanchester  attrition  rate  functional  forms 
which  may  be  used  in  the  DAMSTAC  or  in  future  generation 
models.   Chapter  5  details  some  desired  enhancements  of  the 
model,  while  Chapter  6  discusses  the  work  required  to  tran- 
sition from  the  current  homogeneous  force  DAMSTAC  model  to  a 
more  realistic,  albeit  more  complicated  heterogeneous  model. 
Chapter  7  contains  the  concluding  remarks. 

Appendix  A  lists  the  variables  used  in  the  computer  pro- 
gram for  the  DAMSTAC  model.   Appendix  B  lists  some  notes  of 
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interest  concerning  the  program  and  model,  which  may  not  be 
apparent  to  the  reader.   Appendix  C  discusses  the  utility 
programs  developed  to  support  the  model.   Appendix  D  is  a 
discussion  of  additional  work  done  by  the  authors  in  the 
form  of  a  computer  model  that  explicitly  plays  the  terrain 
effect.   Appendix  E  consists  of  work  done  by  two  other  stu- 
dents on  a  postprocessor  which  supports  the  computational 
requirements  of  the  model  and  acts  as  an  interface  between 
STAR  and  the  model.   Finally,  Appendix  F  contains  the  compu- 
ter listings  of  all  programs  written  for  this  thesis. 
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II.   HISTORICAL  BACKGROUND 

Lanchester-type  combat  models  are  analytical  models  in 
which  differential  equations  are  used  to  express  the  rate  of 
change  of  force  levels  to  represent  the  effects  of  the 
attrition  process.   The  concept  of  representing  attrition 
effects  in  land  warfare  with  such  mathematical  equations  was 
originated  by  F.  W.  Lanchester  in  1914  in  order  to  quantita- 
tively justify  concentration  of  forces. 

Lanchester  hypothesized  that  under  conditions  of  "modern 
warfare"  attrition  between  two  homogeneous  forces,  each 
using  "aimed"  fire,  could  be  modelled  by: 


dx/dt  =  -ay  (t)    with  x(0)=  x« 


dy/dt  =  -bx(t)    with  y(0)=  y^ 


where  x (t)  and  y  (t)  are  the  number  of  X  and  Y  firers  alive 
at  time  t,  and  a  and  b  are  positive  constants  that  are 
called  Lanchester  attrition- rate  coefficients.   The  constant 
•a*  represents  the  effectiveness  of  an  x  firer  in  killing  y 
targets  per  unit  of  time  and  • b*  represents  the  equivalent 
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effectiveness  of  the  y  firers.   The  equations  written  above 
are  usually  called  "Lanchester • s  equations  for  modern  war- 
fare".  Clearly,  these  functional  forms  for  attrition  rates 
that  Lanchester  developed  are  directly  proportional  to  the 
attrition  rate  coefficients  and  the  number  of  firers  from 
the  opposing  force.   The  assumptions  that  must  be  made  in 
order  to  use  Lanchester's  aimed  fire  equations  to  model  com- 
bat attrition  are  very  restrictive.   These  assumptions  can 
be  stated  as  follows  [Ref.  2]: 

1.  Two  homogeneous  forces  are  engaged  in  combat; 

2.  Each  unit  on  either  side  is  within  weapon  range  of  all 
units  on  the  other  side; 

3.  The  effects  of  successive  rounds  in  the  target  area 
are  independent; 

4.  Each  unit  is  sufficiently  well  aware  of  the  location 
and  condition  of  all  enemy  units  so  that  it  engages  only 
live  enemy  units  (one  at  a  time)  and  kills  them  at  a 
constant  rate,  which  does  not  depend  on  the  enemy  force 
level.   When  an  enemy  target  is  Jcilled,  search  begins  for 
a  new  target,  with  the  rate  of  acquiring  a  new  enemy  tar- 
get being  independent  of  the  enemy  force  level; 

5.  Fire  is  uniformly  distributed  over  surviving  enemy 
units. 
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Lanchester  also  hypothesized  that  when  both  sides  used 
"area"  fire,  the  attrition  of  forces  could  be  modelled  by: 

dx/dt  =  -axy 
dy/dt  =  -bxy 
The  first  three  conditions  under  which  Lanchester 's 
equations  for  area  fire  apply  are  identical  to  those  for  the 
equations  for  modern  warfare,  but  assumptions  4  and  5  are 
replaced  by  [ Ref .  3]: 

4.  Each  firing  unit  is  aware  only  of  the  general  area  in 
which  enemy  forces  are  located  and  fires  into  this  area 
without  feedback  about  the  consequences  of  its  fire; 

5.  Fire  from  surviving  units  is  uniformly  distributed 
over  the  area  in  which  enemy  forces  are  located; 

6.  Each  unit  presents  the  same  vulnerable  area  to  enemy 
fire.   This  vulnerable  area  is  much  larger  than  the  effec- 
tive area  of  a  single  round  of  enemy  fire.   Additionally, 
the  number  of  hits  required  for  a  kill  obeys  a  geometric 
probability  law. 

The  state  equation  that  relates  the  force  levels  and 
does  not  explicitly  contain  time  for  the  aimed-fire  assump- 
tions is  known  as  Lanchester 's  Square  Law, 

b(X/  -  X2)   =  a(j;2    -  y2)  . 
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The   state   equation    for  the  area   fire   model   yields 
Lanchester's   Linear.  Law, 

b  (x^    -   X)    =  a(yo    -    y)  . 

From   these   state   equations   many   important   results    can   be 
deduced,    for   example,    battle   outcome  prediction  conditions. 
Further    analysis   of  Lanchester's  classical   equations   and 
their   implications   will  not   be   done   here.      Suffice  it   to  say 
that   Lanchester*s   original   work   forms   the  basis   for  the  ana- 
lytical  combat    model,    with   his    simple   model   serving   as  a 
point   of  departure   for  the   more   complicated   and   elaborate 
analytical   models  discussed   later  in   this  thesis.      Readers 
of   this   thesis    and   students   who   intend   to  continue   with   work 
initiated   here   should   have  a   sound   understanding   of   Lanches- 
ter's   basic  models   before   continuing   further. 

Lanchester*s  equations   for    modern   warfare  and  his   equa- 
tions  for   area    fire   have,    over   the   years,    been  thoroughly 
dissected   and   analyzed  by   students   of   differential   equations 
and   students   of   analytical   combat  models.      One  such   student 
was   Sordon   Clark,    a   doctoral   student   at   Ohio    State  Univer- 
sity  in   the  late    1960*s.      He   was   well   aware    of  the  power 
that   simulation    models   such   as    Carmonette,    a    high   resolution 
combat    model,    had   for   describing   the   complex    process    of  com- 
bat.     However,    the   major   drawback  to   simulation,    as   he 
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perceived  it,  was  that  the  method  was  slow  and  expensive. 
.Clark  proposed  that  the  use  of  an  additional  model  to  inter- 
pret simulation  results  would  greatly  benefit  Army  studies 
which  required  numerous  runs  with  the  simulation  model.   He 
went  on  to  develop  this  analytical  model  which  he  called  the 
Combat  Analysis  Model  (COMAN) .   This  much  simpler  and  faster 
running  model  would  greatly  reduce  computer  execution  time 
and  costs,  and  could  serve  as  a  filter  model  for  Carmonette. 
A  similar  approach  to  ClarJc^s  methodology  in  developing 
COMAN  will  be  used  here  to  develop  the  analytical  model  for 
STAR.   The  next  chapter  details  the  first  steps  taken  to 
build  the  analytical  combat  model  for  STAR. 
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III.   DAaSTAC  MODEL 

The  aggregated  force-on- force  analytical  model  is  an 
integral  part  of  combat  modelling  and  serves  several 
significant  purposes.   Hith  the  high  cost  of  computer  time 
and  the  limited  computer  core  storage  available  to  most  mil- 
itary users,  the  analytical  combat  model  offers  a  viable 
alternative  to  detailed  high  resolution  simulations. 

Analytical  combat  models  can  be  used  in  a  stand  alone 
configuration  or  in  conjunction  with  a  high  resolution  simu- 
lation.  In  the  stand  alone  mode,  the  models  are  usually 
"transparent",  or,  in  other  words,  easy  to  understand.   They 
are  also  small,  guick  but  most  importantly,  they  increase 
the  user's  understanding  of  the  process  of  combat.   In  con- 
junction with  high  resolution  combat  simulations,  the  ana- 
lytical model  offers,  as  Gordon  Clark  states  in  his  doctoral 
thesis,  the  opportunity  to  "extrapolate  the  results  of  the 
combat  simulation  to  predict  the  outcome  from  force  mixes 
that  were  not  explicitly  simulated."   In  this  way  it  can  act 
as  a  filter  to  eliminate  certain  force  mixes  from  further 
investigation.   Also,  Clark  identifies  another  use  for  it  as 
"an  integral  component  of  a  large  unit  model  to  predict  the 
outcome  from  individual  unit  actions."  [Eef.  U] 
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A  major  objective  of  this  thesis  was  to  build  a  simple 
aggregated  combat  model  on  the  Apple  II  microcomputer.   The 
input  parameters  for  the  model  were  estimated  using  the  out- 
put data  from  a  STAR  run.   This  concept  of  estimating  the 
attrition  rate  coefficients  involved  additional  wcrlc  in 
developing  a  data  postprocessor  compatible  with  STAR,  and 
significant  effort  in  the  statistical  analysis  of  the  infor- 
mation provided  by  the  postprocessor. 

There  are  several  reasons  why  the  Apple  II  computer  was 
chosen  for  the  development  of  the  model  rather  than  a  main 
frame  computer.   First,  the  need  exists  for  a  'portable' 
model  that  is  easily  accessible  to  the  Army  analyst.   A  desk 
top  type  system  that  can  run  a  model  quickly  and  with  mini- 
mum setup  time  will  be  a  valuable  asset  to  the  analyst. 
With  the  Apple  II  computer  gaining  wider  acceptance  in  the 
Army  for  various  uses,  the  accessibility  of  a  model  on  such 
a  system  is  greatly  increased. 

Second,  the  nature  of  the  model  is  such  that  it  lends 
itself  to  a  small  system  format.   The  aggregated  form  of  the 
model  reduces  the  required  core  storage  from  that  required 
by  a  high  resolution  model.   The  concept  of  OAJHSTAC  is  to 
serve  as  a  filter  model  to  weed  out  unacceptable  force  mizes 
or  scenarios  before  turning  to  a  high  resolution  model.   The 
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idea  here  is  to  reduce  the  time  and  cost  of  running  numerous 
iterations  of  a  large,  high  resolution  model  on  a  main  frame 
computer. 

Third,  the  cost  of  a  computer  system  to  support  the 
model  is  certainly  economical.   At  a  cost  of  under  $3000, 
the  Apple  II  computer  system  with  a  model  such  as  DAMSTAC 
provides  a  powerful  analytical  tool  for  a  minimal  cost,  cer- 
tainly in  relation  to  the  cost  of  a  main  frame  computer. 

Fourth,  the  Apple  II  has  the  capability  to  provide  addi- 
tional analysis  functions  on  the  output  of  the  model.   With 
such  programs  as  Apple  Plot  and  available  commercial  statis- 
tical and  numerical  analysis  packages,  most  of  the  needed 
analytical  worlc  can  be  accomplished  on  the  Apple  II  compu- 
ter.  This  can  result  in  additional  savings  in  time  by 
eliminating  the  requirement  to  transfer  any  data  from  one 
system  to  another. 

The  laurels  cast  on  the  Apple  II  should  not  go  without 
some  discussion  of  its  limitations.   The  Apple  II  is  limited 
in  capacity  to  6U  kilobytes  of  Random  Access  Memory  storage 
space.   So,  essentially  the  model  must  be  able  to  work 
within  the  amount  of  core  storage.  There  are  methods  of 
working  around  this  problem,  if  it  arises,  but.  this  will 
usually  create  a  new  problem  or  limitation.   For  example. 
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some  of  the  data  or  a  part  of  the  program  can  be  kept  on  the 
storage  diskette  and  put  into  coniputer  memory  only  when  it 
is  needed.   This  saves  space  in  the  computer's  memory  but 
increases  the  run  time  of  the  model  due  to  the  time  required 
to  access  the  diskette  to  retrieve  the  required  information. 
This  raises  another  limitation  of  the  Apple  II — its  run 
time.   Compared  to  a  main  frame  computer,  the  Apple  II  is 
very  slow.   This  is,  however,  relatively  speaking  since  com- 
puter time  is  normally  measured  in  milliseconds  and  nanose- 
conds.  For  the  analyst  waiting  for  the  output  from  DAMSTAC, 
the  model  on  the  Apple  II  would  take  about  four  to  ten 
minutes  to  run.   This  is  probably  acceptable  for  most 
analysts, 

A.   SCENARIO 

The  model  portrays  a  battle  ijetween  a  RED  tank  battalion 
attacking  a  BLOE  tank  company  in  the  10  x  10  KM  Pulda  Gap 
(Keyhole)  terrain.   See  Figure  2  for  the  terrain  map. 

The  BLUE  defensive  positions  and  the  RED  positions  along 
with  their  attack  routes  {coordinating  points)  are  provided 
in  Table  1.   The  RED  units  attack  along  three  preplanned 
attack  routes  toward  the  BLUE  defensive  positions.   Grid 
coordinates  for  unit  locations  and  coordinating  points  were 
computed  as  the  centroid  of  the  elements  of  the  unit.   In 
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Figure    2:      Terrain   map 
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star's  data  base,  each  tank  has  its  own  individual  location. 
For  example,  BLUE  unit  1  has  four  different  locations  asso- 
ciated with  the  it,  each  one  corresponding  to  a  tank  in  the 
unit.   The  centroid  of  the  unit  was  taken  to  be  the  average 
of  these  grid  coordinates.   Thus,  the  initial  location  for 
BLUE  unit  1  in  Table  1  is  the  average  of  these  coordinates: 
54150  100690         54150  100620 
54190  100520         54160  100790 

The  rate  of  RED  advance  is  a  constant  but  specified  by  the 
user.  BLOE  units  can  move  to  alternate  positions  at  a  rate 
designated  by  the  user.   However,  for  this  scenario,  BLUE 
units  remain  stationary  for  the  entire  battle. 

To  obtain  the  raw  output  data  that  would  correctly 
reflect  the  above  scenario,  it  was  necessary  to  "turn  off" 
many  of  the  submodels  of  STAR.   For  example,  the  Field 
Artillery  module  was  shut  off  by  scheduling  the  first  fire 
mission  well  after  the  battle  ended.   Weapon  systems  other 
than  tanks  within  the  fighting  units  were  rendered  useless 
by  setting  their  ammunition  counts  (basic  loads)  to  zero. 
Additionally,  the  communications  and  the  Air  Defense  modules 
were  not  played.   After  all  of  these  preliminary  actions 
were  taken,  the  result  was  a  battle  exclusively  between 
tanks.   This  effort  worked  well  to  reduce  a  complicated, 
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heterogeneous  force  battle  to  a  simple  homogeneous  force 
battle.   One  minor  problem  arose,  however,  when  a  weapon 
system  other  than  a  tank  was  detected  and  fired  on  when,  by 
the  scenario,  it  did  not  even  exist.   The  effects  of  this 
target  acquisition  and  engagement  were  considered  insignifi- 
cant and  were  disregarded. 

Table  1:   UNIT  LOCATIONS/COORDINATING  POINTS 

IzCOOED  SIZE  REMARKS 

100655  4  initial  location 

98120  4  initial  location 

97893  4  initial  location 

IzQOQRD  SIZE         REMARKS 

101491      10     initial  location 
96000  coordinating  pt 

101696      10     initial  location 
96000  coordinating  pt 

101491      10     initial  location 
96000  coordinating  px 


BLUE 

X-COORD 

ONIT 

1 

54163 

UNIT 

2 

54788 

UNIT 

3 

57448 

RED 

X-COORD 

ONIT 

1 

58727 
52000 

UNIT 

2 

581  10 
52000 

UNIT 

3 

58727 
52000 

B.   ASSUMPTIONS 

The  attrition  rate  functional  form  chosen  for  the 

DAMSTAC  model,  which  will  be  discussed  later  in  this  thesis, 

requires  the  following  assumptions. 

1.   The  conditional  kill  rates,  a  and  /3  ,  are  constant  and 

independent  of  time  if  an  individual  firer  has  at  least 
ona  acquired  target.   If  there  are  no  acquired  targets, 
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the   Jcill   rates  are   zero. 

2.  P   and  q,    the   probabilities  of   nondetection,    are 
assumed  to   be   constant    functions  of   time,    that   they   are 
independent   of  the   force  size   auid  the   number  of   acquired 
targets,    and   that  the   probabilities   p   and   q  are  equal  for 
all    x-firer/y-target   and  y-f irer/x-target   combinations 
respectively. 

3.  Synergistic  effects   are  neglected,    therefore,    the 
actual   force   attrition   rate  is  equal   to  the  sum  of   the 
individual   kill   rates   for  the   opposing   forces. 

Additionally,    the  assumptions   listed   below   are  designed  to 
remove    the   distractions   of   a   complicated   scenario   and   to 
make   the   initial  model   easier   to  understand: 

1.  All   units   have   point    locations;    ie.    a   battalion    is 
represented   by   a   grid   coordinate   and   all   calculations   for 
the    battalion   are   made    from   that    location. 

2.  Units   are   homogeneous.    That  is,    they  are   composed  of 
identical   firing   elements/weapons  systems. 

3.  Fire   is   uniformly  distributed  over   surviving  enemy  tar- 
gets   within    range. 

U.    Sufficient   supplies   of   ammunition    for    both   sides   exist 

for   the   duration   of   the   battle. 

5.    Units   move   at   a   constant   rate   of   speed. 
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C.   SYNOPSIS  OF  THE  COMPUTER  MODEL 

The  DAMSTAC  computer  program  is  a  modular  type  program 
written  in  the  Applesoft  BASIC  language.   The  program 
accesses  external  data  files  as  a  source  of  input  of  initial 
data  for  the  two  forces.   These  data  files  are  produced  by 
the  user  with  the  help  of  the  Force  Maker  utility  program 
described  in  Appendix  C.   Additionally,  the  program  has  the 
capability  to  make  output  files  that  can  be  stored  on  a 
diskette  and  accessed  by  the  user  with  the  Results  Reader 
utility  program,  also  described  in  Appendix  C.   Results 
Reader  will  print  the  output  to  the  screen  or  to  a  printer. 

A  copy  of  the  program  for  the  model  is  included  as  List- 
ing A  in  Appendix  F.   A  sample  run  of  the  model  including 
user  inputs  is  included  as  Figure  3.   Flow  charts  of  the 
main  program  and  subroutines  are  included  as  Figure  U,  Fig- 
ure 6,  and  Figure  7. 

''  •   Initialization  of  Arrays  and  Variables 

Lines  30  thru  170  dimension  the  arrays  and  matrices 
and  initialize  hardwired  data.   Dimensions  are  set  to  handle 
up  to  five  units  in  each  force  and  up  to  ten  coordinating 
points,  including  the  initial  location,  for  each  unit.   The 
user  can  modify  these  capacities  by  changing  the  dimensions 
of  the  appropriate  variables. 
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Lines  180  thru  350  read  BLUE  force  data  from  a  file 
to  include  the  number  of  units,  initial  locations,  any  coor- 
dinating points,  and  rates  of  advance,  and  compute  total 
force  level  based  on  user  input  of  the  unit  force  levels. 
The  time  increment  (DT)  is  set  to  30  seconds  while  the  maxi- 
mum time  of  the  battle  (MT)  is  set  to  1500  seconds.   These 
values  are,  of  course,  adjustable  by  changing  them  in  line 
80.   For  both  forces,  the  breakpoint  criteria  is  set  at  0.3 
of  the  force  level  and  movement  criteria  is  set  at  0.5  of 
the  unit  strength,  even  though  the  movement  criteria  was  not 
used  for  any  runs  of  the  model. 

Lines  360  thru  520  perform  the  same  operations  for 
the  P.ED  forc9s  as  lines  180  thru  350  perform  for  the  BLUE 
forces. 

Lines  530  and  600  ask  if  the  printer  should  be 
turned  on  or  off  and  calls  the  printout  subroutine  to  print 
the  initial  data. 

2.   Maia  Program  (lines  6IO-78Q) 

The  main  program  performs  these  operations: 
increments  the  time, 

checks  if  the  time  exceeds  the  max  time, 
calls  the  movement  and  attrition  subroutines, 
recomputes  current  total  force  levels  after  attrition,  and 
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Figura   H:      Flow    charii    of  the   Main   Program 
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checks  force  levels  against  the  designated  breakpoints  for 
each  force,  to  determine  if  the  battle  has  ended. 
If  the  breakpoints  or  the  max  time  are  exceeded,  the  battle 
terminates  and  the  program  prints  the  force  levels.   Addi- 
tionally the  program  will  write  the  force  levels  for  each 
time  interval  onto  the  diskette,  (lines  1470  thru  1580) 
3.   Movement  Subroutine  (lines  790-1  160) 

Both  RED  and  BLOB  forces  are  capable  of  moving. 
Movement  of  the  BLUE  force  is  nullified  by  setting  its  move- 
ment rate  to  zero  in  the  initialization  part  of  the  program. 
The  movement  subroutine  uses  basic  geometry  to  compute  new 
coordinates  for  the  moving  unit  (see  Figure  5) .   The  move- 
ment process  is  based  en  moving  a  unit  from  one  coordinating 
point  to  the  next  in  a  straight  line.   The  amount  of  move- 
ment depends  on  the  rate  of  movement  (BH  (I)  or  RH  (I)  )  and 

»€XT 
CooAD.  PT. 


CUtfRtHT 
iiMITLOC 


Figure   5:      Geometric   Representation    of   Unit    i^ovemeni 
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Figure   6:      Flow   Ciiarr   of    the    Movement    Subroutine 
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the  time   interval   length    (DT) .      If   the   distance   computed   is 
farther   than   the  next    coordinating   point,   only  movement  to 
the   coordinating   point   is    accomplished.      Lines  800  to    970 
perform   computations   to   move   RED  forces    while  Lines   980  to 
1150   perform   the  same   computations   for    BLUE    forces. 

'♦•       Attrition   Subroutine    (lines    1120^480) 

A  range  check   is   made   for  each   firer/target   combina- 
tion  to   insure   the   units   are   within   firing   range.      If   so, 
attrition  caused   by  that   firer    on  the   target    is   computed. 
Next  total  attrition   for   each   unit  is  computed  based   on  the 
partial   attritions  of   the   firers  as  a   fraction  of  their 
force   levels.      In   other   words,    when   a    firer    is  firing   at   3 
enemy    units,    each   target    unit   receives    1/3   of  the   firer»s 
fire   power. 

-••       Printout   Subroutine    (lines    1490-1530    &  lines 
1760-1940) 

There   are   two    different   output   formats   used   in   the 
program.      The    first    (lines    1490    to    1530)    prints   only    a  sum- 
mary of    the  total   force    levels    for   both    forces   for  each  time 
interval.      This   is   the  same   data  that  is   saved  to   a   file   by 
the   final   portion   of   the   program.      The   second   format    (lines 
1760   to    1940)    prints    cut    the    units   of   each   force    with   their 
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curren-t  locations  and  strengths.   This  format  allows  the 
user  to  see  which  unit,  specifically,  is  being  attritted  and 
the  grid  coordinates  of  that  unit  at  the  time  of  attrition. 
6.   Output  File  Making  Subroutine  (lines  154Q-1750) 
This  subroutine  allows  the  user  to  make  a  file  of 
the  output  from  the  current  run  of  the  program.   This 
routine  actually  produces  two  files,  one  consisting  of  the 
attrition  rate  coefficients  and  probabilities  of  non-acqui- 
sition, and  the  second  containing  the  time  and  force  levels 
of  both  forces  at  that  time.   These  files  are  seperated  so 
that  the  second  file  can  be  used  with  the  APPLEPLOT  program 
to  produce  a  graph  of  the  two  force  levels  as  a  function  of 
time.   This  file  has  the  postscript,  ".RESULTS",  added  to 
its  file  name,  while  the  first  file  with  the  attrition  rate 
coefficients  has  the  postscript  ".COEFS". 

D.   METHODOLOGY  FOR  COMPaTING  PARAMETER  ESTIMATES 

The  process  of  building  a  computer  simulation  combat 
model  requires  the  development  of  methodologies  for 
modelling  the  various  aspects  of  combat.   Constructing  a 
f orce-on-f orce  attrition  model  requires  methods  for  consoli- 
dating combat  phenomena  into  simplified  algorithms  with 
coefficients  representing  one  or  more  facets  of  combar. 
This  section  discusses  one  significant  part  of  the 
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forc9-on-f orce  attrition  model:  the  methodology  for  estimat- 
ing the  attrition  rate  coefficients  in  the  Lanchester  type 
attrition  rate  equations. 

In  deriving  the  equations  for  the  attrition  rates  (dx/dt 
and  dy/dt)  ,  it  is  necessary  to  understand  the  processes  that 
must  occur  to  cause  a  firer  to  inflict  attrition.  These  key 
processes  can  be  summarized  as: 

the  environmental  process; 

the  target  acquisition  process; 

the  target  engagement  process. 
The  environmental  process  involves  the  terrain  and  weather 
factors  that  affect  f irer/target  line  of  sight  (LOS) .   The 
target  acquisition  process  represents  the  visible  capability 
of  the  firer  to  detect  a  target  on  the  battlefield,  while 
the  target  engagement  process  describes  zhe    firer 's  ability 
to  hit  and  kill  an  acquired  target.   These  three  processes 
will  all  be  represented  by  the  attrition  rate  coefficient, 
A,  in  the  attrition  rate  equation.   This  coefficient  will  be 
broken  down  into  its  three  components  which  describe  mathe- 
matically the  phenomena  of  the  three  processes.   The  three 
components  are: 

P^  =  the  probability  that  a  target  is  visible  to  the  firer 

(environmental  process) 
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P_  =  thG  probability  that  a  visible  target  is  acqaired 

(target  acquisition  process) 

q;=  the  casualty  rate  at  which  a  firer  attrits  the 

acquired  target  (target  engagement  process) 

a  should  not  be  confused  with  A,  the  attrition  rate 
coefficient,  of  which  a  is  a  factor.   The  attrition  rate 
equation  may  be  represented  in  words  by: 

dx/dt  =  -  [P(tgt  is  visible)] [P(tgt  is  acquired)]  [casualty  rate] 
The  second  probability  above  can  be  broken  down  further  to 
P(tgt  is  acquired  |  tgt  is  visible).   Note  that  this  is  a 
conditional  probability  since  it  is  assumed  that  the  target 
must  be  visible  to  be  acquired.   By  identifying  the 
mathematical  expressions  for  each  of  these  components  of  the 
attrition  rate  coefficient,  the  attrition  rate  equation  can 
be  obtained. 

1.   The  Environmental  Process 

There  are  generally  two  ways  to  represent  line  of 
sight  in  a  model: 

as  a  mathematical  simulation  of  actual  terrain; 

as  a  stochastic  process. 
The  mathematical  simulation  method  normally  involves  a 
digital  representation  of  terrain  elevation  over  a  segment 
of  the  terrain.   The  elevations  of  the  segments  between  the 
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firer  and  the  target  are  used  to  determine  if  line  of  sight 
exists.   Virtually  all  high  resolution  models  use  this 
method,  or  variations  of  it,  to  determine  line  of  sight. 
STAR  uses  a  different  form  of  mathematical  simulation  to 
represent  terrain.   This  is  called  functional  terrain,  in 
which  hills  and  forests  are  represented  by  elliptical  con- 
tours of  elevation.   Hills  are  modelled  by  describing  them 
in  terms  of  bivariate  normal  equations  with  varying  parame- 
ters for  each  hill.   Forests  or  built-up  areas  are  repre- 
sented by  an  elliptical  pattern  on  the  ground  again  with 
varying  parameters  for  each  forest. 

For  a  low  resolution  or  highly  aggregated  model,  the 
mathematical  simulation  method  for  line  of  sight  is  usually 
not  practical  or  desireable  in  terms  of  core  space  and  com- 
putation time  required.   Therefore,  the  second  method,  in 
which  LOS  is  represented  as  a  stochastic  process,  is  used. 
The  concept  of  this  process  is  as  follows.   The  effect  of 
terrain  on  the  process  of  line  of  sight  is  that  a  target 
falls  into  one  of  two  categories  or  states:   either  it  is 
visible  or  it  is  invisible. 

Assuming  a  steady-state  situation  for  the  terrain, 
probabilities  can  be  assigned  to  each  of  these  two  states. 
These  probabilities  represent  the  percentage  of  time  that 
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the  target  will  be  in  that  state  (visible  or  invisible)  for 
that  terrain  over  a  reasonable  amount  of  time.   Graphically, 
the  states  of  the  target  can  be  represented  by  the  diagram 
in  Figure  8. 


TAR&eT 

rr 

Tl 

t: 

v»stsi.ff' 

r 

n 

1        1 
1 

1       ' 

1 

TAUGffT 
Z>IV  IS  181.6 

1 

1 

'^    TIME 


T? 


Figure  3:      Two-stats  Stochastic  Graph   of   Intervisibility 

Since   this   is   a   stochastic   process,    there   is    an 
exponential   distribution    associated   with   each  of    the   zvo 
states.      This   distribution  represents   the  time  between  two 
occurrences  of   the  same  state.      On   the   graph   in   Figure  8, 
the   ^ime   between   the    and    ox   T^    and   the   start   of   T^    is    expo- 
nentially  distributed    with   mean,  i  .      LiJcewise,   rhe  time  bet- 
ween  the   end    of    T,     and  the   start  of   tJ  is   exponentially 
distributed  with   mean,  ^  ,      Graphically,    this   can    be   repre- 
sented   by   Figure   9,    where   the   circles    represent   the   two 
states,    visibility   and   invisibility,    and   the    arrows 
represent   the   transition   from   one   state    zo   the  other    at   the 
given    rate. 
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Figure   9:      Two-state  Stochaistic  Model   of  Intervisibility 

The  attrition   caused   by   a  firer   on   a   target   can 
occur  only   when   the  targe::   is   visible    (assuming   no  artrition 
froi  wild,    blind  shots  or   indirect  fire)  .      Therefore,    the 
attrition  rate    is   a   function  of   the   probability   that   the 
target    is   visible  to   the   firer.      This   probability   will  be 
deno-ed   as   P    ,    and  can  be   expressed  as   the  proportion   of 
time  that   the   target   is  in  the   visible   state   in   a  steady 
state   situation,    p     =     -  ^  ,,     . 

2.      The  Target   Acquisition   Process 

Assuaing   now   that   line   of  sight    exists  between  a 
firer   and  a  target,   the  process   of   target  acquisition   must 
now   be   represented.      There  are    two    primary  modes   of   target 
acquisition: 

Serial  acquisition,    and 

Parallel   acquisition. 
In   the    serial    acquisition   process,    targets    are  acquired 
between    engagements.      That    is,    in  concept,    the  firer    is 
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concant rating   on  the    engaged  target    and   does   not    or   is  not 
able  to   acquire   other   targets   on  the   battlefield.      Any  addi- 
tional  acquisitions   made   prior   to  engaging   the  current  tar- 
get  are   lost    and  must    be   re-acquired  after   the  current 
engagement    has   ended.      Examples   of   firers   who   use  serial 
acquisition  are: 

any   firer  who   is   "buttoned   up"  and   has  his    field  of   view 

narrowed   to   a   thin   band; 

wire   guided   missiles/anti-tank  systems. 

In   the   parallel   acquisition   mode,   the   firer   main- 
tains  his   list   of   acquired  targets   throughout   an   engagement 
(except   those   lost   through   loss   of   line    of   sight) .      Thus, 
when   an    engagement    has   ended,    by  the   target    being   killed  or 
loss  of   line   of    sight,    the   firer   may   immediately    engage  a 
target    that   was    previously   acquired   and    is   still    a  valid 
target.      All   systems    which   do   not   use   serial    acquisition   are 
assumed   to    use   the   parallel   mode. 

For   purposes   of  this   thesis,    all   firers   are   assumed 
to   use    parallel    acquisition.      Therefore,    serial   acquisition 
will   not   be   discussed   any    further.      The    reader  should   note 
that   the   serial   acquisition   process    does   have   validity   in 
modelling   action   on   the   battlefield,    and   should   be  consid- 
ered  in    future   enhancements  of    this   model. 
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TO  explain  the  derivation  or  the  target  acquisition 
portion  of  the  attrition  rate  equation,  the  diagram  in  Fig- 
ure  9   must    be   modified  to   that   in  Figure    10. 
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Figure    10:      Three-state  Stochastic  Model  of   Inter visibility 

The  two-state   stochastic  process   has    been  changed   ro 
a  three-state   process    by   dividing  the   visible  state  into   the 
visible-not  acquired  srate    (YNA)    and  the   visible- acquired 
state    (VA) .      What   has    occurred   is  that   a   target  can   move 
from  the  visible-unacquired   state  to   the   visible   acquired 
state   at   a   rate    \  .      Since   both   of   these   states   comprise   the 
visibility   state   in  total,    the    rate    of    moving   from  either   of 
them   to    the   invisibility   state   remains    at  n .      Notice   that 
there   is   no   direct    transition   from   the   invisible   state  to 
the   visible-acquired   state.      This  is   only  a    matter  of 
representation,    since    it   is   assumed   that    ao    more    than   one 
event,    ie.    moving   to   a   new   state,    will   occur    during   a    time 
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increment.   Thus,  the  situation  described  can  still  be 
represented  by  Figure  10. 

Another  point  of  interest  in  this  diagram  is  the 
fact  that  there  is  no  transition  from  visible-acguired  to 
visible-not  acquired.   This  implies  the  assumption  that  once 
a  target  is  acquired,  it  remains  acquired  unless  it  is 
killed  or  line  of  sight  is  lost,  in  which  case  it  transfers 
to  the  invisible  state. 

To  determine  the  effect  of  the  target  acquisition 
process  on  the  the  attrition  rate,  the  probability  that  a 
target  is  acquired  must  be  determined.   From  Figure  10,  and 
using  the  same  concept  of  a  stochastic  process  as  that  used 
in  section  III-D-1,  the  equation  can  be  written 

P  (tgt  acquired)  =  — r— ^^ .        (1) 

This   should  more   properly   be   stated   as   the   probability  that 
a  target   is   acquired   given   that   it   is   visible,   since    this    is 
actually  a   conditional  probability   with   the   state   of    acqui- 
sition  being   dependent  on   the   target   being   visible. 

Therefore,    the  probability    that   a  target    is   visible 
and  acquired   in   the  steady  state  is: 


P„,    =    P  (tgt    vis   &   acq   ) 

=    P^  (tgt    acq    I  tgt    vis)    P,,(tgt   vis)  (2) 


"Vx+mJ       (^Tj+|ij 
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The  implications  of  this  equation  will  be  discussed  later  in 
Chapter  V. 

3.   Formulating  the  Attrition-.rate  Equation 

Osing  equation  (2)  in  the  attrition- rate  equation, 
the  attrition  rate  of  a  single  firer  on  a  single  target  can 
be  summarized  using  expected  values  as: 


dx/dt  =  -P^^  oi  (3) 


That  is,  the  attrition  rate  equals  the  rate  of  killing  an 
acquired  target  times  the  probability  that  a  target  is  visi- 
ble and  acquired.   However,  since  there  are  likely  to  be 
more  than  one  firer  and  more  than  one  target  on  the  battle- 
field, this  equation  musr  be  modified  to  reflect  this. 

If  a  firer  is  going  to  cause  attrition  on  his  tar- 
gets, he  must  be  able  to  acquire  one  or  more  of  these  tar- 
gets.  For  firer  y1  with  x  targets,  the  probability  that  y1 
acquires  one  or  more  targets  is  just  one  minus  the  probabil- 
ity that  he  acquires  no  targets.   Assuming  the  same 
probability  of  acquiring  any  of  the  targets,  the  probability 
of  acquiring  none  of  the  targets  is  (1  -  P„^ )^  .   Therefore, 

V  a. 


45 


The  attrition  rate  now  becomes 

dx/dt  =  -P^  a  (5) 

for  a  single  firer  against  x  targets.   Since  there  are  more 
than  one  firar  on  the  battlefield,  each  assumed  to  be  acting 
indapendent  in  this  process,  the  attrition  rate  in  equation 
(5)  must  be  multiplied  by  the  number  of  firers(y)  and  the 
equation  is  now  complete: 

dx/dt  =  -P^   Qt  y  .  (6) 

In   terms  of  the   probability  of    acquiring   a  target,   this 
equation  can   be   written  as: 

dx/dt  =   -(1   -    (1  -   P^^)""   )  ay  .         (7) 

This  is   the  same   equation   which   Clark   uses   in   his  COMAN 
model,    that   uses  results    from   a   high   resolution   model    to 
estimate  the    value   for  the   conditional   casualty  rate   and  the 
probability  that   a  target    is  not   acquired.       Henceforth,    (1    - 
P^J    will   be   referred   to  as   p   for  RED   firer/BLUE   target  com- 
binations   (or   q   for  a    BLUE   firer/RED   target    combination) , 
U.       The  Target    Engagement   Process    (Bonder   vs   Clark) 

The   target   engagement   process   can   be    modelled    in   its 
most   detailed   form,    as   in   the   high   resolution   model,    or   in 
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its  simplest  form,  as  in  the  aggregated,  analytical  model. 
The  high  resolution  model,  such  as  STAR,  includes  explicit 
representation  of  every  factor  that  is  involved  in  the  pro- 
cess of  aiming  at,  firing  at,  hitting,  and  killing  a  target. 
These  factors  include: 

weapon  system  characteristics; 

weapon/ammunition  type  characteristics; 

target  characteristics  such  as  size,  shape,  vulnerability, 

etc. 

In  contrast,  the  low  resolution,  aggregated  model 
represents  the  engagement  process  by  one  value — the  attri- 
tion rate  coefficient.   This  is  the  method  to  be  employed  in 
the  model  presented  here.   The  task  now  is  to  develop  a 
method  to  accurately  estimate  this  coefficient  for  use  in 
the  attrition  rate  equation.   Note  that  we  are  actually 
talking  about  two  attrition  rate  coefficients  and  two  equa- 
tions.  The  derivation  and  method  for  estimating  one  is  the 
same  as  the  other,  so  for  purposes  of  the  discussion,  refer- 
ences to  an  attrition  rate  coefficient  will  be  in  the 
singular. 

There  are  presently  two  basic  approaches  to  the 
problem  of  estimating  attrition  rate  coefficients  for  paral- 
lel acquisition  and  stochastic  line  of  sight.   The  first  is 
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the  use   of   the   independent   analytical   model,    involving 
weapon   system   characteristics  and  the   use  of   algorithms  to 
determine  the   expected  time   for   a  firer-  to   kill   a  single 
target.      The   attrition  coefficients  are   then   the   reciprocal 
of  these  expected  times.      This    method  will   be  referred  to  as 
Bonder's  method  after   Seth   Bonder  who  has  used  it  exten- 
sively   in   his   combat   models   such  as    VECTOR    II. 

The  second   method   is  the  fitted- parameter  analytical 
model,    involving  the   use   of   a   separate   high   resolution, 
Monte   Carlo  combat   simulation  to  produce   data  to   determine  a 
maximum   likelihood  estimate   of    the   attrition   coefficients. 
This  method   will   be   referred  to   as   Clark's   method,   since 
Gordon   Clark   popularized   it   in   his   doctoral   thesis  and  in 
the   development    of   the  Combat   Analysis   Model    (COMAN) . 

Most    of   this   section   involves   a    discussion  of 
Clark's   method   and  the  practical  application   of   it.      How- 
ever,   in  order   to   provide   a   contrast   to   that   method,    a   short 
description   of    Bonder's  method   will   be   provided. 

In   general,    the  concept   of    Bonder's   method  is    to 
produce   an   algorithm  representing  the   expected  time  that   it 
takes  a   firer  to   kill    a  single   target.      The   elements  of  this 
algorithm   include   computations    for   line   of  sight,    target 
acquisition,    target   selection,    and   the   killing  process. 
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Each  element   accounts   for   a   certain   amount  of  time  in   the 
expected  time  to   kill   a  target.      This   entire   process    is  seen 
by   Bonder  as   a    Markov   renewal   process    where   the   outcome  of   a 
round    fired  is   dependent   only   on  the   outcome    of   the   previous 
round.      By   using  weapon  systems  characteristics   data    as  raw 
input,    times   to    acquire,    select,   and   kill  a    target  can   be 
determined   and   used   in   the   algorithm   for   expected  time   to 
kill  a   target. 

Bonder's   method   uses   an   independent,    analytical 
model   requiring   a   large  data   base  of    factors,    such   as   range 
dependent   weapon  system   characteristics,    quantifications   of 
line  of   sight   and   target   acquisition   parameters,    to   support 
the  computation   of  the  attrition  rate   coefficient. 

A   detailed  explanation   of  Bonder's   method   for    esti- 
mating  the   attrition   rate   coefficients  can   be  found   in  the 
basic   documentation    for   the   VECTOR  II  model  [  Ref .    5].       Basi- 
cally,   Bonder   says   that   the   attrition   rate   coefficient   A . . is 
represented   in   the   parallel   acquisition    mode    by   A- ■   =    Q-  ■  af 
where   Q^^is  the    probability   that  a   given   group-i    weapon   is 
firing    at   a  group- j   target.      Q.  .    can   be    expressed   in    terms 
of   the    probability   that   a   group-j   target    is    visible   and  has 
been   detected   by   a   grcup-i   firer.      This    probabiltiy   is 
denoted    as   S-.  ,    and  the   expression    for   Q..    is   : 
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"fi' (1   -   S..)    for    j   >   2    . 
Assuming   homogeneous    forces,    this  equation   reduces  to 


Q.  .=   S 
ID         11 


Q..    =    S. .    ,    and  the   i-j    subscripts  can   be    eliminated 
i:  ID 

since  there  is  no  distinction  among  types  of  force  elements. 
His  resulting  expression  for  S  is  identical  to  that  of  equa- 
tion  H    ,    where    S- • and    P,    are   defined  the   same. 

Bonder's   expression   for   the   conditional   casualty 
rate,    a   ,    is   written   as  the   inverse   of   the   expected   time  to 
kill   an   acquired  target.      This    expression,    when   further  bro- 
ken  down   contains   elements  dependent   on    the    weapon  system 
characteristics,   target  characteristics,    and   a  number   of 
other   subsystem   performance   parameters.      These  parameters 
were  determined   empirically  outside  the   operational   environ- 
ment  of    the   battlefield  and   molded   together    in  this    pre- 
determined  structure   -co   produce   the    value  of    the    attrition 
coefficient.      This   derivation   of  the   attrition  coefficient 
is   the    point    where   Bonder's   method   and   Clark's   method 
differ. 

Clark's   methodology   is   based   on    the    assumption   that 
the   process  of   attrition   in    combat    is   a    nonhomogeneous,    two 
independent   type   Poisson   Process   where   the   time   between 
casualties   has   an    exponential   distribution    with    parameter  \  , 
called    the   total   conditional   casualty   rate.      Likewise,    the 
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time  between  RED  casualties  and  the  time  between  BLOE 
casualties  each  has  an  exponential  distribution  with  parame- 
ters x^  and  \,  respectively,  where  X     =  '^r"*"^b'  ^^°^   "the 
theory  of  Poisson  Processes,  the  probability  of  a  RED  casu- 
alty occurring  is^r/x  #  and  the  probability  of  a  BLUE  casu- 
alty occurring  is^b/^   .   With  this  basis,  Clark  is  able  to 
develop  a  likelihood  functon  in  terms  of  the  conditional 
casualty  rates  a    and  /3  ,  of  which  the  casualty  rates  are 
functions,  ie,,  X^=   f  (a)  and  \j-=  f  (/?)  •   [Ref.  6]   It  should 
be  noted  that,  in  order  to  maintain  consistency  of  terminol- 
ogy, the  equations  attributable  to  Clark •s  work  have  been 
adjusted  to  reflect  the  definitions  presented  in  this  the- 
sis. 

Clark  also  uses  the  concept  that  attrition  is  a 
function  of  target  acquisition,  which  depends  on  the  prob- 
ability of  a  firer  acquiring  a  target.   This  concept  pro- 
duces two  more  parameters  which  must  be  estimated:  p,  the 
probability  that  a  BLUE  target  is  not  acquired  by  a  RED 
firer;  q,  the  probability  that  a  RED  target  is  not  acquired 
by  a  BLUE  firer. 

For  the  functional  form  in  Clark's  COWAN  model,  the 
equations  can  be  written  as: 
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£=  -a[i-p'"]   with  m(o)  =  m^ 

^  (8) 

fr=  -^[1-  q"  ]   with  n(o)=  no 
where  the  parameters  are  defined  by: 

f^=  the  conditional  attrition  rate  function  for  BLDE 
forces 

f ,.  =  the  conditional  attrition  rate  function  for  RED  forces 
Qr=  conditional  kill  rate  for  a  RED  weapon  against  acquired 

BLUE  targets 
/3=  BLUE  weapon  conditional  Icill  rate  corresponding  to  a 

p  =  probability  that  a  BLOE  target  is  unacquired  by  an 
individual  RED  firer 

q  =  probability  that  a  RED  target  is  unacquired  by  an 
individual  BLUE  firer 

m  =  number  of  BLUE  firers  surviving 

n  =  number  of  RED  firers  surviving 

These  functional  forms  depend  directly  on  the  abil- 
ity of  firers  to  detect  targets.   To  see  this,  consider  just 
the  conditional  attrition  rate  function  for  the  BLUE  force, 
f ^j .   The  definition  of  p  implies  that  the  value  of  (1-?"*)  is 
the  probability  that  a  RED  firer  has  at  least  one  acquired 
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BLUE  target.      Therefore,    the  attrition   rate    for   the   SLOE 
force   is   the   kill   rate  for   acquired   targets   times  the    prob- 
ability  of  having  an   acquired   target   rimes   the  number   of 
firers. 

It   should   be   noted  that   Lanchester's    classical 
attrition   rate   functional   forms   are  special   cases  of   Clark's 
attrition  rate   functional   forms,   and,    under   appropriate  con- 
ditions,  Clark's   equations  reduce  to   Lanchester's   equations. 
When  p   and   q   are   zero,    or   very   nearly   zero. 


f.     =   -    a    (I-O"*)  n  =   -    Of  n  (9) 


and. 


f^    =    -  13  (1-0'')  m    =   -     /3  m  (10) 


which   are   Lanchester's  classical  aimed  fire    equations.      When 

individual   targets   become    increasingly  difficult    to   acquire, 

p  and   q   assume   values    close   to   one   and  it  can   be    proven  that 
in   the    limit: 
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f^    ^  '    a  (1-p)  an  (11) 


and. 


fr    =   -    /3    (1-q)a«a  (12) 


If  the   attrition  rata   coefficient  a  is   set   equal   to     a(1    - 
p)    and   the  coefficient  b   is   set  equal  to     /3  ( 1   -   q)  »    then   zhe 
two  equations  above  are  identical  to   Lanchester's   equations 
for  area   fire. 

In   his   doctoral  dissertation,   Clark    developed  -he 
likelihood   function   for   a   heterogeneous   model,   however   since 
the   OAHSTAC  nodel   assumes   hooogeneous   forces,    the  necessary 
equations  have   been  reduced  to    a  simpler   form  than  that 
derived   by  Clark.      By   differentiating  the  logarithm   of  the 
likelihood   function   with   respect  to    p,    the   following    aqua- 
tion  is    derived: 

where 
C,      -      0    if   casualty   k   is   RED 
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1  if  casualty  !c  is  BLUE 

n       -     number  of   BLUE  survivors  prior   to   Icth  casualty 

a,       =     number  of   RED  survivors  prior   to  Icth  casualty 

t,    -   t,    ,  -   time  between   casualty   Ic  &   casualty   Jc-1 
k    r     k-l 

Likewise,  by  differentiating  the  logarithm  of  the  likelihood 
function  with  respect  to  a  ,  and  by  solving  for  a  ,  a  second 
equation   in  terms  of   ?  and   S  is   produced: 

-1 


£(1   -   ^'^"^)n^..(t^   -   t^.     ) 


Differentiating  with  respect  to  q  and  (3   produces  equivalent 

A  A 

type  equations  which  will  allow  the  estimates  for  q  and  (3   to 
be  determined  in  the  same  manner  as  follows. 

Generally,  the  two  equations  (13)  and  (14)  can  now 
be  solved    simultaneously  to  determine  the  values  for  p 
and  a   .   anfortunately,  this  is  easier  said  than  done.   The 
complexity  of  these  equations  (more  so  in  the  heterogeneous 
case)  does  not  allow  explicit  equations  for  p  and  a:  .    The 
remainder  of  this  chapter  will  discuss-  a  method  for  deter- 
mining the  estimates  of  p  and  a   from  these  equations  and  a 
way  to  extract  the  necessary  data  from  a  high  resolution 
aont2  Carlo  simulation  run  to  compute  these  estimates. 
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E,   ESTIMATING  PARAMETERS  FROM  A  HIGH  RESOLUTION  MODEL 

As  mentioned  before,  equations  (13)  and  (14)  must  be 
solved  simultaneously  to  produce  values  for  a  and  p.   Gener- 
ally, these  equations  are  too  complicated  to  produce  such 
values  easily  so  other  techniques  must  be  used.   The  method 
suggested  by  Clark  is  to  enter  a  value  for  p  in  equation 
(14)  and  solve  for  o;  .   This  value  and  the  selected  value 
for  ?  are  then  put  in  equation  (13)  and  the  equation  is 
evaluated.   Notice  that  both  values  are  estimates  since  we 
are  computing  the  maximum  likelihood  estimates  of  p  and  a  . 
The  values  of  p  and  a    which  give  a  value  for  equation  (13) 
closest  to  zero  are  the  best  maximum  likelihood  estimates 
for  p  and  a   .   If  equation  (13)  equalled  zero  exactly,  the 
true  values  of  the  maximum  likelihood  estimates  would  be 
known. 

In  using  these  equations,  the  input  received  from  the 
high  resolution,  Monte  Carlo  simulation  must  include: 

the  time  between  casualties; 

the  number  of  RED  and  BLUE  survivors  prior  to  the  kth 

casualty; 

the  total  number  of  casualties  for  RED  and  BLUE  for  the 

battle  ; 

the  type  of  casualty  (RED  or  BLUE)  of  the  kth  casualty. 
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In  Clark's  theory  for  a  homogeneous  force  battle,  the 
conditional  attrition  rates  and  the  probabilities  that  a 
target  is  not  acguired  are  assumed  to  be  constant  over  the 
entire  battle.   This  means  that  equations  (13)  and  (14)  are 
summed  over  all  casualties  in  the  battle  and  the  estimates 
for  p  and  a  ,    which  are  determined  from  them,  are  also 
assumed  to  be  constant  over  the  entire  battle. 

In  applying  this  theory  to  STAR  output,  the  authors 
obtained  questionable  results  for  the  values  of  the  parame- 
ter estimates.   The  estimates,  p  and  q,  were  computed  to 
have  values  of  0.99  and  0.93  respectively.   These  values 
imply  that  both  BLUE  and  RED  forces  experienced  difficulty 
in  acquiring  targets.   Analysis  of  the  Fulda  Gap  terrain  and 
consideration  of  the  respective  locations  of  the  opposing 
forces  support  the  estimated  values  for  p  and  q.   However, 
the  stated  value  for  p,  0.99,  is  a  truncation  of  the  actual 
value  computed.   If  proper  rounding  procedures  had  been 
used,  p  would  have  been  set  equal  to  1.0,  which  means  that 
the  RED  forces  would  not  detect  any  BLUE  targets.   Nondetec- 
tion  implies  that  no  BLDE  attrition  will  occur.   But  this 
contradicts  rhe  STAR  battle  which  did,  in  fact,  produce  BLUE 
casualties  amounting  to  about  two-thirds  of  the  total  BLUE 
fores. 


The  estimates,  aand  (3,    were  computed  to  be  2.4  and  0.03, 
respectively.   The  relatively  large  value  of  a   implies  that 
when  the  RED  forces  are  able  to  detect  BLUE  targets,  the 
BLUE  force  will  suffer  heavy  casualties.   There  appears  to 
be  an  unusually  large  difference  between  the  values  for  a 
and  J  which  are  not  reflected  by  the  values  for  p  and  g* 
This  disparity  is  most  likely  attributable  to  the  restric- 
tion that  all  parameter  values  be  held  constant  over  the 
entire  battle.   Experimentation  with  different  parameter 
values  in  DAHSTAC  indicates  that  a  value  of  about  0.40  more 
closely  approximates  the  attrition  observed  in  the  STAR  bat- 
tle.  Experience  with  these  computations  reveals  that  this 
is  not  necessarily  the  case  and  that  a  more  likely  assump- 
tion is  that  the  parameters  are  constant  over  different  time 
intervals  of  the  battle.   In  general,  a  simple  battle  can  be 
divided  into  phases  where  the  intensity  of  the  battle  is 
constant  within  each  phase  but  different  among  the  phases. 
Basically,  the  phases  are: 

movement  to  contact  {low  intensity) 
initial  contact  (mid  intensity) 
main  battle  (high  intensity) 

withdrawal  from  contact  or  annihilation  (mid/low  inten- 
sity) 
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Moving  from  one  phase  of  the  battle  to  the  next  should  cause 
a  change  in  conditional  casualty  rates  (  aand/S  )  as  well  as 
the  probability  of  not  acquiring  a  target  (p  and  g) . 

The  question  now  becomes  how  to  determine  the  breakdown 
of  the  battle  by  phases  using  the  output  from  a  high  resolu- 
tion, Monte  Carlo  simulation.   That  is,  where  do  the  changes 
in  phases  occur  in  the  simulation?  One  possible  method  is 
to  look  at  the  number  of  casualties  as  a  function  of  time. 
Another  is  to  look  at  the  number  of  shots  fired  as  a  func- 
tion of  time  (see  Figure  11) .   In  both  cases,  time  is  seg- 
mented into  time  intervals  (eg.  every  100  seconds)  and  a 
histogram  of  the  casualties  or  shots  is  plotted  for  each  of 
these  intervals.   Having  done  this  for  several  battles,  the 
results  identify  fairly  clear  choices  for  the  end  of  one 
phase  of  the  battle  and  the  beginning  of  rhe  next  phase. 
However,  mor^  analysis  of  these  two  methods  reveals  that 
they  light  not  truly  identify  the  phases  of  the  battle  and 
that  m  alternate  method  should  be  considered,  primarily  one 
based  on  range.   This  concept  is  discussed  in  detail  in 
Chapter  VI. 

Using  the  number  of  casualties  per  time  interval,  a  dis- 
tortei  picture  of  the  intensity  of  the  battle  might  be 
obtained.   The  interaction  of  the  probability  of  acquiring  a 
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Figura    11:      Histogram  of   Shots   per  Time  Interval 

target   with  the   attrition   rate   of  a   firer  mighz   offset  each 
other  and  produce   results    which   conflict   with  the  rrue  bat- 
tle  intensity.      In   other   words,    if   the   acquisition  probabil- 
ity is   high  but   the  attrition  rate   is  low,    the  number   of 
casualties    in   a   time    interval   mighn   indicate    a  low  intensity 
phase   of   the   battle   when,    in   fact,    it    is    a   high   intensity 
phase.      Therefore,    this   method    should   be   used   with  caution, 
if   at   all.      The    method  of   using   the   number   of   shots   per  time 
interval   has    similar   faults.      Twenty    shots   by   one   firer   in   a 
time   interval   is   much    different   in   terms  of    intensity   than 
one   shot   by  twenty   firers.      Again,    this    method  should    be 
used    4izh   caution. 
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The  question  is  then:  what  method  should  be  used?   One 
method,  which  is  used  by  the  Combined  Arms  Combat  Develop- 
ment Activity  (CACDA/CASAA) ,  is  measuring  the  number  of 
shots  per  firer  per  time  interval, ^  Although  it  also  has 
some  faults,  it  does  provide  a  better  indication  of  the 
intensity  of  the  battle.   It  is,  for  the  most  part,  indepen- 
dent of  the  conditional  casualty  rates,  although  the  number 
of  shots  would  indicate  a  dependency  on  target  acquisition. 
This  dependency  is  reduced  by  averaging  the  number  of  shots 
over  the  number  of  firers.   In  effect,  this  measure  can  be 
considered  to  be  unbiased  by  the  changes  in  conditional 
casualty  rates  and  probabilities  of  not  acquiring  targets- 
Thus  one  is  able  to  get  a  feel  for  the  intensity  of  the  bat- 
tle from  the  histogram  of  the  number  of  shots  per  firer  per 
time  interval.   This  is  not  necessarily  scientific  reason- 
ing, but  it  seems  to  work. 

Once  the  phases  of  the  battle  are  delineated  for  a  high 
resolution,  Honte  Carlo  simulation  run,  the  next  step  is  to 
determine  how  to  treat  the  'dead*  time  from  the  beginning  of 
a  phase  to  the  first  casualty  and  from  the  time  of  the  last 
casualty  to  the  end  of  the  phase.   For  phases  in  the  middle 


^Conversation  with  Mr.  Ernest  Boehner  from  CASAA  on  5 
December  1981 
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of  tha  battle,  each  of  these  two  times  is  a  portion  of  the 
time  between  two  casualties  and  should  thus  have  a  bearing 
on  our  maximum  likelihood  estimates  from  equations (13)  and 

(14). 

To  resolve  this  situation,  the  time  of  the  last  casualty 
within  a  phase  is  used  as  the  end  of  the  phase.   For  exam- 
ple, if  a  phase  is  supposed  to  end  at  1200  seconds  and  casu- 
alties occur  at  1140  and  1225  seconds,  the  phase  will  be 
terminated  at  1140  and  the  new  phase  begun  at  1141.   This 
also  allows  this  time  between  casualties  to  be  accounted  for 
in  the  next  phase.   This  method  produces  some  variation  in 
phase  lengths  when  running  several  replications  of  the  bat- 
tle, but  this  variation  is  minimal,  especially  since  the 
beginning  and  end  of  phases  are  being  judged  from  something 
less  than  a  scientific  method. 

Several  replications,  as  many  as  affordable,  should  be 
used  to  compute  the  maximum  likelihood  estimates  for  the 
four  parameters.   By  using  a  large  sample  population  (the 
casualties) ,  a  more  accurate  estimate  of  the  parameter  can 
be  obtained.   Likewise,  when  the  heterogeneous  model  is 
used,  a  larger  sample  population  is  desirable,  since  the 
number  of  engagements  between  two  particular  weapons  systems 
will  be  much  less  than  the  number  of  engagements  in  the 
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homogeneous  case.   This  could  produce  values  that  are 
statistically  insignificant. 

A  postprocessor  to  the  high  resolution  model  was  devel- 
oped to  make  the  numerous  calculations  required  to  produce 
these  estimators.   The  postprocessor  uses  an  optimization 
routine  to  determine  the  maximum  likelihood  estimates  from 
the  given  equations.   Detailed  information  on  the  postpro- 
cessor is  provided  in  Appendix  E. 
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IV.       ALTERNATE    FUNCTIONAL    FORMS 
Different    military  situations  have   been    hypothesized   to 
yield    different    functional   forms   for   Lanchester-type    equa- 
tions.     The  term   "Lanchester-type  equations"    refers   to  any 
differential-equation    model   of   combat   and  the   term  "func- 
tional  form"    relates   to  the   structure   of    the    attrition  rate, 
i.e.,    how  the   attrition  rate   equation    is   written   and    which 
parameters   are   considered   explicitly   in   the    equation. 
Henceforth,    a   standard  shorthand  notation   will  be  adopted   to 
refer    to   the   functional   forms   used   in   homogeneous-force  Lan- 
chester-type  combat   models.      For  example,    Lanchester*s   equa- 
tions  of   modern    warfare  are   expressed  as    F|F    functional 
forms    (or,    simply   F|F   attrition)    because  they  are   propor- 
tional  to   only   the   number   of   enemy   firers.      Taylor  [ Ref .    7] 
provides  a   good   explanation  of   this   convenienr  shorthand 
notation. 

Generally,    the   classical   Lanchester    functional  forms   for 
expressing   the   attrition   rate   for  forces    using  aimed    fire   or 
area   fire   are   too   simple   to    be   used   "as    is"    in  a    detailed 
analytical   combat   model.      The  assumptions  required   for   these 
basic   equations   to   hold  are   too   limiting   and    do   not   allow 
for   all   of  the    processes   that  are  known    to   occur    in   combat 
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to    be   considered   in   the   force   attrition   model.      Some    of  the 
more   common   functional  forms   for  attrition   rates   that    have 
been  considered   and   used   in  Lanchester-combat  models    are 
shown    in   Table    2. 

Of    special    interest  is  the    T|T   functional   form.      This 
functional   form    is    particularly   useful   in  describing    attri- 
tion  in    target    rich   environments.      Peterson   [ Ref .    8]    hypoth- 
esized  that  the    equations, 

dx/dt   =   -ax    and   dy/dt   =  -by, 
characterize   the   early  stages   of  a   small   unit  engagement   in 
which   the   vulnerability   of  a   force   dominates   its    ability   to 
acquire   enemy   targets.      It   is   recommended  that  Peterson's 
T|T   functional    form  be  considered  in   the   heterogeneous  model 
described   in   Chapter   VI  to  model  attrition  in  the   initial 
phase   of  the   battle.      The   remaining   sections   of   this    chapter 
will   discuss   several   not-so-common   functional   forms   that 
have   been   formulated   and   used   in   "recently"    developed   combat 
models. 

A.       HELHBOLD-TYPE    EQUATIONS 

An    alternative   functional   form   for   homogeneous-force 
attrition   rates    was   proposed   by   R.    Helmbold.      Helmbold 
believed   that   when   opposing   force   sizes    were    grossly    une- 
qual,   the   larger  force   experienced    "inefficiencies  of   scale" 
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Table  2:   FUNCTIONAL  FORM  NOTATION 
FUNCTIONAL  FORM  DIFFERENTIAL  EQUATION 

F|F  dx/dt  =  -ay 

dy/dt  =  -bx 

FT! FT  dx/dt  =  -axy 

dy/dt  =  -bxy 

F| FT  dx/dt  =  -ay 

dy/dt  =  -bxy 

TIT  dx/dt  =  -ax 

dy/dt  =  -by 

(F+T)  I  (F+T)  dx/dt  =  -ay  -  ex 

dy/dt  =  -bx  *  dY 


in  producing  casualties.   Limitations  on  space  available  for 
maneuver,  reduced  reaction  time,  increased  command,  control, 
and  communication  problems,  etc.  result  in  the  reduction  of 
casualty  producing  effectiveness.   These  limitations  may 
well  prevent  the  larger  force  from  using  its  full  destruc- 
tive capability  against  the  smaller  force.   There  is  no 
definition  for  "grossly"  unequal  force  sizes,  but  the  fol- 
lowing hypothesis  is  proposed.   The  validity  of  using  Helm- 
bold'  s  equations  in  a  combat  model  is  directly  related  not 
only  to  the  force  ratios  (i.e.,  x/y  or  y/x)  but  to  the  rela- 
tive size  of  the  two  opposing  forces  with  respect  to  their 
surrounding  terrain.   For  example,  a  5  to  1  force  ratio  may 
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well  be   a   situation   where   Helmbold^s   equations  apply    if  both 
of  the   opposing   force   sizes   are  relatively   large    with   res- 
pect  to   the  terrain  available   to  them   for  movement,    conceal- 
ment,   etc.      Conversely,   a    force   with   a   30  to    1   force   ratio 
may   not   exhibit    inef f iciences   of  scale   if  the  opposing  force 
sizes   are  small   in   relative    sizes  and   with   respect  to   their 
surrounding  terrain.      Consequently   the   decision   to  use  Helm- 
bold»s    functional   form  for   attrition   rates   rather  than  some 
other   form    (like  Lanchester*s  equation   for   modern  combat) 
really   depends   on   the   tactical    "situation". 

Helmbold's    functional   form    introduces  a   modification 
that   accounts    for  the   reduced   fire   effectiveness    of  the  lar- 
ger force.      In    mathematical   form   his   equations  would    read: 
dx/dt   =  -a  g  (x/y)    y  with    x(0)     =  x^ 

dy/dt   =  -b   h(y/x)    x  with    y(0)    =   y^ 

where   a   and  b  represent  time-dependent   attrition-rate  coef- 
ficients,   and   g  (x/y)    and   h  (y/x)    represent  functions   that 
modify    fire   effectiveness   of   individual    x  and  y    firers 
respectively.      According   to   Helmbold*s   hypothesis   the   two 
functions   g(x/y)    and   h  (y/x)    must  satisfy: 

g(1)    =   h(1)    =    1 

g(u)    =   h(u)    =    E(u) 

E(u)  is  a  strictly  increasing  function  of  u. 
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So,    the    generalized  Helmbold-type  combat   equations  can   be 
written   as: 

dx/dt   =   -a    E(x/y)    y 
dy/dt   =   -b   E(y/x)    x 
For   the  case   in   which   E(u)    =   u      with   c  a   constant 
greater   than   or   equal   to   zero,    the   following   equations  for 
Helmbold-type   combat   are   obtained: 

dx/dt   -  -a    (x/y)''^  y        with    x(0)    =   x^ 
dy/dt   =  -b    (y/x)^"^  X        with    y(0)    =  y^ 
M  is   commonly   called   the   "Weiss   parameter"    with   W   =    1-c. 

The   above   equations  represent  the   special  exponential 
case   of    Helmbold   combat.    If   a(t)    and   b(t)    represent   constant 
attrition   rate   coefficients   denoted   by   a   and    b  respectively, 
then,    when    W   =    1,    the   above   equations   reduce    to   Lanchester's 
Equations   for   Modern    Combat.      When   H   =    1/2   Lanchester's 
Area-fire   Equations  are   obtained. 

B.       HOW    TO    CHOOSE    THE    RIGHT    FUNCTIONAL    FORM 

A   very  logical  question   that   readers   of   this    study  may 
ask  at    this  time   is,    "Now    that   I've   seen    some  of    the    Lan- 
chester-type   attrition   rate   functional   forms,    which  one 
should    I   choose    for   my  analytical  model?"      The  answer    to 
this   question    is   that   there    is    no  answer   to    this    question. 
The   choice  of   functional   form   is  purely    arbitrary  and   is 

68 


usually  left  to  the  discretion  of  the  model  builder.   The 
following  steps,  however,  may  help  in  this  choice. 

1.  Carefully  consider  the  scenario  that  the  combat  model 
is  intended  to  portray.   The  scenario  may  well  indicate 
that  some  particular  subset  of  the  functional  forms  is 
preferable. 

2.  List  the  assumptions  that  must  be  made  for  each  of  the 
functional  forms  in  the  subset  to  hold  and  see  which 
assumptions  are  acceptable  to  you  and  your  model  and  which 
are  not. 

3.  Observe  the  parameters  involved  in  each  functional 
form  and  eliminate  those  forms  for  which  the  parameter 
values  or  estimates  are  unobtainable  or  for  which  you  are 
not  willing  to  plug  in  arbitrary  values. 

At  this  poin*,  the  final  selection  of  functional  form  is 
purely  one  of  subjective  judgement. 

There  is  no  requirement,  however,  that  a  model  builder 
must  accept  one  of  the  already  developed  functional  forms. 
There  is  no  doubt  that  other  forms  can  be  formulated,  but 
the  formulation  will  require  extensive  research,  analysis 
and  work.   A  good  check  to  see  if  an  attrition  rate  equation 
is  reasonable  is  to  insure  that  the  Lanchester  classical 
equations  are  special  cases  of  the  new  functional  form. 
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Once   the    final   functional   form   has   been    selected   or 
formulated,    it    may  be   enhanced   as  necessary    to  reflect   more 
of   the   complexity   and    realism  of  combat.      Some  additional 
operational   factors   that    can   be   considered  are: 

Onit   breaJcpoints   and  battle  termination  criteria; 

Suppression ; 

C3; 

Reinforcements  or  replacements. 
The  main  assumption  that  allows  these  operational  factors  to 
be  "added"  to  the  attrition  rate  equation  is  that  synergis- 
tic effects  are  not  considered  in  Lanchester-type  force-on- 
force  combat  models.  No  synergistic  effects  imply  that  the 
Principle  of  Superposition  holds,  therefore,  the  effects  of 
the  operational  factors  are  independent  and  additive. 
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V.   FUTURE  ENHANCEMENTS 

There  are  numerous  enhancements  to  the  analytical  model 
that  can  be  made.   The  enhancements  that  are  recommended, 
for  the  most  part,  are  not  aimed  at  improving  the  inner 
workings  of  the  current  model,  but  at  building  successive 
models  that  are  more  detailed,  require  fewer  limiting 
assumptions,  and  are  generally  more  complicated. 

For  follow-on  models  that  retain  the  homogeneous  force 
concept,  it  is  recommended  that  the  model  have  the  ability 
to  allow  the  user  to  choose  alternative  functional  forms  for 
the  force  attrition  rates.   This  should  be  done  to  better 
describe  the  situation  of  the  battle  at  a  particular  time. 
A  likely  method  would  be  to  change  from  one  functional  form 
to  another  during  zhe   battle  in  order  to  more  realistically 
portray  the  attrition  occurring  then.   This  enhancement 
would  require  considerable  effort  in  developing  the  maximum 
likelihood  estimates  for  the  various  parameters  used  in  each 
alternative  functional  form,  additional  work  with  the  STAR 
postprocessor,  and  significant  change  to  the  current  DAMSTAC 
model.   Questions  which  arise  from  this  suggestion  are: 

Which  forra(s)  should  be  chosen? 
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How  do  you  know  when  to  change  functional  form? 

These  questions  are  partially  answered  in  Section  I7-C. 
Further  research  into  this  area  is  required. 

The  movement  subroutine  of  the  model  needs  to  be  modi- 
fied to  more  realistically  portray  the  movement  of  units  on 
the  battlefield.   Even  though  a  relatively  low-resolution 
time-step  operation  is  occurring,  the  resolution  of  the 
movement  of  units  can  be  refined  to  vary  the  rate  of  move- 
ment from  one  time  step  to  the  next,  if  necessary.   Cur- 
rently, a  change  in  rate  of  movement  occurs  only  at  the 
initial  contact  between  forces  and  at  breakpoints.   The  rate 
of  advance  of  the  RED  forces  should  be  allowed  to  vary  based 
on  decision  logic  which  considers  attrition  rate,  terrain, 
etc.   BLUE  forces  should  be  given  the  capability  to  move  ro 
alternate  or  fall  back  positions  when  the  batnle  situation 
dictates.   Decision  logic  must  be  enhanced  to  allow  BLUE 
forces  to  move  to  alternate  or  fall-back  positions.   It 
should  be  noted  that  this  "pace  of  battle"  question  has 
never  been  fully  explained  or  understood,  however,  the  move- 
ment routine  logic  in  DAMSTAC  should  parallel  that  of  STAR. 

The  concept  of  units  having  different  formations  is 
ignored  in  the  current  model.   Units  are  identified  by  a 


72 


point  location  only.   Some  work  needs  to  be  done  to  allow 
the  units  to  take  on  different  formations  for  road  marches, 
attacks,  etc.   In  order  to  do  this,  additional  attributes 
such  as  frontage,  depth,  orientation  and  geometric  configu- 
ration should  be  considered  to  describe  the  formation  of  the 
unit.  k   postprocessor  for  STAR  could  take  the  individual 
tank  locations  from  the  STAR  data  base  and  determine  the 
above  mentioned  attributes  for  input  into  the  DAMSTAC  model. 
Attrition  rates  on  both  sides  must  be  adjusted  for  the  vari- 
ous formations  of  an  attacking  unit. 

The  inclusion  of  indirect  fire,  suppression,  smoke, 
intelligence,  and  close  air  support  are  all  important  to 
making  the  model  a  usable  analytical  tool.   Additions  of 
communications  and  logistics,  primarily  ammunition  and 
petroleum,  oil,  and  lubricants  (POL) ,  are  also  important 
enhancements  that  must  be  included  in  the  model. 

As  was  seen  in  Section  III-E-3,  P^_^,the  probability  that 
a  target  is  visible  and  acguired  in  the  steady  state  is 


P   =   P   P 


(-xhr)   (-.f^  o^> 


Both  of  the  factors  represent  a  particular  element  of  the 
attrition  process,  one  the  terrain  affects  on  line  of  sight. 
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and  the   other,    the  acquisition    process.      It   should  be   appa- 
rent  that   these    two   factors   can   be   separated    from   each  other 
with  the  only   common   variable,  ^  ,.  the  rate  of  transfer  from 
the   visible  state  to   the   invisible   state.      With  this    factor- 
ization,    if  the   line  of  sight   parameters   for    a  particular 
terrain    type  are   known  and  the   attrition  rate  coefficients 
have   been   determined,    the   target  acquisition   capability  of   a 
weapon    system   against    another   weapon   system   can   be  deter- 
mined.     This    would  allow    the   target   acquisition   parameter   to 
be   used   for  other  terrain   types   to   compute   attrition   rate 
coefficients   without   having   to   estimate   the    coefficients 
from   a    high  resolution  simulation. 

The   concept   works   in  this   manner,      fiuns   are   made   on  the 
high  resolution   model    (STAR)    using   a   particular   terrain. 
The  attrition   coefficients  are   determined   using   the   maximum 
likelihood   estimates   described   in  Section  III-D.       With  pre- 
determined  data    for  the  terrain    coefficients,    the  line   of 
sight    parameter   can  be  factored   out   of    the   attrition   rate 
coefficient,    A.      This    would   leave  the   value    of  the  acquisi- 
tion  factor   which  should  theoretically   be  constant    for  the 
firer-target   type  combination   used   in   the  scenario.       With 
this   value   and    the   terrain   coefficient   for   other   types  of 
terrain,   the   attrition  rate   coefficient    can    be  determined 
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for  a  new  scenario  with  the  same  firer-target  combination, 
without  running  the  high  resolution  model  again.   This 
method  assumes  that  the  terrain  affects  and  the  target 
acquisition  process  are  independent  of  each  other.   The 
impact  of  this  concept  is  enormous  in  terms  of  time  and  cost 
saved  from  not  having  to  make  multiple  runs  of  the  high 
resolution  model. 

One  concept  which  is  not  strictly  an  enhancement  but  is 
an  important  consideration  is  the  idea  of  using  simulated 
combat  data  from  a  source  other  than  STAR.   Specifically, 
data  generated  at  the  National  Training  Center  from  actual 
field  exercises  could  conceivably  be  used  to  compute  attri- 
tion coefficients  in  much  the  same  manner  as  STAH.   What 
this  provides  is  not  only  an  alternate  source  of  data,  but 
also  a  different  model,  since  it  would  include  many  human 
factors  that  would  not  likely  be  accounted  for  in  STAR. 
Depending  on  the  intended  use  of  the  model,  one  source  or 
the  other  may  be  desirable. 

It  is  necessary  that  a  progression  from  homogeneous 
forces  to  heterogeneous  forces  be  made  if  DAMSTAC  is  ever  no 
be  really  useful.   This  ennancement  entails  a  significant 
amount  of  research,  model  building,  and  statistical  analysis 
in  order  to  implement  the  concept  of  heterogeneous  forces 
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into  the   analytical   model.      Target   priorities,   fire 
allocation   algorithms,    and   a  library   of   attrition  rate  coef- 
ficients  must    be  developed   for   various  scenarios    and    force 
mixes.      Changing  the   model  to  a   heterogeneous   type  requires 
some  changes   in   the  equations    (13)    and    (14),    and   thus   some 
significant  modifications   to   the  postprocessor  for  the 
model. 
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VI.   HETEROGENEOUS  MODEL 
DAMSTAC,  in  its  current  configuration,  models  combat 
between  two  opposing  homogeneous  forces,  but  actual  combat 
consists  of  many  different  weapon-system  types  operating 
together  as  combined  arms  teams.   For  example,  there  may  be 
infantry,  tanks,  antitank  weapons,  artillery,  etc.  on  each 
side.  k   necessary  extension  or  modification,  therefore,  to 
DAMSTAC  is  the  consideration  of  combat  between  heterogeneous 
forces.   This  chapter  briefly  indicates  how  the  basic  ideas 
and  methodologies  used  previously  to  model  combat  attrition 
for  homogeneous  forces  can  be  extended  and  adapted  to  the 
heterogeneous  force  structure. 

A  key  assumption  that  must  be  made  in  order  to  model 
combat  between  heterogeneous  forces,  where  Lanchester-type 
differential  equations  are  used,  is  that  the  attrition 
effects  of  various  different  enemy  weapon  systems  types 
against  a  particular  friendly  target  type  are  additive.   In 
other  words,  synergistic  effects  between  weapon  system  types 
are  not  considered  [Ref.  9],   Also  inherent  in  the  modelling 
of  heterogeneous  force  combat  is  the  need  to  introduce  an 
allocation  factor.   This  allocation  factor  may  be  defined  as 
the  fraction  of  the  firers  of  a  particular  type  that  engage 

77 


a  particular  target  type.   The  allocation  factor  concept 
will  be  discussed  in  greater  detail  later  in  this  chapter. 

One  theoretical  concept  for  modelling  the  attrition  pro- 
cess between  two  opposing  forces  is  to  consider  attrition  as 
a  nonstationary  Markov  process  where  the  states  of  the  pro- 
cess are  defined  as  the  numbers  and  types  of  surviving  com- 
batants and  the  ranges  between  opposing  combatants.   The 
transition  probabilities  between  the  Marlcov  states,  under 
this  concept,  depend  on  the  ratio  at  which  surviving  combat- 
ants or  weapon  systems  inflict  casualties  on  a  particular 
type  of  opposing  combatant.   Thus  the  nonstationary  Markov 
process  model  provides  a  means  of  relating  the  sequence  of 
states  and  the  times  of  the  state  transitions  observed  dur- 
ing a  battle  to  attrition  rate  coefficients.   In  turn,  this 
allows  for  the  estimation  of  attrition  rate  coefficients  on 
the  basis  of  the  observed  sequence  of  states  and  transition 
times.   Andrighetti  [ 3ef .  10  ]  developed  such  a  nonstationary 
Markov  model  that  related  weapon  system  effectiveness  to  the 
time  sequence  of  casualties  observed  in  a  two-sided,  hetero- 
geneous force  land  combat  simulation.   There  is  a  flaw,  how- 
ever, in  Andrighetti* s  work  that  readers  should  note.   He 
implies  that  his  methodology  is  easily  extended  to  time-de- 
pendent attrition  rate  coefficients.   This  is  not  true. 
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Andrighetti 's   thesis   does   contain   a    good    discussion  on 
the  theoretical   background   which  allows    for    the   maximum 
likelihood   estimation   of   parameters   used   in    attrition   rate 
equations.      An    important   observation   made  by    Andrighetti   in 
his  thesis   is  that   the  difference-differential  equations 
obtained   from   the  Chapman-Kolmogorov  equations  of  the 
nonstationary   Markov    model  are    susceptible  to  a   general 
solution  only   in   special   cases   as  Gordon   Clark  showed   in   his 
doctoral  dissertation.      Clark  used  the   Chapman-Kolmogorov 
equations  to   develop  expressions  for   time  state   probabili- 
ties and  the   expected   survivors  in   a   heterogeneous  force 
battle. 

The    extended   version   of  the   DAMSTAC    model  which   shall   be 
called   DAMSTAC    II    should   be   a   heterogeneous   force  analytical 
model    that   uses    STAR   output   to    compute   parameter    values. 
These    parameter   values  should   be  allowed   to    change  as    the 
distance  or  range   between   opposing   forces   vary  during   the 
battle.      In   DAMSTAC,    an   attempt    was   made   to    allow  the   param- 
eter  values  to   change    with   respect    to   time.      However,    it   was 
discovered   that   the   maximum   liklihood   estimators    developed 
for   the   critical   parameters    in   the   attrition   rate   equations 
could   not   be    used   for    time   dependent   variables.      Two    feasi- 
ble  alternatives   that    allowed   for   variable   parameter    values 
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were  considered.   Both  alternatives  are  recommended  for 
future  evaluation.   The  first  alternative  looks  at  dividing 
the  battle  into  phases  based  on  battle  intensity  as  measured 
by  the  number  of  shots  divided  by  the  number  of  live  firers. 
(Refer  to  Section  III-E)   The  second  alternative  for  divid- 
ing the  battle  into  phases  is  based  on  the  range  between 
combatants.   Here  also  the  parameters  can  change  from  phase 
to  phase.   This  technique  will  work  assuming  that  the  com- 
batants on  each  side  move  as  a  single  force  and  that  units 
are  located  by  their  centers  of  mass.   For  example,  a  battle 
can  be  divided  into  three  distinct  time  intervals  or  phases 
of  combat:  an  artillery  prep  phase  during  which  all  combat 
vehicles  and  weapon  systems,  other  than  artillery,  are  not 
participating;  a  moderate  to  long  range  battle  phase;  and,  a 
terminal  shorter  range  battle  phase.   Parameter  values  can 
then  be  computed,  using  procedures  to  be  discussed  later, 
for  each  battle  phase.   This  second  alternative  is  used  in 
the  COMANEW  model  and  is  the  one  recommended  for  DAMSTAC  II. 

A.   THE  ATTRITION-RATE  FUNCTIONAL  FORM 

The  attrition  rate  functional  form  that  will  ultimately 
be  chosen  for  the  heterogeneous  force  model  is,  of  course, 
the  prerogative  of  the  modeller.   The  authors  of  this  the- 
sis, however,  would  like  to  recommend  two  related  attrition 
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rate  functional  forms  for  consideration.   The  first,  which 
shall  be  labeled  as  Functional  Form  1  (FF1) ,  is  based  on  rhe 
assumption  that  during  the  course  of  a  battle  fire  is 
directed  at  the  highest  priority  target  surviving  and 
acquired.   Also  it  is  assumed  that  the  relative  priority  of 
target  weapon  types,  as  considered  by  each  firer,  remains 
constant  throughout  the  engagement  and  that  this  relative 
priority  is  identical  for  each  firing  weapon  type.   Func- 
tional Form  1,  developed  by  Gordon  Clark  for  use  in  the 
COMAN  model,  can  be  written  as  follows: 


f*(mn)  =  t     °!iUn  p..'"/(l-p'"i) 
J  lui 


(16) 


where 

f^  (m,n)  =  conditional  casualty  rate  of  BLUE  weapons  of 
type  j  during  phase  i  of  the  battle  given  strengths  of  m 
and  n. 


f^  (m,n)  =  analog  of  f^.  for  RED  weapons  of  type  k. 
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n.    =   the   number  of   surviving   RED   weapons  of  type  Ic 


m-     =   the   number  of   surviving    BLUE   weapons   of  type    j 

%j  =    kill   rate   of  a    BLUE  firer  of  type   j   for  an  acquired 
RED   target   of   type   Ic  in   the   i»th   time    interval. 


^Uk"    kill  rate   for   a  red-f irer/blue-target   combination 
corresponding   to  jSj^ 


q; =  the  probability  that  a  specific  RED  weapon  is  unac- 
quired by  an  individual  BLUE  firer  in  the  i«th  phase  of 
the   battle. 


p;  =    the   probability    for   a   specified   blue-target/red-firer 
combination   corresponding   to    q; . 


mj    =   the   number  of   surviving    BLUE   weapons   of  higher    prior- 
ity  than   a   type    j   weapon. 
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n',^  =  the  number  of  surviving  RED  weapons  of  higher  prior- 
ity than  a  type  k  weapon. 

The  reader  should  note  that  since  the  assumption  is  made 
that  a  firer  will  engage  the  highest  priority  target  he  can 
acquire,  the  probability  of  a  BLUE  firer  engaging  a  RED 
weapon  of  type  Ic  is  qj'*'(1-q]S  .   Similarly,  the  probability  of 
a  RED  firer  engaging  a  BLUE  weapon  of  type  j  is  p'^CI-j^"^  . 
These  two  probabilities  may  be  regarded  as  the  allocation 
factors  for  Functional  Form  1. 

All  of  the  critical  parameters  necessary  to  compute  the 
conditional  casualty  rates,  f  j^.  (m,n)  and  f^j^  (m,n)  ,  can  be 
obtained  either  directly  or  indirectly  (using  maximum  like- 
lihood procedures)  from  the  output  of  a  high  resolution 
simulation  such  as  STAR. 

If  one  is  not  willing  to  assume  that  target  priorities 
remain  constant  throughout  the  battle  then  FF1  will  not  be 
sufficient.   The  terms  in  Clark's  functional  form  which 
result  from  the  assumption  of  constant  target  priorities  are 
qj^and  p{^i  where  q-.^is  the  probability  that  all  surviving  RED 
weapons  of  higher  priority  than  type  k  have  not  been 
acquired  and  pp^  is  the  corresponding  probability  for  type  j 
blue  weapons. 
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One    way  to   extend   Clark's   formulation  in   order  that  the 
fixed   priority   may   be   dynamically  overridden   is   by  replacing 
the   terms   q."''    withGjj^    where  ej^^   is  defined  as   the    probability 
that   BLUE   weapon  type    j  engages   RED    weapon   type   k,    given 
that   BLUE   weapon   type    j  acquires  at    least  one   RED   weapon 
type   k.      Note    that 


^^^C"?***^ )    -    P(j   engages   k|  at  least   one  k   acquisition) 
P (acquire   at  least  one   k)    =   P(j   engages    k) 


Functional   Form    2    (FF2)    may    be   written   as; 


.w        .         .         Ju  ..       m: 


.f^.      (m,n)    =  2;  aj^ic    n^  5^  (1-P;  M  (17) 


.irt  ^ 


f^j^     (m,n)     =S/3,Vj    mj  e:i,  (l-qi**)  (18) 


where  6^^■    is  the  RED  firer  BLUE  target  probability  corre- 
sponding to  ej^^  and  all  other  parameters  are  as  defined  for 
FF1.   Just  as  PrqfQ!,/3  are  assumed  constant  over  the  entire 
phase  of  the  battle,  5     and  9  ,    for  each  j  and  k,  are  also 
assumed  constant  over  the  entire  interval. 
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Serious  consideration  should  be  given  to  incorporating 
FF2  or  some  similar  functional  form  into  DAttSTAC  II  because 
STAR  does  allow  target  priorities  to  change  during  a  battle. 
Target  priorities  in  STAB  are  based  on  the  types  of  targets 
on  the  firer's  detected  list  and  the  type  of  ammunition  the 
firer  has  available  for  use. 

It  may  also  be  useful  to  use  a  different  functional  form 
during  the  initial  phase  of  combat  where  attrition  rates  are 
more  a  function  of  the  exposure  or  vulnerability  of  targets 
than  they  ara  a  function  of  the  firer* s  capability  to 
acquire  targets.   Clark  recommended  the  use  of  Peterson's 
Logarithmic  Law  to  model  attrition  during  the  initial  phase 
of  the  battle.  [Ref.  11] 

B.   MAXIMOM  LIKELIHOOD  ESTIMATORS 

The  ability  of  the  DAMSTAC  II  model  to  provide  insighr 
into  the  interactions  represented  in  the  STAR  simulation  is 
based  on  the  estimation  of  conditional  casualty  rates  (  a 
and /3  )  and  other  parameters  from  STAR  data.   These  estimated 
parameters  should  be  constant  within  a  phase  of  the  battle, 
but  independent  and  unrelated  to  results  in  neighboring  time 
intervals.   The  independence  of  parameter  values  between 
time  intervals  implies  that  the  estimators  for  a  specific 
time  interval,  or  phase  of  the  battle,  can  be  defined  by 
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analyzing  a  sample  of   data   from  the  hi^h   resolution  simula- 
tion observed   during   that  single  -time   interval. 

The    maximum   likelihood   estimators    (MLS's)    for   the 
unknown   parameters  a, (3,   p  and   g   used   in   FF1    were   developed 
by  Gordon   Clark.      It   may   be  of    interest    to  the  reader   to 
note  that  simplified   versions   of  these   MLS»s    were  used   in 
the   DAMSTAC  modal   for    homogeneous  forces,    however,   since   we 
are  now   concerned   with  heterogeneous   forces,    the    "una- 
bridged"  versions   are    presented.      The   maximum   likelihood 


estimate   for    a    as   a    function   of  p   is: 


\.(D)      = 


n^o"-**    (1   - 


;[;n^D"'<*  (1  -  p'"^*  )(t^  -  ty., )         *^''    •*» 


(19) 

0  otherwise 

where     o^Api-  ^i^®   estimator   for   the  parameter  a    for  a   given   p. 

The   MLE   for   p   is   found  by   solving: 


aln   L(d)  g  (1   -   Cv)(m;^   -(m'^^   ^   m^^^)^^^^ 


'^*y 


(20) 


>(*l  jal   l^sl 

The   aquations   for  |^.    (q)    and   g   are  similar. 

The   es-imators   obtained   from  the   above   equations    result 
from   the   data   from   a   single  replication   of  the  combat   simu- 
lation.     Jlora   precise    estimates  are   obtainable  by  using  a 
larger    sample   consisting   of   a   number   of    observations    during 
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the  interval  in  question  from  a  number  of  independent 
replications  of  the  STAB  simulation.   This  modification 
entails  an  additional  summation  (^   ,    where  x  =  the  total 
number  of  replications)  in  each  of  the  estimating  equations. 
The  data  needed  for  the  computations  of  the  MLE*s  include; 

y  =  total  number  of  casualties 

t^  =  elapsed  time  since  the  beginning  of  the  interval 

until  the  y'th  casualty 

ty^,=  the  end  of  the  battle  phase  or  the  elapsed  time  since 

the  beginning  of  the  phase  until  the  battle  ended  if  it  is 

less. 

m^  =  the  number  of  BLUE  survivors  of  type  j  just  prior  to 

the  y'th  casualty 

nY^  =  the  number  of  RED  survivors  of  type  k  in  the  battle 

jus-  prior  to  the  y»th  casualty 

C^  =    1  if  the  y«th  casualty  is  BLUE  ;  0  otherwise 

f^  =  firer  weapon  type  for  the  y'th  casualty 

gy  =  target  weapon  type  for  the  y^th  casualty 

myj  =  the  number  of  BLUE  survivors  of  higher  priority  than 

■'^ype  J  just  prior  to  the  y'th  casualty 
b 

n^  =  RED  analog  of  ntyj 
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B..  =  the  number  of  BLUE  weapon-type  j  casualties  due  to 

RED  weapons  of  type  k 

One  technique  of  computing  the  MLE's  for  a   and  p  is  to 
numerically  determine  p  using  equation  (20)  and  then  substi- 
tute this  value  into  equation  (19)  for  j  =  1,2,...,b  and  for 
k  =  ^ ,2, .. , ,T    (b  =  the  total  number  of  BLUE  weapon  types  and 
r  =  the  total  number  of  RED  weapon  types).   A  recommended 
technique  for  numerically  solving  for  p,  the  estimate  for  p, 
can  be  found  in  Clark's  work  on  the  Dyntacs  combat  model 
[Ref.  12].   This  technique  is  similar  to  the  one  developed 
for  the  homogeneous  case  discussed  in  Chapter  III.   Basi- 
cally the  procedure  is  as  follows.   Since  p  represents  a 
probability,  it  must  satisfy  the  inequality  0<p<1.   There- 
fore, the  function -^  Ln  L(p)  can  be  evaluated  at  points  in 
the  interval  0<p<1  until  two  points  ,Pj  and  p^  ,  are  found 
such  that  the  product  of  -Jr  Ln  L  (Pj  )  and  -^  Ln  L  (p^)  is  less 
than  or  equal  to  zero.   Since  -^  Ln  L (p)  is  a  continuous 
function,-^  Ln  L  (p)  =0  for  some  p  satisfying  p,  ^  p  ^  P2  • 
This  interval  is  reduced  by  evaluating  -J —  Ln  L(p)  at  the 
midpoint  of  the  interval  (p,  ^p^^).   Repeating  this  process 
converges  to  a  solution  for  -^  Ln  L(p)  and  this  solution  is 
a  local  maximum.   This  technique  is  similar  to  the  bisection 
method  for  solving  nonlinear  problems. 
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with   respect   to  Functional   Form   2,    the  introduction  of 
the     0j.     allocation   terms   adds  considerable   complexity   to   the 
maximization   of  the  likelihood   function.      One  way  to   circum- 
vent this   problem   is   to  compute   a  "best    guess"   for  the  ej^^  's 
from  STAR   samples    (for  the   given  interval) •    This   best    guess 
can   then  be    substituted  into  the  ej|^  's   prior   to   maximization 
of  the    likelihood   function,    thereby   reducing   thee^j^'s   to 
constants.      The    maximization   can  then   be   performed,    as   in 
Clark's    original  equations,    with  respect   to    p,   g,a  ,  /3     . 

C.       THE    MODEL 

The  computational  procedures  that  will  constitute 
DAMSTAC  II  must  be  developed  by  the  modeller  who  chooses  to 
carry  on  the  work  initiated  in  this  thesis.   As  a  starting 
point,  however,  the  techniques  used  by  Clark  in  his  COMAN 
model  would  be  excellent  choices.   The  adequacy  of  Clark's 
procedures  in  DAMSTAC  II  given  a  particular  scenario  and 
force  mix,  to  predict  attrition  in  a  STAB  simulation  with 
the  same  scenario  and  force  mix  is  questionable  and  would 
have  to  be  tested.   Detailed  comparisons  and  analysis  of  the 
output  from  both  models  will  be  required  to  answer  the  ques- 
tion of  model  comparability  and  modifications  to  the  com- 
puter code  and/or  computational  procedures  in  DAMSTAC  II 
must  be  accomplished  as  necessary.   The  computational  cycle 
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employed  in  COMAN  is  based  on  the  assumption  that  the  dis- 
tribution times  between  casualties  in  the  high  resolution 
simulation  is  exponentially  distributed.   It  uses  Monte 
Carlo  procedures  to  generate  the  time  until  the  next  casu- 
alty and  also  to  find  the  casualty  weapon  type. 

The  following  steps  constitute  the  computational  cycle 
recommended  for  DAMSTAC  II. 

Step  1.   The  first  step  in  the  cycle  sets  rhe  time  equal 
to  zero  and  sets  the  battle  phase,  designated  by  the  letter 
i,  equal  to  1. 

Step  2.   Next  the  total  conditional  casualty  rate, 
\^(m,n)  is  computed  from  the  equation 
^W)  =|:ft-Kn)+  |;^f'i;^(m,n^  for  i=l,2,...,  I   (21) 
where  fj^.  (m,n)  and  fjtj^  (m,n)  are  the  conditional  casualty 
rates  in  phase  i  pf  the  battle,  given  force  strength  m  and 
n,  for  BLDE  weapon  type  j  and  RED  weapon  type  k,  respec- 
tively. X'  (m,n)  is  the  parameter  of  the  exponential  distri- 
bution describing  the  times  to  the  next  casualty. 

Step  3.   Now,  using  the  exponential  distribution  with 
parameter  \'(m,n),  one  Monte  Carlo's  for  the  time  interval 
to  the  next  casualty.   Designate  this  time  interval  as  t* . 
If  the  next  casualty  occurs  in  another  time  interval,  ie.  if 
t  +  t»  >  t;  where  t-,  represents  the  end  of  the  i'th  battle 
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phase,  then  i  is  increased  by  1,  the  time  is  set  equal  to  t 
and  the  program  is  directed  back  to  step  2  where  \'**(m,n)  is 
selected  as  the  distribution  parameter.   If,  however,  t  ♦  t* 
<  t;  then  the  next  step  is  to  Monte  Carlo  for  the  casualty 
weapon  type. 

Step  4.   The  weapon  type  for  the  next  casualty  is  repre- 
sented by  the  random  variables  C^  and  C^  which  are  defined 
as 

C^  =  j  if  the  casualty  weapon  is  BLUE  weapon  type  j 

0  if  the  casualty  is  RED 
C,.  =  k  if  the  casualty  weapon  is  RED  weapon  type  k 
0  if  the  casualty  is  BLUE 
The  Monte  Carlo  procedure  uses  the  conditional  casualty 
weapon  probability  distributions  for  C^    and  C^  .   These 
probabilities  are  computed  by 

P'  (Cb  =  J|C  =  1,m,n)  =  f^*  (inrn)/x(rn,n^ 
p'  (Cr  =  k|C  =  1,m,n)  =  ftjj  (m,n)/ A(m,n) 

for     i=    1,2,...,    I   and  when   variable  C   is  set   equal   to   1    if 
a  casualty    is  determined   to   occur  in   interval   i   from    the 
previous   steps.      The    probability  that   the  casualty   is   of   a 
particular    type    is   computed   for   each   RED   and    BLUE  type 
weapon.      The   sum  of  these   probabilities    should  equal    1. 
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step  5.  The  next  step  is  to  Monte  Carlo  using  the  Oni- 
form  (0,1)  distribution  to  determine  the  casualty  type.   An 
example  may  help  in  understanding  the  above  procedure.   Sup- 
pose there  are  two  types  of  BLUE  weapons  and  two  types  of 
RED  weapons.   Assume  the  conditional  casualty  probabilities 
for  each  weapon  are  computed  using  the  equations  shown  pre- 
viously and  they  are  0.25  for  each  weapon  type,  b1,  b2,  r1, 
r2.   If  one  orders  the  weapon  types  in  some  arbitrary  fash- 
ion such  as  b1,  r1,  b2,  r2,  and  then  accumulates  the  prob- 
abilities on  the  interval  (0,1),  one  can  assign  an  interval 
to  each  weapon  type.   The  example  would  yield  the  following 
int  arvals: 

b1  -  (0,0.25)  ; 

r1  -  (0.26,0.50)  ; 

b2  -  (0.51,0.75)  ; 

r2  -  (0.76,1 .00)  . 
The  uniform  random  number  can  then  specify  the  casualty 
weapon  type.   If  the  random  number  is  0.60  then  the  casualty 
weapon  type  in  our  example  is  b2. 

Step  6.   After  the  casualty  weapon  type  has  been  found 
the  next  step  is  to  decrement  the  number  of  that  type  weapon 
by  1.   Now  increase  the  value  of  t  by  t' ,  the  ^ime  to  the 
next  casualty.   If  t  is  greater  than  or  equal  to  t^  which 
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represents   the    ending   time  of  the  last    phase   cf   the   battle, 
or   if    either   side   is   annihilated   ,    the   battle   has   ended  and 
the   program  terminates.      If  t   is  less   than   t^    then  go   to 
step   2. 

Observe  that   the   above   computational   procedures   are  sto- 
chastic  in  nature   and   different   replications   of   an  engage- 
ment  will   likely   yield  similar    but    different   results.      If 
this  method  is   chosen   for   DAMSTAC  II,    variance  reducing 
replications   for   each   engagement  should   be   performed. 

D.       THE    STAR    POSTPROCESSOB 

It    is   very   important   that   anyone   who   attempts  to    develop 
the   DAMSTAC   II    model    understand   and   be   able    to  manipulate 
the   STAR   simulation   model.      In    the   initial  attempts   to 
extract    data    from   STAR,    certain   subroutines    must    be   turned 
off.      These  routines    include  the  SmoJce,    Engineer,    Close-air 
support,    etc.    whose   effects   on    parameter    estimates  used  in 
DAMSTAC    II   have   not  been   thoroughly   analyzed.      The   number   of 
types   of  combatants   played   in   the  simulation    will   have  to   be 
regulated   so    that   they  do   not   overpower    the    capabilities   of 
the   analytical    model.      It   is   recommended   that   no    more   than 
three    different    types   of    weapon   sysrems    for    each   opposing 
force   be   attempted   when   initially   formulating  the   heteroge- 
neous   force  model. 
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It    is   important   that   the   finished   version  of   the 
analytical   model   be  as  simple   and  transparent  as    possible 
and  not    encumbered   with    parameter  estimating   algorithms  and 
computations.      This   consideration   plus   the   amount   of    data 
that    must   be    extracted   from   the   STAB   output    require  the 
development  of    a   STAR    postprocessor    program.      As    an   illus- 
tration,   if  one   considers   a   simple   battle  divided   into  three 
battle    phases   with   only   three   types   of    weapon   systems   on 
each   side,    the    DAMSTAC  II   model,    as    it    has   been   described, 
requires   the   estimation   of    114    parameter   values   and  the  com- 
putation  of  six   attrition   rates. 

The  postprocessor  at  Appendix  E,  developed  for  the  homo- 
geneous model,  should  provide  a  good  foundation  for  the  more 
complicated  postprocessor   needed  for   DAHSTAC    II. 
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VII.   CONCLUDING  REMABKS 

The  concept  of  using  differential  equations  to  represent 
combat  began  with  Lanchester's  classical  eguations  and  has 
developed  into  a  very  dynamic  field  of  study.   The  theoreti- 
cal background  of  the  fitted  parameter  method  has  been  well 
documented  by  Clark  and  Taylor.   However,  the  practical 
application  of  the  theory  lacks  this  same  degree  of  discus- 
sion and  documentation.   This  paper  has  focussed  on  some  of 
the  techniques  of  applying  Clark's  theory  to  an  actual  high 
resolution,  Monte  Carlo  simulation.   The  problems  are  real 
and,  as  has  been  seen,  the  solutions  presented  are  not 
doctrine.   It  is  hoped  that  further  investigation,  elabora- 
tion and  documentation  of  this  modelling  technique  will 
occur  in  the  near  future. 

The  choice  of  which  functional  form  to  use  in  a  particu- 
lar analytical  combat  model  is  entirely  up  to  the  discretion 
of  the  model  builder.  In  choosing  an  attrition  rate  func- 
tional form  the  model  builder  should  consider  the  scenario 
portrayed  in  the  model  and  the  assumptions  that  are  compati- 
ble with  it. 

Finally,  it  is  recommended  that  the  entire  program  for 
the  model  be  converted  to  the  UCSD  Pascal  language  on  the 
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Apple  II.   The  basic  structure  of  Pascal  is  more  conducive 
to  combat  modelling  simulation  than  the  structure  of  the 
BASIC  language.   The  concept  of  entities  and  attributes 
found  in  SIMSCRIPT  II. 5  can  be  closely  emulated  in  Pascal 
while  this  is  not  as  easily  accomplished  or  clearly  defined 
in  BASIC. 

The  modular  structure  of  procedures  and  functions  in 
Pascal  is  also  an  advantage  in  the  enhancement  process  of 
model  building.   Modellers  can  use  one  procedure  to  repre- 
sent a  particular  action  or  phenomenon  occurring  on  the  bat- 
tlefield.  Any  action  which  is  unable  to  be  modelled  can  be 
left  in  the  program  as  a  dummy  procedure  for  later  work. 

Additionally,  the  Pascal  language  allows  more  descrip- 
tive variables  to  be  used  than  the  two-character  variables 
used  in  Applesoft.   This  would  make  the  program  itself  more 
transparent  and  easier  to  modify  and  update. 

The  DAMSTAC  model  is  clearly  in  its  infancy  stage. 
Through  continued  enhancement,  this  model  can  contribute  to 
the  ongoing  work  in  the  field  of  combat  analysis. 
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APPENDIX    A 

LIST    OF    VARIABLES 
A   =  casualty  rate    for   RED   attritting   acquired  BLUE   targets 
AN  =   angle   in    radians   measured   counter-clockwise   from  due 
east    made   by   the  line   between    the    current   unit 
location   and  the  next   coordinating   point. 
B  =   casualty   rate   for  BLUE  attritting   acquired  RED   targets 
BE  =   the   breakpoint    for   the   BLUE  force 
BM  =    movement    criteria   for    BLUE  units    as  a    fraction   of 

unit   strength 
BF (I)    =   number   of   RED  targets    being   fired  at   by   BLUE    unit   I 
BP  (I)    =   most   recently  passed  coordinating   point    for    BLUE 

unit  I 
BR  (I)    =   rate   of   movement   for   BLUE   unit    I  in    M/SEC 
BT(I,J)    =    attrition(casualties)    caused   by   RED  unit   J   on 

BLUE    unit    I   in   one   on  one   battle 
BX(I,J)    =    x-coord   of    coordinating   point   J   for  BLUE    unit   I 
BY (I, J)    =    y-coord   of   coordinating    point   J   for   BLUE   unit   I 
BO (I)     =  current   total   BLUE   force   strength   for   time   period   I 
B0(0)     =   initial   BLUE    force   strength 
B1    =    number   of   coordinating    points   for    BLUE    (including 
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initial  location) 
B1  (I)  =  initial  strength  of  BLUE  unit  I 
DB(I)  =  number  of  casualties  or  amount  of  attrition  for 

BLOE  unit  I 
DH  =  distance  travelled  by  a  unit  in  time  increment  DT,  at 

speed  ai 
DI  =  distance  between  current  unit  location  and  next  turn 

point 
DR ( J)  =  number  of  casualties  or  amount  of  attrition  for  RED 

unit  J 
DT  =  time  increment 
DX  =  X  component  of  distance  between  current  unit  location 

and  next  coordinating  point 
DY  =  y  component  of  distance  between  current  unit  location 

and  next  coordinating  point 
DI  =  X  component  of  distance  moved  in  time  increment  DT,  at 

speed  RR  (I)  or  BR (I) 
D2  =  y  component  of  distance  moved  in  time  increment  DT,  at 

speed  RRtI)  or  BR  (I) 
M  =s  number  of  BLUE  units 
MT  =  max  time  of  simulation 
N  =  number  of  RED  units 
N$  =  file  name  (used  recursively) 
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P  =   probability  that    a   BLUE  target   is   nor   acquired   by 

a    RED    anit 
PI    =    3. 141592654 
PR   =   Status   of   printer;    0  =   printer   is   off 

1   -   printer   is   on 
Q  =   probability  that   a   RED  target   is  not  acquired  by 

a    BLUE   unit 
RE  =    the   breakpoint    for   the   RED  force 

RF  (J)    =   number   of   BLUE   targets   being  fired    at  by   RED    unit   J 
RG   =   range   between   a    BLUE   unit   and   a   RED  unit 
RM   =    movement    criteria   for   RED   units   as   a   fraction   of 

unit    strength 
RP(I)    =   most   recently  passed   coordinating   point    for    RED 

unit   I 
RR(iy    =   rate  of   movement    for   RED  unix    I   in    M/SEC 
RT(J,I)    =   attrition    (casualties)    caused  by   BLUE    unit   I  on 

RED    unit  J    in    one   on   one   battle 
RX(I,J)    =    x-coord   of    coordinating   point   J   for  RED   unit  I 
RY(I,J)    =   y-coord   of    coordinating   point  J   for  RED   unit   I 
RO(I)     =   current  total  RED   force  strength  for   time  period   I 
RO (0)    =   initial  total   RED   force  strength 
R1    =    number   of   coordinating   points   for   RED    (including 
initial   location) 
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R1(I)  =   initial   strength   of   RED   unit   I 

SB  (I)  =   size   of   BLUE   unit   I 

SR(I)  =   size  of  RED    unit   I 
T   =   simulation    time 

XB(I)  =   current   x-coord   for   BLUE  unit   I 

XR (I)  =  current   x-coord   for   RED  unit  I 

YB (I)  =  current   y-coord   for   BLUE  unit   I 

YR  (I)  =   current   y-coord   for   RED  unit   I 
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APPENDIJC  B 

NOTES  ON  THE  COMPUTER  MODEL 

1.  Data  which  have  been  hardwired  into  the  model  are; 
A  =  2.4 

B  =  0.03 

BE  =  0.3 

BT  =  0 

DT  =  30 

MT  =  1500 

P  =  0.99 

PI  =  3. 14159265U 

Q  =  0.93 

RE  =  0.3 

2.  Once  contact  is  made  (range  <  2700) ,  the  rate  of  advance 
decreases  to  2  meters/second. 

3.  The  option  exists  within  the  program  to  allow  the  move- 
ment of  units  based  on  a  criteria  of  fractional  force 
strength.   When  unit  strength  goes  to  50*  (movement  crite- 
ria) ,  the  unit  will  move  to  an  alternate  location  (for  BLUE 
units)  or  retreat  along  its  path  of  advance  (for  RED  units) , 
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APPENDIX   C 

UllLlll    PROGRAMS 

This   appendix   provides   a  brief   description  of   the    util- 
ity  programs   associated   with  the   Aggregated   Model  of    STAR   on 
the   Apple   II   Computer.      These    programs   provide  the  capabil- 
ity  to    the   user    to   transfer   the   necessary   data   from   the  STAR 
model   to  this    model   with   the   minimum   effort. 

The    utility   programs   available   are: 
Force   Maker 
Results   Reader 
Both   programs    are   interactive   with   the   user,    who    answers 
questions   provided  by   the   computer.      These   programs   use  the 
commands   for   the    Apple  II    Disk   Operating   System    (DOS)     which 
allow    them    to   read/write   Sequential   Text    Files   and  Random 
Access    Files   onto   the   disk.      Familiarity    with  these   types   of 
files   and   the    DOS    would   improve   the   understanding   of    these 
utility    programs.    [Ref.    13] 

In    the   following    program  descriptions,    the   variable  name 
used   in   the  program   is  listed   in  parentheses    after  its 
description. 


102 


A.   FORCE  MAKES  PROGRAM 

This  program  will  write  data  files  containing  informa- 
tion about  the  two  opposing  forces. 

Force  Maker  produces  a  Sequential  Text  File  containing: 

Force  size  (number  of  units)  (H) 

The  strength  of  each  unit  (SB  (I)) 

The  number  of  coordinating  points  for  movement  of  each 

unit  (the  first  coordinating  point  is  always  the  initial 

location  of  the  unit)  (B1(I)) 

The  grid  coordinates  of  each  coordinating  point  for  each 

unit  (BX(I,J),  BY  (I,  J)) 

The  user  has  the  option  to  make  a  file  for  the  BLUE 
force  or  the  RED  force.   After  all  the  data  has  been 
entered,  a  data  check  routine  will  print  the  information  for 
all  units,  so  that  a  check  for  errors  can  be  made.   After 
the  user  has  checked  the  data  and  identified  the  name  for 
the  file,  the  program  will  record  the  file  on  the  diskette. 

A  copy  of  the  program  is  included  as  Listing  B  in  Appen- 
dix F  and  a  sample  session  using  this  program  is  attached  as 
Figure  12. 
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]RLJN    FGRCE    ^^AK£R 

1    i^AKE    ~:L£    ?^GR    5L^£    FCRCE5 
i    •'-AKE    FILE    FCF: 


;y!-!AT    :5    THE    FGRCE    SIZE    CF    3LJE    3 

FCR    UMIT    i 

hCW  MANY  CGGRDIMATIMG  PGINTS  (INCLUDING  START  LGCATIGN)'^  1 

ENTER  SRID  GGCRDS  FGR  CGCRDINATING  POINT  *  1 

v™!-'-!.-!?;. PS  T  (z     er/Li  -"r 

'.    K^^lmiF'   s^  ^  .^    ^w  ^  1  <3  ■_ 

T   _  w  v^  ~.  w   ^  ^   ^  .' '.'  O  «^  %^ 


HGW  i^ANV  GGGRDINATING  PGINTS  'INCLUDING  START  LGCATIQN)  ^  I 

ENTER  GRID  GGGRDS  FGR  GGCRDINATING  ='GINT  #  1 
X-CGCRD  IS  5^753 
Y-CGGRD  IS  98120 

-GR  UNIT  3 

HCVJ  r^ANY  GGCRDINATING  POINTS  (INCLUDING  START  lGCATIGN/-^  1 

ENTER  GRID  GGGRDS  FGR  CGGRDINATINQ  POINT  #  i 
X-GGGRD  IS  37^43 
Y— GGGRD  IS  ^"'S93 

It.^S  DATA  CHECK  -K''^^ 


UNIT  FT  A-GGCRD  ^-GGCRD 
i    L    34ij3    10Ci:S5 

2  1    347SS    931 GO 

3  L    5"'443    ^7393 
ANY  CHANGES^  (Y/N)  N 
NAi-E  GF  ^Z'^E^    THESIS 

;  r';AKE  ~I:_E  FGR  3LUE  FCRC 
Z  r^AKE  "ILE  ~GR  RED  FGRCE 


iiMTr'-C^   ,". 


Figure    12:      Sample   Session   of   Force   aaJter  Program 


104 


B.       RESULTS    READEB    PROGRAM 

This  program   reads  the   data   file  of   time    intervals  and 
corresponding   BLUE   and  RED   force   levels    produced    by   the 
model.       Results    Reader  produces   two   sequential  text   files: 
the   first   file   with   suffix   ".COEFS"   contains: 

Conditional  Casualty  rate,   a.    (A) 

Conditional  Casualty  rate,  /?     (B) 

Probability  that   a    ELOE  target  is   not    acquired    by   a    RED 

firer    (P) 

Probability  that  a  RED  target  is  not  acquired  by  a  BLUE 

firer  (Q) 
The  second  file  with  suffix  ". RESULTS"  contains: 

The  number  of  data  points  for  BLUE  force  level  versus  time 

(N) 

The   number   of    data    points   for   RED   force  level    versus   time 

(N) 

Tims  of  the  battle  at  time  interval  I  (T) 

BLUE  force  level  corresponding  to  time  interval  I  (BO(I)) 

Time  of  the  battle  at  time  interval  I  (T) 

RED  force  level  corresponding  to  time  interval  I  (RO(I)) 
The  ".RESULTS"  file  contains  some  duplication  of  data  for 
the  fDllowing  reasons.   With  the  ".RESULTS"  file  written  in 
this  format,  and  free  of  additional  data  such  as  that  in  rhe 
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".COEFS"  file,  it  can  be  used  directly  by  the  Apple  Plot 
Program  to  produce  a  plot  of  the  force  levels  versus  time. 
A  listing  of  the  program  is  included  as  Listing  C  in  Appen- 
dix P.   A  sample  plot  of  the  DAaSTAC  output  data  produced  by 
the  Apple  Plot  Program  is  included  as  Figure  13.   A  sample 
of  th2  output  of  DAaSTAC  from  Hesults  Header  is  included  as 
Figurs  iu. 
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Figure    13:      Sample    Plot   of   aesults   Header    Program  Outpur 
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Figure  14:   Sample  Output  from  Hesulrs  Reader  Program 
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APPENDIX  D 

THE  TERRAIN  MODEL 

In  addition  to  the  basic  computer  model,  work  was  done 
to  input  the  STAR  terrain  data  base  into  the  Apple  II  compu- 
ter.  This  effort  was  done  during  the  initial  phase  of  the 
thesis  work  based  on  an  erroneous  assumption  that  the 
authors  made  concerning  the  functional  form  of  the  attrition 
rates.   This  assumption  was  that  the  terrain  of  the  battle- 
field and  line  of  sight  computations  must  be  explicitly 
played  in  the  model.   In  the  given  functional  forms  for  the 
attrition  rates,  line  of  sight  and  terrain  are  accounted  for 
by  the  probability  values  of  p  and  g. 

Before  this  assumption  was  found  to  be  false,  four  util- 
ity programs,  and  a  model  that  computed  the  elevations  of 
the  units  were  developed  and  are  included  in  this  document 
for  possible  future  use.   If  further  work  is  done  with  a 
model  on  the  Apple  II  computer,  the  ability  to  input  the 
STAR  terrain  into  the  Apple  II  may  become  an  important  part 
of  the  model.   Examples  of  models  that  might  use  the  STAR 
terrain  in  the  Apple  II  are  small  aggregated  models  with 
functional  forms  that  explicitly  play  line  of  sight,  or 
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small,  high-resolution  models  that  are  suitable  on  the  Apple 
II  computer. 

A.   THE  TERRAIN  UTILITY  PROGRAMS 

The  utility  programs  that  allow  the  STAR  terrain  to  be 
input  in  the  Apple  II  computer  are  listed  below. 

Map  Maker 

Hill  Maker 

Woods  Maker 

Forest  Maker 

These  four  programs  have  two  portions  each:  a  file  maker 
and  a  file  reader.   In  general,  the  file  maker  portion 
allows  the  user  to  input  a  new  file  or  update/change  an  old 
file.   The  file  reader  allows  the  user  to  print  out  on  the 
screen  or  to  a  printer  the  selected  file. 

?or  the  Apple  II  combat  model,  a  STAR  terrain  file  (eg. 
FULTERR  file)  is  divided  into  its  four  main  parts: 

the  list  of  hill  reference  numbers  in  each  grid  square; 

the  hill  parameters; 

the  list  of  forest  reference  numbers  in  each  grid  square; 

the  forest  ellipse  fitting  parameters. 
Each  part  is  created,  respectively,  by  one  of  the  four  util- 
ity programs:  Map  Maker,  Hill  Maker,  Woods  Maker,  and  Forest 
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Maker.      When   these   parts   are  being  created  on  the   Apple 
Disk,    the   same   general  file   name  must   be   used  to   identify 
the   data's   terrain    (eg.   FOLDA) .      Each   of   the   utility    pro- 
grams   will   add   its   own  suffix  to  the  general   file   name.      For 
example,    for   the   Fulda  terrain,    the   four   utility    programs 
will   produce    files   named: 

FULDA. MAP 

FULDA. HILLS 

FULDA. WOODS 

FULDA. FORESTS 
This   suffix  identifies  the   file   as   a   particular   portion  of 
the   terrain   data,    and    will  allow   the   Apple   computer   to 
recognize   it    from   the   other    files   with   the   same   terrain  name 
when   reading   the   file    for   the   combat  model. 
1.      Ma£  Maker   Program 

This   program    will   read   and   write   data   files   contain- 
ing  a   catalog   by  grid   sguare  of   the   hills   which   influence 
the   terrain   in   that   grid   sguare.      This   program   produces   an 
eguivalent   file    to   that   in   STAfi   called  LIST. H.      Map   Maker 
produces   a   Seguential   Text   File   containing: 

The   size   of   the   battlefield    (normally    10   by    10)     (L,M) 

The    lower   left   grid   square   in    1000    meter   units    (LX,Lir) 

The   base   elevation   of  the   battlefield    (HL  (0,0,0)) 
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For  each  grid  square  I, J: 

(1)  The  number  of  hills  influencing  the  terrain 
(HL(I,J,0)) 

(2)  The  reference  number  of  these  hills  (HL(I,J,K)) 

The  user  has  the  option  to  maXe  a  new  file  or  access 
an  old  file  for  updating  or  changes.   For  new  files,  the 
program  will  start  with  the  lower  left  grid  square  and  pro- 
ceed from  west  to  east  then  south  to  north.  The  data  can  be 
saved  after  the  information  for  a  grid  square  has  been 
entered.   For  old  files,  the  program  will  access  it,  print 
out  the  file  on  the  screen  and  ask  for  the  starting  grid 
square.   The  grid  square  index  is  used  to  identify  in  which 
grid  square  to  start.   The  grid  square  index  is  an  integer 
from  1  to  L  in  the  east-west  direction,  and  1  to  W  in  the 
south-north  direction  of  an  L  by  M  sized  battlefield.   The 
lower  left  grid  square  is  1,1;  the  lower  right  is  L,1;  the 
upper  left  is  1,H;  and  the  upper  right  is  L, W. 

A  copy  of  the  program  is  attached  as  Listing  D  in 
Appendix  F  and  a  sample  file  is  attached  as  Figure  15. 
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Figure  15:   Sampla  File  for  Map  Maker  Program 
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2.   Hill  Maker  Program 

This  program  reads/writes  data  files  containing  the 
fitting  parameters  for  each  hill  on  the  battlefield.   Hill 
Maker  produces  a  Random  Access  Text  File  containing: 
The  number  of  hills  on  the  battlefield  (N) 
For  each  hill  I: 

(1)  X  coord  of  the  center  location  of  the  hill  (Ei(I,1)) 

(2)  I  coord  of  the  center  location  of  the  hill  (H(I,1)) 

(3)  The  elevation  of  the  hill  top  measured  from 
zero  =  sea  level  (H(I,3)) 

(4)  The  orientation  angle  of  the  ellipse  measured  in 
degrees  counter-clockwise  from  east  to  the  major  axis. 
(H(I,4)) 

(5)  The  eccentricity,  defined  as  the  ratio  of  major 
axis  length  to  minor  axis  length  (H(I,5)) 

(6)  The  spread,  defined  as  the  distance  in  meters  meas- 
ured along  the  major  axis  from  hill  center  to  the  contour 
line  which  is  50  meters  down  from  the  peak  (H(I,6)) 

(7)  The  maximum  height  of  the  normal  curve  describing 
this  hill  mass  (H(I,7)) 

(8)  The  vertical  distance  measured  from  the  peak  beyond 
which  this  hill  is  not  considered  in  computations  (H(I,8)) 
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A  copy  of  the  program  is  included  as  Listing  E  in 
Appendix  F,  and  a  sample  file  is  attached  as  Figure  16. 
3.   Woods  Maker  Program 

This  program  will  read  and  write  data  files  contain- 
ing a  catalog  by  grid  square  of  the  woods  which  influence 
the  terrain  in  that  grid  square.   This  program  produces  an 
equivalent  file  to  that  of  Map  Maker.   Woods  Maker  produces 
a  Sequential  Text  File  containing: 

The  size  of  the  battlefield  (normally  10  by  10)  (L,W) 
The  lower  left  grid  square  in  1000  meter  units  (LX,LY) 
For  each  grid  square  I, J: 

(1)  The  number  of  forests  influencing  the  terrain 
(FL(I,J,0)) 

(2)  The  reference  number  of  these  forests  (FL(I,J,K)) 
Woods  Maker  contains  the  same  options  as  those  found 

in  Map  Maker.   A  copy  of  the  program  is  attached  as  Listing 
F  in  Appendix  F,  and  a  sample  file  is  attached  as  Figure  17. 
^«   Forest  Maker  Program 

This  program  reads/writes  data  files  containing  the 
fitting  parameters  for  each  forest  on  rhe  battlefield.   For- 
est Maker  produces  a  Random  Access  Text  File  containing: 
The  number  of  forests  on  the  battlefield  (N) 
For  each  forest  I: 
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(1)  X  coord  of  the  center  location  of  the  forest 
(F(I,1)) 

(2)  Y  coord  of  the  center  location  of  the  forest 

(F(I,1)) 

(3)  The    height  of  the  trees  in   the    forest    (F(I,3)) 

(U)    The   orientation   angle   of  the   ellipse   describing   the 
forest   measured   in   degrees  counter-clockwise   from   east   to 
the   major  axis    (F(I,U)) 

(5)  The    length  of   the   semi-major  axis   of  the  ellipse 
(F(I,5)) 

(6)  The    length   of  the   semi-minor   axis   of  the   ellipse 

(F(Ir6)) 

Forest  Maker  contains  the  same  options  as  those 
found  in  Hill  Maker.   A  copy  of  the  program  is  included  as 
Listing  G  in  Appendix  F,  and  a  sample  file  is  attached  as 
Figure  18. 

B.   THE  TERRAIN  MODEL  PROGRAM 

A  copy  of  the  model  which  explicitly  plays  terrain  is 
included  in  Appendix  F  as  Listing  H.   A  line  drawn  along  the 
edge  of  the  program  delineates  the  terrain  model  portion 
which  was  added  to  the  basic  model.   The  additional  varia- 
bles required  to  run  the  terrain  model  are  also  included  at 
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Figure    17:      Sample    File   from    Moods   daksr  Program 
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the    end   of  this   appendix.      Professor  James   K.    Hartman   [Ref. 

14]   provides   the   details   of   the   methodology    for   line    of 

sight   and  elevation   using   the  STAR   terrain, 
equivalent   STAR   variables  are   in     (  ) 
Al    =   computing   parameter   for   a   hill    (A) 
B1    =   computing   parameter   for  a   hill    (B) 
CR   =   critical    value   parameter   for   a   hill    (CRIT.H) 
CI    =    conversion  of   orientation   angle   to   radians 
C2   =    (spread  of  a    hill)  squared 

FI    =    value   of   the   elevation   of   a  point    due   to  a    hill    (FI) 
H(I,J)    =  the  J-th    parameter  for  hill   I;    the   parameters  are; 

1  =   x-coord   of  the  center    of   the  hill    (XC.H) 

2  =   y-coord  of  the  center    of   the   hill    (YC.H) 

3  =   peak  elevation  of   the    hill    (PEAK.H) 

4  =   orientation   angle   of   the   hill    (ANG. H) 

5  =   eccentricity   of  the   hill    (ECC.H) 

6  =   spread    of    the    hill    (SPED.H) 

7  =   max    height   of   the    hill    (HT.H) 

8  =   cut    height   of   the    hill    (CDT.H) 
H(I,0)    =   status  of   hill   I; 

0  =   hill  I    has   not   been   checked   for   LOS 

1  =   hill  I   has   been  checked  for   LOS 
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HL(I,J,K)    =   the   number   of  the   K-th   hill   in    grid   square  I, J 
which   influences   the  terrain  in    that   grid 
square   I   =    1    to   L;    J  =    1    to   W;    K  =    1 
to    HL(IrJ,0) 
HL(IrJ#0)    =   the  total  number  of  hills   in  grid  square   I, J 

which   influence  the  terrain   in   that    grid   square 
HL  (0,0,0)    =   the  base   elevation  of   the   terrain  map 
H2  =   temporary   variable   identifying  the   hill  being   checked 

for   line   of   sight 
L  =   number  of   grid  squares   on    the   terrain   map   in  the 

X   direction    (NGRIDX) 
LX  =   x-coord   of  lower  left   grid  square   of   terrain  map  in 

1000    meter    units    (X.LO.BDRY) 
LY  =    y-coord   of   lower  left    grid  square   of   terrain  map   in 

1000    meter    units    (Y.LO.BDRY) 
MS   =   name   of   files  containing   terrain    data 
NH  =   total   number   of    hills   on    the   terrain   map    (NHILLS) 
PI    =   computing   parameter   for  a    hill    (PXX.H) 
P2   =   computing    parameter   for   a   hill    (PXY.H) 
P3   =   computing   parameter  for  a   hill    (PYY.H) 
QI   =    quadratic   function   of   the   ellipse    represnxing   the 
hill    (QI) 


120 


W  =   number  of   grid  squares   on   the  terrain   map  in  the   Y 
direction    (NGRIDY) 

X  =  current    x-coord   of  unit   whose  elevation    is   being 

computed    (in    lOQ   meter   units)     (X) 
XB(I)    =  current   x-coord   for   BLOE  unit  I 
xa(I)    =  current   x-coord   for   RED  unit   I 
XS  =   distance   in  x  direction   between   unit  location   and 
center   of   hill   mass    (XS) 

XI  =   relative    number   of   grid  squares   in   x   direction    of  a 

BLOE   unit   location   with   respect   to   the   lower  left 

grid  square 
X2  =   relative   number   of   grid   squares   in   x  direction   of  a 

RED   unit   location   with  respect   to   the   lower   left 

grid  square 
X3   =   number  of   grid   squares   in  x  direction    between   a   BLUE 

unit   and   a   RED    unit 
X«»,YU   =  current   grid   square  being   checked   for  line   of 

sight    in   LOS    subroutine 
Y   =  current   y    coord   of   unit   whose   elevation   is   being 

computed    (in    100    meter   units)     (Y) 
YB  (I)    =   current   y-coord   for   BLUE  unit   I 
YR  (I)     =  current   y-coord   for   RED   unit   I 
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YS  =  distance  in  y  direction  between  unit  location  and 
center  of  hill  mass  (IS) 

11  =  relative  number  of  grid  squares  in  y  direction  of  a 

BLDE  unit  location  with  respect  to  the  lower  left 
grid  square 

12  =  relative  number  of  grid  squares  in  y  direction  of  a 

RED  unit  location  with  respect  to  the  lower  left 
grid  square 

13  =  number  of  grid  squares  in  y  direction  between  a  BLUE 

unit  and  a  RED  unit 
Z  =  temporary  variable  for  the  elevation  at  a  point  on  a 

hill  (Z) 
ZB  (I)  =  elevation  of  BLOE  unit  I 
ZR  (J)  =  elevation  of  RED  unit  J 
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SYNOPSIS  OF  THE  POSTPROCESSOR 

The  following  appendix  is  a  summary  of  the  work  done  by 
other  students  as  a  part  of  a  joint  project  for  OA  4655. 
The  authors  gratefully  acknowledge  the  work  done  by  Captains 
Ambrose  R.  Hock  and  Steven  L.  fladdox  in  developing  this 
postprocessor  to  reduce  the  amount  of  work  required  in 
extracting  the  necessary  data  from  STAR.   Their  work  is 
included  here  in  this  appendix  since  it  is  an  integral  part 
of  the  parameter  estimation  process. 

The  postprocessor  developed  for  this  thesis  is  a  Sim- 
script  language  program  that  provides  to  the  user  a  summary 
of  the  data  generated  from  the  Simulation  of  Tactical  Alter- 
native Responses  (STAR)  Model.   Additionally,  the  postproces- 
sor performs  some  routines  that  assist  the  user  in  the  data 
analysis  of  the  model.  The  STAR  Model  is  a  Simscript  lan- 
guage program  that  simulates  combat  between  two  combined 
arms  teams  in  a  combat  environment  that  includes  field 
artillery  fire  and  electronic  warfare. 

Since  the  model  assumed  in  the  project  only  required 
data  generated  by  a  tank  vs  tank  battle,  it  was  necessary  to 
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change   certain   aspects  of   the   STAB   model.      Thus   it  was 

necessary  to   "turn   off"   the   field  artillery   and   Electronic 

Warfare   modules   within  the  STAR   model.    Also    the  non-tank 

units    needed   to    be   deleted.      Referring   to  figure    19,    the 

Electonic   Warfare   code(EWCODE)    was   deleted   through  the  use 

of   a   comment   card,    and  a    Electonic   Warfare  off  (EWOFF) 

subroutine   was   inserted. 

//STARTHES    JOB     (3102, 0234) ,• COMB   ARMS    BATTLE •, CLASSIC 
//SIM.SYSPRINT    DD      DUMMY 

//SIM.SYSIN      DD      DISP=SHR, DSN=MSS.S3102. THESIS (WRKPREAM) 
//  DD      DISP=SHR, DSN=MSS.S3102. THESIS (SNAPDOTR) 

//*  DD      DISP  =  SHE,DSN=MSS.S3102. THESIS  (EWCODE) 

//  DD      DISP=SHR,DSN=MSS.S3102. THESIS (MOVCOORD) 

//  DD       DISP=SHR,DSN=MSS.S3102.THESIS(COMMWORK) 

//  DD      DISP=SHR,DSN=MSS.S3102. THESIS (TEMPCOMM) 

//  DD       DISP=SHR,DSN=MSS. S3102. THESIS (TEMPGRND) 

//  DD       DISP=SHH,DSN=MSS.S3102.THESIS(EWOFF) 

//  DD      DISP=SHR,DSN=MSS.S3102. THESIS (FAEXPER) 

//  DD      DISP=SHR,DSN=MSS.S3102. THESIS (ARTYMESG) 

Figure    19:      Electronic  Warfare   Deletion 

Thus,    whenever   the   main   program   called   for    the  usage   of  EW, 
the   EWOFF   subroutine   automatically   returned   the   program  bade 
to   the    point    where  the   EW   routine  was   called   and   the    simula- 
tion  continued    with  no   Electronic   Warfare  generated.      The 
field   artillery   within  the   STAR   model   was  cancelled   by 
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having    the   simulation   create   FA   fire   at    time   equal  to    5000 
time   units  (see    Figure    20). 

Schedule   a   FABEGIN  in   5000.0001    time  units 

Figure   20:      Field   Artillery   Deletion 

Since   the   simulation    was   scheduled    to   run  only   2500   time 
units,    no   field    artillery   was   generated.      Finally,   the 
internal  logic   of   the   STAR  model  was   used  to    "destroy"   the 
non-tank   units   within   the    simulation   by   changing    their   basic 
load   to   zero (see   Figure  21).      This   action  created  situa- 
tions,   where   the  tank    elements    did   not   consider   the   non-tanJc 
elements  as   lethal   or   dangerous   targets   and   thus    were 
ignored. 
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Figure   21:      Alteration   of   Basic   Load 

Output   from   the   STAB   model   can    be   divided   into  three 
categories: 
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Opiate  list  of  attacker/defender  status. 

Shot  data 

Significant  event  listing  (i.e.  movement,  EH  and  FA 

events) 
Only  the  shot  data  from  the  STAR  model  is  needed  to  provide 
the  necessay  data  for  computations.   Figure  22  is  a  sample 
of  this  output.   In  order  to  utilize  this  data,  the  internal 
output  control  of  the  model  was  used  to  place  the  shot  data 
in  a  mass  storage  location  within  the  computer  to  be  later 
manipulated  by  the  postprocessor. 

The  flow  chart  for  the  postprocessor  is  contained  in 
Figure  23;  a  copy  of  the  program  itself  can  be  found  in 
Appendix  F  as  Listing  I.   The  program  flow  is  very  basic  and 
needs  little  explanation.   From  the  mass  storage  location  of 
the  computer,  the  postprocessor  reads  pertinent  information 
into  the  SHOTLIST.   Once  all  of  the  data  has  been  read,  the 
program  processes  the  data  into  a  casualty  list (CASLIST) 
keying  from  a  status  of  "dead"  from  the  Shotlist.   Figure  24 
is  a  sample  of  the  casualty  list  output.   Once  CASLIST  is 
created  the  postprocessor  performs  histogram  computations 
and  graphs.  The  postprocessor  can  print  two  types  of 
histograms.   Figure  25  contains  the  Simcsript  Library  pro- 
gram for  the  histogram.  Additionally,  the  HISTG  routine  in 
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Figure  22:   STAR  Output  Shotlisr 
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FORTRAN  can  be  accessed  to  produce  a  histogram  of  the  data. 
After  the  histogram  data  is  generated,  the  postprocessor 
calculates  the  attrition  rate  coefficients  using  the  formu- 
las discussed  in  Section  III-E  of  this  thesis. 

Brief  description  of  the  postprocessor  (see  this  appen- 
dix for  variable  description) : 
Preamble : 

lines   1  -   9   definition  of  permanent  entities. 

lines   10  -  39   definition  of  temporary  entities. 

lines  40  -  54   definition  of  data  to  be  collected  for 

histograms. 
Main : 

line    2        variable  definition. 

line    3-4   system  parameters  read. 

lines   13  -  31   creation  of  shotlist,  reading  data  from 

file, 

lines  43  -  89   creation  of  casualty  list,  file  "dead". 

lines  90  -  125   processing  of  histogram  data 

lines  125  -  186   processing  of  attrition  rate  coeffi- 
cients. 
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VARIABLES, 


REfiiD 
DATA 


CREATE 
SHOTLlST 


NO 


INiriALliE 
DATA 


POfl    Ps.Oi 

TO  .9^  Si  .01 

DATA 


*   PROCESS 


Compute 

^}ST06RAM 


NO 


'igura    23:      Plow    Chart    of  -he   Postprocessor 
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DO 

CALCULATIONS 


CALCULATE 
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Print 
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SHOT    LIST    DATA    FOR    CASUALTIES   FOB    THIS    STAR    SIMULATION 


TIME    BETWEEN    CASUALTIES 
0 

12 

134 

106 

7 

10 


6 
28 
27 
38 
20 
98 


NOM 
1 

TIME 
322 

TARGET 
59 

SYS-WPN    TYPE 
1-7 

2 

334 

46 

1-7 

3 

46  8 

58 

1-7 

4 

57  4 

28 

1-1 

5 

581 

42 

1-7 

6 

591 

27 

1-1 

107 

1210 

72 

1-7 

108 

1238 

69 

1-7 

109 

126  5 

73 

1-7 

110 

130  3 

75 

1-7 

1  11 

1323 

19 

1-1 

112 

1421 

76 

1-7 

Figure   24:      Casualty   Lisx 
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NUMBER  OF    BLUE    KILLED    =  28 

BLMU   =         47.25000    BLSIG   =  299.259 

NUMBER  OF    RED    KILLED    =  84 

RDMU    =          13.24096    RDSIG   =  219,106 


HISTOGRAM   OF    BLUE    AND    RED    CASULTIES 
(DATA    BASED    ON    TIME    BETWEEN    BLUE    CASULTIES 
AND   TIME    BETWEEN   RED    CASULTIES) 

INTERVAL=10    UNITS  #BLUE    CAS       #RED   CAS 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 


Figure  25:   Simscript  Histogram 

Lisi  of  Variables  For  Postprocessor 
AHAT  =  estimate  of  the  attrition  coefficient  for  RED 

attriting  BLUE 
BHAT  =  estimate  of  the  attrition  coefficient  for  BLUE 

attriting  RED 
BLHST  =  Simscript  routine  for  the  histogram  of  BLUE  time 

between  casualties 
BLMU  =  Simscript  routine  for  calculating  the  mean  time 
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3 

11 

4 

3 

4 

20 

0 

4 

0 

8 

0 

4 

0 

0 

0 

8 

0 

0 

0 

0 

0 

4 

4 

8 

0 

4 

4 

4 

0 

0 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

5 

0 

between  casualties 
BLSHTHST  =  Simscript  routine  for  the  histogram  of  BLUE 

shot  times 
BLSHTIG  =  Simscript  routine  for  calculation  of  the 

standard  deviation  of  BLUE  shot  times 
BLSHTMU  =  Simscript  routine  for  calculation  of  the  mean 

of  BLUE  shot  times 
BLSIG  =  Simscript  routine  for  calculation  of  the 

standard  deviation  of  BLUE  time  between  casualties 
BLUE. SHOTS  =  time  of  the  previous  BLUE  shot 
BLUE. TIME  =  time  of  the  previous  BLUE  casualty 
CAS  =  casualty 
CASHST  =  Simscript  routine  for  the  histogram  of 

total  time  between  casualties 
CASIG  =  Simscript  routine  for  calculation  of  the 

standard  deviation  of  the  total  time  between 
casualties 
CASLIST  =  array  containing  list  of  casualties 
CASMU  =  Simscript  routine  for  calcualtion  of  the  mean 

for  the  total  time  between  casualties 
CBLUE  =  number  of  surviving  BLUE  force  elements 
CK  =  integer  used  in  the  calculation  of  the  attrition 
coefficient 
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1    =   BLUE   casualty 

0   =  RED   casualty 
CRED  =    number   of  surviving   RED   force   elements 
CORTIME   =   time   of   previous  casualty 
FNAME   =   firer's   name 
FSWTYPE   =   firer's   system-weapon   type 

J.BL  =    number   of   BLUE   casualties (used   in   BLHST   routine) 
J.RD  =    number   of   RED   casualties (used   in    RDHST   routine) 
JBLHST    =   number   of  BLUE   shots(used   in   BLSHTHST  routine) 
JCAS   =    number   of  casualties (used  in   CASHST   routine) 
JRDSHT   =   number   of   RED  shots (used  in   RDSHTHST  routine) 
JSHOT   =   number   of  shots(    used   in  SHOTHST   routine) 
JTOT  =    number   of  casualties (used  in   TOTHST   routine) 
LIM  =   array   of   time  periods  that  represent   phases  of 

the   battle 
LIMA  =    lower    time   limit   of   phase 
LIMB  =    upper   limit   of    phase 
LP  =   partial   derivative  of   p 
LQ   =   partial    derivative   of   q 
MAXBLDE    =    initial    BLUE   force   size 
MAXRED   =   initial   RED    force   size 
MINP  =    minimum    of   p 
MINQ   =    minimum    of   q 
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MK1   =   number   of   RED   elements   in   time   interval 

NK1    =    number   of    BLUE    elements   in  the, time  interval 

NLIM   =    number   of  time   intervals   to  be   evaluated 

NSHOT   =   number   of   shots  in   the  total   battle 

NTABLE   =   number   of   tables (battles)    to   be   evaluated 

NOM  =   counter    fo   number  of   casualties 

NUMBL.CAS   =   number   of   BLUE  casualties 

NOMRD.CAS   =   number   of    RED   casualties 

NUMCAS    =   casualty   entity    number 

P  =   probability  of   not   acqiring   a  target 

PHAT   =    estimated   p 

PROVAL   =   process   variable 

0  =  stop 

1  =  process 
QHAT  =  estimate  of  g 

RANGE  =  range  of  firer  to  target 

RDHST  =  Simscript  routine  for  the  histogram  of  RED 

casualties 
RDMU  =  Simscript  routine  for  calulation  of  the  average 

time  between  RED  casualties 
RDSHTHST  =  Simscript  routine  for  the  histgram  of  total 

red  shots 
RDSHTMU  =  Simscrip-L  routine  for  calcularion  of  the  average 
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rad  shot  time 
RDSHTSIG  =  Simscript  routine  for  the  calculation  of  the 

standard  deviation  of  the  RED  shot  time 
RDSIG  =  Simscript  routine  for  calculation  of  the 

standard  deviation  for  RED  time  between  casualties 
RED. SHOTS  =  time  of  previous  RED  shot 
RED. TIME  =  time  of  previous  RED  casualty 
RUNNO  =  run  number 
SHOT  =   shot 

SHOTHST  =  Simscript  routine  for  total  shots  fired 
SHOTLIST  =  array  of  shot  times 
SHOTMU  =  Simscript  routine  for  calculation  of  the  average 

time  of  shots 
SHOTSIG  =  Simscript  routine  for  calculation  of  the 

standard  deviation  of  total  shots 
STATUS  =  status  of  shot 
STATUS. CAS  =  status  of  casualty 
T.BTWN.CAS  =  time  between  casualty 
T.CAS  =  casualty  number 

TBLUE  =  number  of  BLUE  casualties  within  the  interval 
TGT.CAS  =  target  casualty 
TIME  =  time  of  shot 
TIMECAS  =  time  of  casualty 
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TIMELIM  =  maxtime  of  battle  when  a  casualty  occurred 

TMCAS  =  time  of  previous  casualty 

TMRDCAS  =  time  of  previous  RED  casualty 

TOTAL. CASTIME  =  total  casualty  time 

TOTAL. SHOT. TIME  =  total  shot  time 

TOTAL. TIME  =  total  casualty  time  for  the  creation  of 

T.BTWN.CAS 
TOTHST  =  Simscript  routine  for  the  histrogram  of 

total  casualty  time 
TOTMO  =  Simscript  routine  for  calculation  of  the 

average  casualty  time 
TOTSIG  =  Simscript  routine  for  calculation  of  the 

standard  deviation  of  the  total  casualty  time 
TRED  =  total  RED  casualties  within  the  time  interval 
TSUB  =  time  counter 
TSWTYPE  =  target  system-weapon  type 
TTBLUE  =  total  BLUE  casualties  for  XX  battles 
TTRED  =  total  RED  casualties  for  XX  battles 
TTCAS  =  total  casualties  for  XX  battles 
TTT  =  number  of  shot 

TYPE. CAS  =   system  weapon  type  of  the  casualty 
VAR  =  attrition  coefficient  variables 
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LISTING    OF    C0MP02SH    PHOGSAMS 


A.       DAMSTAC    MODEL    PSOGSAfl 

id     REM     IJPOPTED  24   FEB  :52 

29     REM     OETERMINISTIC  QGi^PEePTEH  HCfOEL  OF  3TQR  nn  THE 
21      REM     mPPLE  computer   v OPMSTflC ) 

38     DIM  )<Q(.  5  ::'/T'e<  5  >  .XRiT  5  ■>,VR<  5  >>eiX<  5.  l>^  >  .3V<  5,  t0>*RX<  5*  li?> 
40     0 1  M  BP<  5  > » RP<.  5  ::• .  S8<  5  > ,  :5R«r  5  > .  ft Vr  5 . 1  y  ^ 
50     DIM  8F<5:'>RF<5:',8T<5.5>,RT<5,5>>n8<'!i's0R<?> 
60     DIM  BR<5>,RR<i=;:' 
79     DIM  81';  ^v^ftl'Sii- 
SS  T  =  19:  DT  =  70:  MT  =    1508 
it&  OS  =     CHRS  '.'4) 
lw30  P  =  0,:B8:Q  =  0.83 
110  PI   =  3. 141582654: PR  =  0 
120  H  =  2.4:B  =   1.2 
130  BE  =  0.3: RE  =  0.3 
140  6M  =   .5:PM  =   ,5 
150  KK   =  0 
160  NK  =  MT  .••••  OT 
iSl      REM     ZZ=FLP8  FOR   TVPE  OUTPUT  ~ 

162  REM     ZZ=1   '3IUES  UNIT 

163  REM     LOCATIONS  :?«  SIZES 

164  REM     ZZ=0  GI'.JES   TIME  ix  TOTQL 

165  REM     RED  :?<  BLUE  FORCE  LEV 'ELS 

166  ZZ  =  0 

170  DIM  80<  NK  >  ,R0<  NK  ;:• 

180  HOME 

181  REM     READ  Ol^TQ  FROM  FILE   FOR 

182  REM     BLUE  FORCE. 

133  REM  OHTfi   INCLUDES: 

134  REM  M^NUMSER   OF   UNITS 

135  REM  81=NUMeER  OF  COORD  RTS 

136  REM  BX,8V=8RID   OF  COORD   PT 

137  REM 

130  INPUT  "NPHE  OF  FILE  =0R  8Li€  OQTQ-'-    "  :NS 

200  PRINT  OS; "OPEN  ";N«;".6LUE" 

2 10  PR I NT  OS i " READ  " ; N« ; " . BL UE " 

220  INPUT  M 

.i30  FOR  I  =  1  TO  M 

240  INPUT  Rl 

250  FOR  J  =  1  TO  Si 

260  INPUT  B:-v:  I  .J't  INPUT  eV':I,J':' 

270  NEXT  :  ^^E:<T 

230  PR  I  NT  us  i  "CLOSE  '* ;  HS  ;  "  ,  SLUE  " 
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281  REH  SET  INITIAL  UNIT 

282  REM  LuCPTTON  SQUfiL  TO 
233  REH   1ST  COORO  PT. 
2S4  REM 

2d€i  FOR  I  =  t  TO  M:xe<r:'  =  9X^' T  >  1.  > :  V8<' I  >  =  SVC  T ,  t  :•  s  RP<  T  >  =  1 
•TSi-J  PRINT  Bl:  FOR  JJ  =  1  TO  ei:  PRINT  BX-:' I  .JJ  :■  .ev-r  I  ,J..r>:  NEXT  JJ 
.?.ly  PRINT  "SIZE  OF  BLUE  UNIT  ".-Tj:  TNPMT  3B<  !• 

320  3U.  I  >  =  SB<  I  ;• 

321  REH 

322  REH  SUM  INITIfiL  FORCE  SIZE 

323  REM 

330  3w  y  >  =  6€K  0  :•  -i-  S8<  I  > 

331  REM 

332  REH  RPTE  OF  BLUE  MOUEHENT 

333  REH  EQUALS  y 

334  REH 

340  BR( I >  =  0 
350  NEXT 

360  PRINT 

361  REH 

T.62  REH  REW  OPTP  FROM  FILE  FOR 

363  REM  RED  FORCE. 

364  REH  SPHE  FORMOT  PS  BLUE 

365  REH 

370  Ih^PUT  "NPHE  OF  FILE  FOR  RED  OPTQ'?  "  :N« 

380  PRINT  OS; "OPEN  ";N;$;".RED" 

390  PRINT  OS;"RPPO  ";HS;",  RED'* 

400   INPUT  N 

410  FOR  I  =  1  TO  N 

420   INPUT  Rl 

■J-30  FOR  J  =  1  TO  Rl 

4419   INPUT  RX<  I,J>:  INPUT  RV<' UJ;:- 

+50  NEXT  :  NEXT 

460  PRINT  OS; "CLOSE  ";NS;".RED" 

470  FOR  I  =  I  TO  N:XR<I;'  =  RX<  I ,  !> :  VP*"  i  :•  =  R'-'^  I .  I  :':RP<  [  ;■•  =  I 

480  PRINT  "SIZE  OF  RED  UNIT  ";!;:  INPUT  SR<  T  > 

■+90  Rl<  I  J  =  SR<  I  :■ 

500  R0«:  0  >  =  R0<  0  :'  +  SR'-  I  ■■ 

510  RR<'  I  >  =  5 

520  NEXT  I 

530  PRINT  :  INPUT  "DO  VOU  MPNT  ri^E   PRINTER  ON'^  •  V.-h  ;,  "  :ps 

535   IF  PS  =  "V"  THEN  PR  =  1 

540  HHJME  :  UTP8  3 

550   IF  PR  =  1.  THEN  PR#  1 

560  PRINT 

564  REH 

565  REH  PRINT  INITIPL  OPTP 

566  9.B^ 

570  PRINT  "TIME";:  POKE  36 p7:  PRINT  "P».MF":;  POKE  36^20:  PRINT  "RED" 
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580  PPM   0 

584  REM 

585  REM  WHICH  OUTPUT  TVPE'?^ 
58b  REM 

590  IF  ZZ  =  1  THEN  r^jri^UB  17B«?  WTQ  R\iyf 

S^(d     i?iJSU6  1490 

SUi  REM  MP  IN  PROGRPH 

820  T  =  T  +  OT 

824  REM 

ri25  REM  CHECK  IF  MPX  TIME 

826  REM  13  EXCEEDED 

827  REM 

tioei   IF  T  <   =  MT  THEh^  840 

n35     PRINT  "END  OF  8PTTLE  DUE  TO  TJHE";  ROSUB  1540:  END 

640  80SU6  790:  REM  MOUEHENT  SUeROUTIHE 

850  I5CJSU8  1170-  REM  PTTRITION  SUPROI»TINE 

854  REM 

855  REM   INCREMENT  THE  TIME 

856  REM  PERIOD  INDEX 

857  REM 

860  KK  =  KK  ■»•  I 

^7^16   e0(KK>  =  0:R0O:)<:'  =  0 

874  REM 

875  REM  SUM  BLUE  FORCE  LE'-'EL 

876  REM  FOR  THIS  TIME  PERIOD 

877  REM 

$30  FOR  I  =  1  TO  H 

890  e0<  KK  :-  =  e0<  KK  >  ■»•  SB<  I  ':• 

?'00  NEXT  I 

70**  REM 

?05  REM  SUM  RED  FORCE  LE'.iEL 

706  REM  FOR  THIS  TIME  PERIOD 

707  REM 

710  FOR  J  =  1  TO  N 

721^  R0<  KK  )  =  R0C  KK  >  >  sp/  J  ;:• 

730  NEXT  J 

734  REM 

735  REM  WHICH  TVPF  OF  OUTPUT'^* 

736  REM 

740  IF  ZZ  =  1  THEN  80SUB  1760:  80Tn  ^i^m 
750  GOSUB  1  «iS0 
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752  REM 

r???.  rt£H   CH6CK  TF  ?H0  OP  BCiTTI  ? 

754  SEM  -  CR ITER IP  POR  PORCE 

755  REN   LEUEL3  13  EXCEEDED 

756  REM  IF  CRITERIm  13  EXCEEDED 

757  REM  I50T0  SPUE  RESULTS 

758  REM 

rsti  IF  B0<kk:>  >  be  *  B0<ei>  then  77c* 

765  PRINT  :  PRINT  "END  OF  BOTTLE  DUE  TO  PTTRTTIOH  OF  BLUE" 

766  80SUB  1540:  END 

770  IF  R0<:KK)  >  RE  «  R0<0:'  THEN  790 

775  PRINT  ;  PRINT  "EHO  OF  BPTTLE  OUF  TO  ATTRITION  OF  RED" 

776  30SU6  1540:  END 
730  i5*IiT0  620 

730  REM  MOUEMEh^T  SUBROUTINE 

734  REM 

.^35  REM  FOR  RED  UNITS 

736  REM 

300  FOR  I  s  I  TO  N 

804  REM 

;::05  REM  SET  LAST  COORD  PT 

806  REM 

310  J  =  RP<  I> 

313  REM 

314  REM  COMP^JTE  COMPONENTS  OF 

315  REM  OISTGNCE  BETWEEN  NEXT 

316  REM  COORD  PT  PNO  CURRENT 

317  REM  UXPTION 
313  REM 

320  DV  =  RV<  I ,  J  -K  I  >  -  VR<  I  :• 

330  OX  =  RX<  I, J  -c  1 )  -  XR<  I  > 

:533  REM 

334  REM  00  3E0METRV  TO  OET 

335  REM  CORRECT  PNGLE  '.PN::- 

336  REM 

340  IF  OX  =  0  MNO  OV  ••  0  THEN  PN  =  PT  •■•  "■;  kOTO  33»^ 

350  IF  OX  =  0  QNO  OV  ■(    0  THEN  QN  =  3  *  PT  -■  :"':  3nTn  330 

360  wN  =  HfN  'DV  -^  OX:' 

370  IF  OX  <  0  THEN  PN  =  PI  +  PN 

:-:73  REM 

374  REM  COMPUTE  OH,  OISTPNCE 

375  REM  TO  BE  TRPUELLED  IN 

376  REM  THIS  TIME  PERIOD 

377  REM 

3'30   01    =  OT   *  RR<  I  >   *     COS   •:  CiN  :< 

390  02  =  OT  *   RRU  :•   *      SIN   -PN':- 
'300  OH  =     SQR  <  0 1    -'^   2  -t-  02    '^   2  > 
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:5e4  REM 

305  REM  COMPUTE  01 »  0I3TPNCE 

306  REM  TO  ^€^<T  COORD  PT 

307  REM  FROM  CURRENT  LOCPTION 

308  REM 

310  01  =  3QR  <  OX  -^  3  +  OV  -^  2> 
911   REM 

312  REM  IF  01  <  OH,  MOUE  UNIT 

313  REM  ONLV  THE  mMOUNT  TO 

314  REM  MEXT  COORD  PT 

915  REM  IF  or  13  USED»  UPDPTE 

316  REM  LPST  COORD  PT 

917  REM 

320   IF  01  >  =  OH  THEN  350 

330  01  =  0'^:03  =  OV 

340  Rp<  I  >  =  J  -K  I 

350  VR<  r>  =   TNT  •:  VR<  I  >  *  02:- 

360  XR<  I  >  =  INT  <  ;<R(  !:'  >  01  > 

370  NEXT 

374  REM 

375  REM  00  THE  3wME  FOR  BLiJE 

376  REM 

380  FOR  I  =  1  TO  M 

385   IF  BR<I>  =  0  THEN  1150 

380  J  =  eP<  I  :• 

1 000  OV  =  8V<  I ,  J  ■♦•   1  :•  -  ve<  I  •> 

1010  OX  =  i5X(  I , J  -^   r>  -  :<B<  i::- 

1020   IF  OX  =  0  i^JNO  OV  ••  0  THEN  QN  =  PT  /  3:  HfiTn  \i^ifj,fi 

1030   IF  OX  =  0  WNO  OV  •:  0  THEN  QN  =  3  *  PT  -'  2:  80T0  1.060 

1040  QN  =  QTN  <0V  •''  OX':' 

1050'  IF  OX  <    0  THEN  PN  =  PI  +  PN 

1 060  01  =  OT  *  8R(  I  :>  *  COS  •:  mN  '> 

1070  02  =  OT  *  B9.<  I  ■<   *  3TN  ^PN) 

1080  OH  =  SQR  '-01  ^  2  +  02  --  2) 

1080  01  =  SQR  ':DX  -  2  +  OV  -^  2::' 

1100   IF  01  >  =  OH  THEN  1130 

1110  01  =  OX: 02  =  OV 

1120  BP<  I  :•  =  -J  >  1 

1130  VB<  I)  =   INT  <'t>B<1:>   +  02' 

1 140  xe(  I  '■'   =   INT  '  XB*  l'.<  ■*-   01; 

1150  NEXT 

1160  RETURN 


142 


n?0      REH        i^TTRITION  SUeROUTIN£ 

1172  REH 

1173  REH  ZERO  OUT  QTTRITTON 

1174  REH  UCiLUES  FOR  wLL  UNtTS 

1175  REH 

1180     FOR   I   =   1    TO  H:eF<I>  =  i??Oe<r:>  =  I-"'''    NEXT 

1190     FOR  J  =   1    TO  N!RF<J)  =  ia:OR<J'J  =  y:    NEXT 

1192     REH 

1133  REH  CHECK  RONijE  8ETWEEH 

1184  REH  OPPONENTS.  IF  RfiNi3E 

1135  REH  <=  2700  PTTRITTON  MILL 

1196  REH  OCCUR.  REDUCE  RPTE  OF 

1137  REH  RED  HOUEHENT  TO  2  H.^' 

1198  REH 

1200  FOR  I  =  1  TO  H 

1210  FOR  J  =  I  TO  N 

1220  8T<I^J>  =  0:RT<J.ir)  =  0 

1230  R6  =  3QR  <  <  XB<  I  >  -  ';<P:<  }>>   >  2  +  ''VP^r-  -  VRCJ>> 

1240   IF  Ri3  >  2700  THEN  1310 

1250  RPXj::-  3  2 

1251  REH 

1252  REH  INCREHENT  THE  NUHPER 

1253  REH  OF  T8T3  FIRED  QT  BV 

1254  REH  EPCH  UNIT,  COHPIJTE 

1255  REH  ATTRITION  CAUSED  8V 

1256  REH  RED  ON  BLUE  ^   MICE 

1257  REH  UERSA 

1258  REH 

1270  8F<  I)  =  BF<  !>  +  1 

1280  RF<J>  =  RF<J)  +  1 

1235   IF  RR^J)  =  -  9  THEN  BT-r  I  ,J  :-  =  0:  iViTO  1295 

1290  BT«.'I,J)  =  w  .-'  N  *  (1  -  (P)   -^  3Bt' I  >  )  *  SR':J-' 

1235      IF  9R<  r>  =  9  THEN  RT«r  ..M  >  =  0-    r^Ojn    nt0 

1300  RT-::j,I::'  =  8    •■'   H  *  <  1   -  <Q>  -^   SR<:  J  >  >  *  SB<  I  > 

1310  NEXT  :  NEXT 

1311  REH 

1312  REH  SUH  TOTAL  ATTRITION 

1313  REH  FOR  EACH  UNIT 

1314  REH 

1320  FOR  I  =  1  TO  H 

1330  FOR  J  =  1  TO  N 

1340   IF  RF«:j)  =  0  OR  rlF<  I)  =  0  TH^N  ^■^7<^ 

1350  OBCI  :.'  =  0B<  I  )  +  8T':I,J->  .•-  RF<;  J  ■> 

1390  OR<J.>  =  DR<J>  -(■  RT*'J,r:-  •-  9F»:r> 

1370  NEXT  :  NEXT 
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1371  9Ei-i 

1372  SEH  CQMPIJTF  HFM  mNTT 

1373  hEH  strength.  SUBTRACT 

1374  kEM  unit  lUTTRTTTi-iN  c^onM 

1375  REM  CURRENT  UNIT  STRENGTH 

1376  P.S}^      IF  UNIT  STRENGTH  i^lFS 

1377  REM  NEGOTIUE>  MPKE  IT 
1373  .  REM  EQUfiL  i9. 

1379  REM 

1380  FOR  I  =  1  TO  M 

1390  SBC  I  :>  =  SBC  r  :'  -  OBC  I  •- 

1400   IF  SB<  I  )  <  0  THEN  '=--R*' T  •'  =  i"* 

1410   IF  SBC  I >  <   =  BM  *  BIC I  )  THEN  SRC  I >  =  S 

1 420  h4£XT 

1430  FOR  J  =  1  TO  N 

1440  SRC  J)  =  SRCj>  -  DRCJ:- 

1450   IF  SRCJ>  <   0  THEN  SRC  J ':•  =  ^"i 

1460   IF  SRC  J)  <   =  RM  *  Pf  I'  THFN  9P<'.  .}>   ==  -  S 

1470  NEXT 

1480  RETURN 

1490  REM  PRINTOUT.'SUMMQRV 

1500   IF  PR  =  1  THEN  PR*  I 

1510  PRINT  T;:  POKE  36. 7 j  PRINT  qpi.fWk  •' ; 

1515  POKE  36,20:  PRINT  Rw3C  KK  ) 

1520  PRI^  0 

1530  RETURN 

1540  REM  SPUE  RESIJLTS  IN  Q  FT1.=? 

1550  PRINT 

1555   INPUT  "DO  VOU  MQNT  Q  FTLF  '"'F  tmtc  RIin^'-'  nqTp'?  -'V/N--  "  :P:$ 

1560   IF  PS  =  "N"  THEN  1750 

!57Q  OS   =  CHRS  -:4> 

1530   INPUT  "h«3ME  OF  FILE  FOR  FOprE  iPntri^  DPTQ'"^  "  ;F.S 

1590  PRINT  OSi'MiPEN  "  ;Ff ;  "  .COEFS" 

1600  PRINT  OS i "WRITE  "  ;FS;",r:OEF?" 

1610  PRINT  P:  PRINT  8:  PRINT  P:  PRINT  Q 

1620  PRINT  OS; "CLOSE  "  ;FS:'MTiFF<^" 

1630  PRINT  OS; "OPEN  " ;FS;", RESULTS" 

1640  PRINT  OS; "WRITE  " ;FS;". RESULTS" 

1650  PRINT  NK  +  1:  PRINT  m   +  1 

1660  FOR  I  =  0  TO  NK 

1670  T  =  I  *  .7 

1630  PRINT  T:  PRINT  90^) 

1690  NEXT  I 

1700  FOR  1=0  TO  t-^y: 

1710  T  =  I  *  3 

1720  PRINT  T:  PRINT  R0C  T- 

1730  Ne<T  I 

1740     PRINT  OS;"CLOSE   " ;F5;". RESULTS" 

1750      RETIjRN 
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1760  REM  PRINTOUT.-'SLIHMPRV 

1770  HOME  :  UTflB  3 

1730  IF  PR  =  1  THEN  PR#  I 

1790  PRINT  "TIHE  IS  "iT 

18€t0  PRINT  :  PRINT 

1310  PRINT  "BLUE  PORCES" 

1320  PRINT  "UNIT";  Tmi    7">i"XCnnR0" ;  TQe-r  t5>;"'T>cnnRD-;  TqB<  23  >; "SIZE" 

1330  FOR  I  =  1  TO  N 

1840  PRINT   "    ";I;   TPe<    7>;XB<  I  ■•;   TOR*'    15:';Ve<  I  :■;   Tpe<    24:.;SB<I> 

1350  NEXT 

I3fci0  PRINT  :  PRINT  "RED  FORCES" 

1 870  PR  I  NT  "  UN  IT";  TP6<  7  > ; "  NCOORD  "  ;  TPR-f  1  «=;  "> ; "  VCOORn  "  ;  TpB<  23  > ; "  S I ZE  " 

1330  FOR  I  =  I  TO  N 

1830  PRINT    "    ";I;    TPB<    7>;XR<I);   TPB<     t5:';VR<r>;    TPB<    24>;SR<I> 

1300  NEXT 

1310  PRINT  ;  PRINT 

1320  PR#  0 

1830  INPUT  "HIT  RETURN  TO  CONTINJJE" ;Zt 

1340  RETURN 
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B.       FOHCS    HAKER    PROGRAM 


Iti      REM     FORCE  OPTQ-FH.E  MQkER 

2vji     DIM  BXC  5*  U3 )  »BV<  5, li3  ::•  ,81':  19  > 

30     GOTO   179 

40     WDME 

50      INRfT   "NaHE   OF  FILS'^   "  ?N^ 

Stri  0*   =     CHRS   •:  4  ' 

TQ  PR  I  NT  US ; "  OPEN  "  ;  HS ',  " ,  "  ;  MS 
30  PRINT  OS; "WRITE  '*  iH^i'\'*  i^^S 
'i^     PRINT  H 

100  FOR  I  =  1  TO  M 

110  PRINT  ei<  I  :• 

120  FOR  J  =  1  TO  Bin::' 

130  PRINT  BXdiiJ):  PRINT  8Va,J> 

140  NEXT  J 

150  NEXT  I 

160  PRINT  OS  J  "CLOSE  '' iHS  i'\ '' i^MS 

170  HOHE 

130  UTwe  5 

190  PRINT  "MENU:" 

:l*00  PRINT  "1  MfiKE  FILE  FOR  BLUE  FOR»"ES" 

210  PRINT  "2  MfiKE  FILE  FOR  RED  FORCES" 

220  PRINT  "0  QUIT" 

230  PRINT  :  U^PUT  "WHICH';'  ";Z 

240  IF  Z  =  0  THEN  END 

250   IF  Z  =  1  THEN  MS  =  "BLUE" 

260   IF  Z  =  2  THEN  MS  =  "RED** 

270  HOHE 

230  PRINT  "WHQT  IS  THE  FORCE  SIZE  OF  ";MS;:  INPUT  "  ";H 

290  FOR  I  =  I  TO  M 

300  PRINT  :  PRINT  "FOR  UNIT  ":I 

310   PRINT  "HOW  MPNV  COORD  I NPT I NO  POINTS  "i 

315   INPUT  "<  INCLUOTNR  STPRT  LOCATIONS  V""  "jPKT' 

320  FOR  .J  =  1  TO  Bl':  I  > 

330  PRINT 

7.40   PRINT  "ENTER  C-fRID  COORDS  FOR  COORDTNClTTNi?  POINT  ^   " 

350   INPUT  "X-COORD  IS  ";8X':I.J-' 
360   INWJT  "V-<::nORD  IS  "  .-PV  I  ^...i  :• 

370  ne;<t  j 

330  NEXT  I 
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390  HONE 

40ti  PRINT   :    PRINT     Tm<     12 >;"***  OOTfi  CHECV 

4Ki  PRINT   :    PRINT   :    PRINT 

4c!y  PRINT   "UNIT  PT";   Tfte<    5  :';"^<-<:i:h"iPO'*  J   TQB<     tF^j'^T^-nXORO" 

430  FOR    I    =    1    TO   M 

440  FOR  J  =    1    TO  SKI  > 

450  PRINT     TP8<    2:)JIi   TpB<    fi^;Ji   Ti:ie<    10  >i:RX<:  T  j^.J  :-;   TP8(    13 ::';8V<  I,J ) 

460  NEXT  J 

•*70  NEXT    I 

*«0  INPUT  "ftNV  CHANGES?  (V/N>  '' iZS 

490  IF  ZS  =  "N"  THEN  40 

500  PRINT 

5 10  INPUT  "WHICH  UNIT-:-  ";I 

520  INPUT  "WHICH  POINT^^  "  ; J 

530  INPUT  "X-<:4D0R0  IS  '' iB^/ I  ,J) 

540  INPUT  "V-COORO  IS  'SBV^I,..!) 

550  GOTO  390 
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C.   RESULTS  BSADSa  PROGRAM 


iti  HUHE  :  UTwB  3 

20  0:5  =  ChRS   (4:. 

3y   INPUT  "NftHE  OF  FILE  POP  FORCE  '.E«'EL  DPTP'^  "  ;Ft 

40  PRINT  OS; "OPEN  "  ;FJ:;".i::OEFS" 

50  PRINT  0:$;"REflO  "  ;FS;'M":nEFS" 

80   INPUT  h:  input  8:  INPUT  P:  INPUT  Q 

P0  PRINT  OS; "CLOSE  " ;FS;". COEFS" 
;30  PRINT  DSi'MjPEN  "  ;FS;". RESULTS" 
30  PRINT  OS;"REPD  " ;FSi". RESULTS" 

100   INPUT  N:  INPUT  N 

110  NK  =  N  -  1 

120     DIM  80<NK>,P0'rNK> 

130     FOR    I    =  0  TO  NK 

140   INPUT  T:  INPUT  B0< I) 

150  NEXT  I 

160  FOR  I  =  0  TO  Nk; 

170   INPUT  T:  INPUT  R0(  I  > 

130  NEXT  I 

130  PRINT  OS; "CLOSE  " ;FS;" . RESULTS" 

200   INPUT  "PRINTER  ON"'  ";wS 

210  NOME  :  UTCiB  3 

220   IF  HS  =  "V"  THEN  PR#  1 

230  PRINT  "fl  =  "  :J4 

240  PRINT  "8  =  ";e 

i50  PRINT  "P  =  " ;P 

260  PRINT  "Q  =  ";Q 

270  PRINT 

ia0  PRINT  "TIHE"::  POKE  38,7?  PRINT  "SLUE"; 

2S5  POKE  3S,20:  PRINT  "RED" 

i90  FOR  KK  =  0  TO  NK 

3y^3  7  =  70  *  kk 

310  PRINT  T;:  POKE  38.7:  PRINT  80<KK:'; 

315  POKE  36,20:  PRINT  R0(:  KK  :- 

320  NEXT  KK 

330  PR#  0 
340  £N0 
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D,   MAP  KAK2H  PROGRAM 


Itt  REH   MPP  CPTRLOS  MPKER.-'REfiDEP 
i0  DIM  HL<  l^:1#l^^*l2> 
30  HOHE  :  MTPB  3 
4e  PRINT  "HENU:" 
50  PRINT  "1  MOP  MfiKER" 
60  PRINT  "2  MwP  READER" 
79  PRINT  "0  QUIT" 
30  PRINT  :  INPUT  "WHICH^  '*;X 
m      IF  X  =  0  THEN  END 
100  ON  :<  SOTO  120»;500 
110  REH   M«P  CPTRLOG  HPKER 
120  HOHE  :  UTfiB  2 

130  PRINT  T«8<  ll:-;'****  ^P   HPKER  **^" 
140  PRINT  :  PRINT 

150  INPUT  "NEH  OR  OLD  HPP-^  (N/O-  "  :ZS 
160  IF  ZS  =  "N"  THEN  II  =  1:JJ  =  1:  i50T0  220 
170  INPUT  "NPHE  OF  HPP^  ";H* 
1S0  QS   =  CHR:$  ';4j 
130  G03UB  620 

iQ0  INPUT  "ENTER  9TPRTIN»5  '5RID  'BOUppc  ''I.-J)  "z'^T.J.J 
210  ijOTO  250 

220   INPUT  "SIZE  OF  HPP  •:■  1000  M  SHUPRE? ''^  "jL^M 
230  INPUT  "ENTER  LOWER  LEFT  GRID  SQUPRE  'EG.  50,33::'  ";L:=:.LV 
240   INPUT  "BP3E  ELEUPTIDN  OF  HPP-^  "  :HL^:  0,0,0) 
250  FOR  I  =  1 1  TO  L 
260  FOR  J  =  JJ  TO  W 
270  PRINT  :  PRINT 

230  PRINT  "FOR  GRID  3QUPRE  %-I;".";J 

290   INPUT  "NUMBER  OF  HILL  MPSSES'^  "  ;HL'::  I  ,J,i-^  >:KK:  =  HLa,J,0; 
300   IF  KK  =   0  THEN  380 

310  PRINT  "ENTER  THE  '';KKi"   HILLS  FOR  nRID  '=:nuPRE  ";!>","  ;J 
320  FOR  K.  =  1  TO  KK 
J0  PRINT  "HILL  ";K;"  IS  ";:  INPUT  H!..-:  I  ,.J,K  > 


■-'•^r 


3'*0  ne;st  k 

350  JJ  =  1 

360  PRINT  :  INPUT  "SPUE'^  .;  V.-N  :-"  ; PS 

370   IF  w$  =  "V"  THEN  400 

380  NEXT  J 

330  NEXT  I 
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4W0  HOME  :  'JTP8  4 

4U3  INPUT  "HPME  OF  M«P  FILt'^  "  ;Hf 

42w  OS  =  CHRS  <4> 

430  PRINT  OS; "OPEN  '•;HS;",MfiP" 

4aiEt  PRINT  DSi "WRITE  *';HS;".MOP- 

45v3  PRINT  L:  PRINT  W 

4b«  PRINT  U<!  PRINT  LV 

479  PRINT  HL<d,id,Q> 

4yei  FOR  I  =  1  TO  L 

490  FOR  J  =  1  TO  M 

500  PRINT  HL< I,J,0> 

510  IF  HL':.I,J,0>  =  0  THEN  550 

520  FOR  K  =  1  TO  HLa*J.0> 

530  PRINT  HLa,J,K> 

540  NEXT  K 

550  NEiT  J 

560  NEXT  I 

570  PRINT  DSi "CLOSE  ";MS;".NPP" 

580  PRINT  :  PRINT  :  PRINT  :  PRINT  "...  DONE":  PRINT  :  ^RINT 

590  INPUT  "DO  you  WfiNT  TO  CONTINUE  INPUTTING  npTP^'  <  V.  N  :•  "  ;ZS 

^^ibid  IF  ZS  =  "V"  THEN  200 

810  iirtliTO  30 

620  0-£   =  CHRS  (4;.. 

1^30  PRINT  OS; -OPEN  ";MS;".MPP" 

640  PRINT  OS;"REPO  ";MS;".NPP" 

650  INPUT  L:  INPUT  W 

660  INPUT  l:<;  input  LV 

670  INPUT  HL<0i.0>0> 

660  FOR  I  =  1  TO  L 

690  FOR  J  =  1  TO  W 

700  INPUT  HL'.:i,J,0::' 

710  IF  HL':.I>J,0>  =  0  THEN  750 

720  FOR  K  =  1  TO  HL(I*J*0' 

:'30  INPUT  HLCI^JuKJ 

740  NEXT  K 

."50  NEXT  J 

760  NEv^T  r 

770  PRINT  OS; "CLOSE  "  ;HS;".^'H5P" 

760  RETURN 

790  REH  liPP  CPTOLOG  FILE  PEPDER 

800  HONE  :  '-'TPe  2 

610  PRINT  TPB<  11  :.;"***  ^^pp  OFQDER  ^^^^" 

620  PRINT  :  PRINT 
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340 

860 
::-i7v3i 
8:50 

800 

810 
920 
830 
840 
.-»50 
860 
870 

890 

1  000 
1010 


FILE' 
<  y.-'N 


1030 
1040 
1050 

1  080 
1070 

1080 

1090 

1  1 00 
1110 

1120 
U30 

1140 
1150 
1160 
1170 
1180 
1180 

I  i00 

1210 


PRINT 
INPUT 
INPUT 
INPIJT 
FOR  I 
FOR  J 
INPIJT 
IF  HL 
FOR  K 
INPUT 
NEXT 
NEXT  J 
NEXT  I 

PRINT 

HOME 

IF  QS 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

FOR  I 

FOR  J 

PRINT 

IF  HL< I 

FOR  K  = 


";MS;".H»:(P" 

M 
LV 


INPUT  "NPME  OF 
INPUT  "PRINTER 
OS  =  CHRS  <4:- 
PRINT  DS;"npEN 

0:lf;"REPD 

L:  INPUT 

LX:  INPUT 

HU  0,0*0  > 

=  1  TO  L 

=  1  TO  W 

HL<I,J,0) 
( I,J,0>  =  0  THEN 

=  1  TO  HL< I*J,0 

HU  I,J,K) 
K 


";P* 


860 


0*;'*CLOSE  '•;N«;",HfiP" 

=  "V"  ThE\^^      PR#  1 
'•NAME  OF  HPP:  ";MS 
"LOWER  LEFT  OR  ID  SQUPRE:  "  ;l.. 
"SIZE  OF  HPP:  ";Li"  BV  ";W 
ELE'-'PTION:  ";HL(  0.0,0 


;LV 


=  1  TO  L 

=  1  TO  N 

I^JiiHU:  I,J,0> 
.J,0)  =  0  THEN 
1  TO  HL< I,J,0 


1160 


PRINT  HL 
NEXT  K 
PRINT 
NEXT  J 
PRINT  : 
PR#  0 
PRINT  : 
INPUT 
bOTO  30 


IfJsK 


NEXT  I 


PRINT 
HIT  RETURN 


TO  CONTINUE"  :><'S 
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E.       HILL    aAKER    PROGHAa 


ly     REfl    -     HILL  Onto   FTL£  MPKEP-'REflPEP 

20     DIM  H<200*;5),K':  8> 

30  K<  1 :.'  =  5:K<2>  =  3:l«3>  =   14!K<4::'  =   Iw 

40  K<5>  =  22:K(6>  =  24:K<7>  =  29:K<3>  =  33 

58     REM  HILL  DPTO  FILE  MPKER 

68     HONE   :    UTPB  2 

79     PRINT   "MENU;" 

i50     PRINT   "l   HILL  MRKER" 

38     PRINT    "2  HILL  REPOER" 

188      PRINT    "8  QUIT" 

118     PRINT   ;    II^^PUT    "WHICH'?   "  ;;< 

128      IF  •;<  =  8  THEN     ENO 

138     ON  X  I50TO   150>320 

140     REM  HILL  OOTfi  FILE  MPKER 

150     HOME   :    UTPB  2 

168  PRINT  TP8<  11>;"^>h**  hILL  MPKER  *^h*." 

178  PRINT  :  PRINT 

138   INPUT  "NEW  OR  OLD  FILE?  •■' N.-O  -  "  ;P? 

138   IF  HS  =  "N"  THEK  238 

288  iSOSUB  560 

1'18   INPUT  "WHICH  HILL  GO  VOU  WPWT  tq  'JTpRT  WTTH'^-  ";II 

228  GOTO  268 

238  WUME  :  UTPB  4 

240  II  =  1 

258  INPUT  "NUMBER  OF  HILLS'^'  ";N 

268  FOR  I  =  11  TO  N 

278  PRINT  "ENTER  THE  DPTP  FOR  HILL  ";I 

2S8   INPUT  "'COORD  13  ";H':.  I,l> 

298   INPUT  "VCOORD  13  ";H<I,2> 

388   INPUT  "PEPK  ELEUPTION  IS  ";H<I,3':' 

318   INPUT  "OR  I ENTPT I  ON  PNGLE  IS  ";H''T.4:' 

328   I NPUT  "  ECCENTR I C I  TV  IS  "  ;  H<  1 , 5  :■ 

338   INPUT  "SPREPO  IS  ";H<I,6:' 

348   INPUT  "MPX  HEIGHT  IS  ";H':'I,7) 

358   INPUT  "CUT  HEIGHT  IS  ";H<I>3> 

368  PRINT  :  PRINT  :  INPUT  "SPUE^  (V.'-N)  ";Z$ 

378   IF  ZS  =  "V"  THEN  338 

338  NEXT  I 

338  HOME  ;  UTPB  4 

400  OS  =  CHRS  ':4::' 

*18   INPUT  "NPME  OF  FILE^  ";NS 

420  PRINT  OS; "OPEN  " ;NS;". HILLS,  L34" 

438  PRINT  OS;"  WRITE  " ;NS;". HILLS,  R0" 

440  PRINT  N 

•+58  FOR  T  =  1  TO  N 

4-68  PRINT  OSi  "WRITE  "  ;NS;"  .HILLS,  P"tI 
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470  FOR  J  =  1  TO  8 

48W  PRINT  H<  I>j> 

490  NEXT  J 

5WW  NEXT  I 

510  PRINT  OS;"CLOSE  " ;NS;". HILLS" 

520  PRINT  :  PRINT  :  PRINT  "...  DONE" 

525  PRINT  :  PRINT 

530  INRJT  "DO  VOU  MPNT  TO  CONTTNUE  I NPI  ITT TN6  DPTR-?-  <V.''N:--  "  jPS 

540  IF  |3;S  =  "V"  THEN  210 

550  i5«IlT0  60 

560  HOME 

570  OS   =  CHR«  •:4) 

580  INPUT  "NPHE  OF  FILE"^  "  ;NS 

530  PRINT  OS;  "OPEN  "  ;NS;",HTLL3.  LI-*" 

800  PRINT  0$;"REPD  ";N«;". HILLS,  R0" 

810  INPUT  N 

820  FOR  I  =  1  TO  N 

830  PRINT  D«i"REPO  ";N«,-". HILLS,  R":T 

840  FCiR  J  =  1  TO  8 

850  INPUT  H<  I,J) 

660  NEXT  J 

870  NEXT  I 

830  PRINT  OS; "CLOSE  " ;NSi". HILLS" 

890  HOHE 

r00  PRINT   "NUMBER  OF  HILLS:    "  ;N:    P-RINT 

710  M  =    I 

r20  FOR    I    =   I    TO   N 

730  IF  M  =  11  THEN  PRINT  :M  =  1 

740  PRINT  li 

750  FOR  J  =  1  TO  8 

780  PR  INT      Tf4e<    K  <  J  >  >  ;  H<  I  , .  J  :• ; 

770  ne^=:t  .J 

730  PRINT 
790  M  =  M  -t-  t 
800  NEXT  I 
810  RETURN 
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32id     REH     HILL  OftTP  FILE  REPDER 

330  PR  =  y 

340     HCiME   :    UTPIB  2 

;550     PRINT     T«B<    10:.;"***  HILL  REPDER  ***" 

360     PRINT   :    PRINT 

;i70  0*  =     CHft$  '■4-> 

380   INPUT  "NPHE  OF  MPP*^  '^N* 

y:30  INPUT  "PRINTER  ON'?  ''V.''N>  "iZ^^ 

900   IF  ZS  =  "V"  THEN  PR  =  1 

310  PRINT  0$,-"OPEN  ";N«;". HILLS.  L34" 

:320  PRINT  OS;"REPO  "  ;N«;". HILLS,  R0" 

330  INPUT  N 

340  FOR  T  =  1  TO  N 

350  PRINT  0*i"REPO  " ;N«;". HILLS,  R";I 

:360  FOR  J  =  1  TO  :=: 

;:»70   INPUT  H<I,J> 

3S0  H£:-<T  J 

bSki     NEXT  I 

1W00  PRINT  OS;"CLOSE  '•  ;N«;". HILLS" 

1010  HOME 

IV320  IF  PR  =  1  THEN  PR#  1 

1030  PRINT  :  PRINT  NS;"  HILL  DflTQ":  PRINT 

1040  PRINT  "NUhSER  OF  HILLS:  ";N 

1050  PRINT 

1060  H  =  1 

1070  FOR  I  =  1  TO  W 

1080  IF  H  =  11  THEN  PRINT  :M  =  1 

1090  PRINT  I; 

1100  FOR  J  =  1  TO  3 

1110  KK  =  K-rj:.'  +  3  *  J 

1120  IF  PR  =  1  THEN  POKE  "^R.KK:  PRINT  H^  T  , j -. ; •  GOTH  tl40 

1 1 30  PR  I  NT  TPB<  K<.  J  )  .:■ ;  H(  I  ,  ...I ) ; 

1140  r^iEXT  J 

U50  PRINT 

1160  H  =  H  *  1 

1170  NEXT  I 

1130  PR#  0 

1130  INPUT  "HIT  RETURN  TO  CONTTNUE"  jX'-S 

1200  8uT0  60 


154 


p.       300DS    MAKER    PROGRAi!! 


IW     REM       WOiJOS  CPTflLOe  MfiKER'-REPOER 

20     DIM  fU:  IS*  113*3 > 

30     HONE   :    UTP8  3 

40     PRINT   "MENU:" 

50     PRINT   "1   WOODS  MfiKER" 

60     PRINT   "2  WOODS  READER" 

re     PRINT    "0  QUIT" 

30  PRINT  :  INP<JT  "WHICH*^  "iK 

m      if   ;•<  =  0  THEN  END 

100     ON  X  60T0   120*770 

110     REM     WOODS  CflTQLDi?  MPIKER 

120     HOME    :    '.'TfiB  2 

130     PRINT     TfiB<     11>;"^H«.  WOODS  h+flKER   ^^^" 

140  PRINT  :  PRINT 

150   INPUT  "NEW  OR  OLD  MfiP-^  <  N-'-O'^  "  :?S 

160   IF  ZS  =  "N"  THEN  IT  =  1:JJ  =  1:  GOTO  220 

170   INP«JT  "Ni:;ME  OF  MOP?  ";M$ 

130  0£   =  CHRS  ';.4> 

190  ijOSUB  600 

i00  INPUT  "ENTER  STfiRTTNG  GRID  SQUQRE  <  I  .J  :•  ";II*JJ 

210  GOTO  240 

220   INPUT  "SIZE  OF  MfiP  <  1900  M  SOUPRFSi'-^  ";L..W 

230   INPUT  "ENTER  LOWER  LEFT  GRID  SQUARE  •:  EG.  50*33)  "*-LX*LV 

240  FOR  I  =  II  TO  L 

250  FOR  .J  =  JJ  TO  W 

260  PRINT  :  PRINT 

270  PRINT  "FOR  L-rRID  3QUPRE  ";r;"*";J 

280   INPUT  "NUMBER  OF  FORESTS'^  "  ;FLa  *J*0>:KK  =  FL':i*J*y> 

290   IF  KK  =  0  THEN  370 

300  PRINT  "ENTER  THE  ";KrK;"  FnPF'=;T'^  POP  i^RlD  SOUORE  ";!;"*";- 

310  FOR  K  =  1  TO  KK 

320  PRINT  "FOREST  "*K*"  IS  ";:  INPUT  FLC  I  *.J*K' > 

330  NEXT  K- 

340  JJ  =  X 

350  PRINT  :  INPUT  "SOUE^  ( V.-'N;-'";p-^ 

360   IF  A;$  =  "V"  THEN  390 

?.70  NEKT  J 

3S0  NEXT  I 
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.?.30  HOME   :    UTPtB  4 

4I30  INPUT   "NfiHE  OF  MPP  FILE?    "  ;M« 

41ti  OS  =     CHRS  <4) 

420  PRINT  OS; "OPEN  ";H$;". WOODS" 

+30  PRINT  OS; ''WRITE  "; MS; ".WOODS'* 

440  PRINT  L:  PRINT  W 

450  PRINT  LA:    PRINT  LV 

460  FOR  I  =  1  TO  L 

470  FOR  J  =  1  TO  W 

460  PRINT  FL<I,J,0::' 

490  IF  FL<I*J*0>  =  0  THEN  5>3 

500  FOR  K  =  1  TO  FUI/iJ,0> 

510  PRINT  FU  l,J,K) 

520  ne:<t  k 

530  ^♦£;^T  J 

540  NElsT  r 

550  PRINT  OS; "CLOSE  "; MS? ".WOODS" 

560  PRINT  :  PRINT  :  PRINT  :  PRINT  "...  DONE":  PRINT  :  PRINT 

570  INPUT  "DO  VOU  WPNT  TO  CONTINUE  TNPUTTIMG  OPTQ'?  -'V.-N:-  '' iZS 

580  IF  ZS   =  "V"  THEN  200 

590  GOTO  30 

800  OS  =  CHRS  ':4::' 

Hi  10  PRINT  OS; "OPEN  "; MS; ".WOODS" 

820  PRINT  OS; "READ  "; MS; ".WOODS" 

830  INPUT  L:  INP^JT  W 

840  INPUT  LX:  INPUT  LV 

850  FOR  I  =  1  TO  L 

860  FOR  J  =  1  TO  W 

870  INPUT  FL<  I,J,0:' 

860  IF  FL(I*J*0:>  =  0  THEN  720 

880  FOR  K  =  1  TO  FL';I.J,0.:' 

700  INPUT  FU:  I*J*K> 

ri0  NEXT  K 

720  NEXT  J 

r30  NEXT  I 

r4vj  PRINT  OS; "CLOSE  "; MS; ",  HOODS" 

750  RETURN 
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I  I 


6iEi  i^EM  WOODS  CPTPLOG  FILE  REPOER 

0  HOME  :  UTkB  2 

7m  PRINT  Tfie<  11  >;"***  W00n*5  REQOEP  ^**" 

790  PRINT  :  PRINT 

i^m  INPUT  "NfiME  OF  MPP  FILE?  " -.H« 

310  INPUT  "PRINTER  ON'^'  <V.--N>  "  ;P*- 
320  OS   =  CHRS  <4.:> 

330  PRINT  G$;"OPeN  " iM$;".NOOOS" 

340  PRINT  OS;"REPO  " ;MS;". WOODS" 

350  INPUT  L:  INPUT  W 

360  INPUT  LX:  INPUT  LV 

370  FOR  I  =  1  TO  L 

330  FOR  J  =  1  TO  W 

390  INPUT  FL< I,J,0> 
300   IF  FL<I*J,0>  =  0  THEN  :^40 

910  FOR  K  =  1  TO  FLa*J,0> 

320    input  fl<i,j,}<> 
330    ne:<t  k 

340  NEaT  J 

350  NEXT  I 

%0  PRINT  0«i "CLOSE  '*; MS; ".WOODS" 

370   HOME 

380   IF  PS  s  "V"  THEN  PR#  1 

390  PRINT  "NPME  OF  MPP:  ";MS 

1000  PRINT  "LOWER  LEFT  f5RID  SQUPRE;  "iLX.;" 

1010  PRINT  "SIZE  OF  MPP:  '^L;"  BV  " ;W 

1020  PRINT 

1030  FOR  I  =  1  TO  L 

1040  FOR  J  =  1  TO  W 

1050  PR  INT  I ,  J  ,FL':  I  ,J  *0  ■' 

1080   IF  FL<I*J*0>  =  0  THEN  1110 

1070  FOR  K  =  1  TO  FLa,J*0) 

1 080  PR I NT  FL< I > J ,K  ) ; "  " ; 

1090   NEXT  K 

1100   PRINT 

1110  NEXT  J 

1120  PRINT  :  NEXT  I 

1130  PR#  0 

1140  PRINT  :  PRINT 

1150   INPUT  "HIT  RETURN  TO  CONTINUE" : XT 

1130   i3uTG  30 
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G.       FOREST    MAKER    PROGRAM 


ly     REM     FOREST  MflKER^-REPDPP 

20     DIM  F<  12iii#i5>»K<i5> 

•30  K<1>  =  5:K<2>  =   12:K<3">  =  2^:K<4>  -  2'=; 

40  K(5>  =  30:K<9>   =  38 

50     HOME   :    UTfie  2 

80     PRINT   "MENU:    " 

r0     PRINT   "1    FOREST  MPKER" 

80     PRINT   "2  FOREST  REPOER" 

30     PRINT   "0  QUIT" 

100  PRINT  :  INPUT  "WHICH"^  "iX 

110   IF  X  =  0  THEN  END 

120  ON  ;<  OOTO  140,800 

130  REM   FOREST  DPTO  FILE  MCiKER 

140  HOME  :  UTfie  2 

I 50  PRINT   rP8<  10 >;"***  FOREST  HQKER 

180   PRINT  :  PRINT 

170   INPUT  "NEW  OR  OLD  FILE*?  -'N^O;'  " ;QS 

180   IF  Qs  =  "N"  THEN  220 

190  G08U8  530 

i00   INPUT  "WHICH  FOREST  DO  VOU  WPNT  TO  -^iTQRT  WI^H"'  ";[! 

210  GOTO  250 

220     HCiME    :    UTfiB  4 

230   11=1 

240   INPUT  "NUMBER  OF  FORESTS'?  "  ;N 

250  FOR  I  =  II  TO  N 

ZS^     PRINT  "ENTER  THE  DPTP  FOR  FOREST  ";T 

270   INPUT  ":<!:OORD  IS  ";F<I,r:' 

230   INPUT  "VCOORD  IS  ";F.'I,2:' 

230   INPUT  "HEIGHT  OF  TREES  IS  ";F<I,3> 

300   INPUT  "ORIENTRTTON  QNGLE  IS  ";F''I.4:. 

310   INPUT  "SEMI-HPJOR  mXIS  LENGTH  ts  ";F';:i,5> 

320   INPUT  "SEMI -MI NOR  fIXIS  LENGTH  IS  "  :F«:- 1 .8  :• 

330  PRINT  :  PRINT  :  INPUT  "ShUE'?  'V.-n;:'  ";ZS 

340   IF  ZS   =  "V"  THEN  380 

350  ne:^<t  I 

360  HOME  :  UTPB  4 
370  OS  =  CHRS  -^i' 


158 


ISid      INPUT  "NCiMP  OF  PTLE'^'  '*  ;NS 

390  PRINT  OS i "OPEN  " ;N«;". FORESTS,  L3S" 

480  PRINT  OS;"  WRITE  "  ;NSi'*  .FORESTS.  Rfr 

410  PRINT  N 

4«i0  FOR  I  =  I  TO  N 

+30  PRINT  D«; "WRITE  ";NS;". FORESTS,  R";I 

440  FOR  J  =  1  TO  S 

450  PRINT  F<  I^J:- 

4»50  NEXT  J 

4?0  hEXT   I 

*y0  PRINT  OS; "CLOSE  " ;NS;'\ FORESTS" 

490  PRINT  :  PRINT  :  PRINT  "...  DONE" 

495  PRINT  :  PRINT 

500   INPUT  "00  you  WPNT  TO  CONTINUE  TNPUTTIN»5  HhTP-^  (V/N:-  ";PS 

510   IF  QS  =  "V"  THEN  280 

520  i50T0  50 

530  HOHE 

540  OS  =  CHRS  <4:. 

550   INPUT  "\4QhS   OF  FILE"^  '' i^£ 

560  PRINT  OS;: "OPEN  ";NS;", FORESTS,  LIS" 

570  PRINT  OS; "READ  " ;NS;", FORESTS,  R0" 

530  INPUT  N 

590  FOR  I  =  1  TO  N 

t?00  PRINT  DS;"REPO  ";NS;". FORESTS,  R";I 

810  FOR  J  =  1  TO  S 

b20  INPUT  F<  I,j:- 

630  NE-'T  J 

640  NEXT  I 

S50  PRINT  OS; "CLOSE  " ;NSj". FORESTS" 
660  HOHE 

670  PRINT  NS;"  FORESTS  LIST":  PRINT 

680  PRINT  "NUHSER  OF  FORESTS:  ";N:  PRINT 

690  N  =  I 

ry0  FOR  I  =  1  TO  N 

710   IF  H  =  11  THEN  PRINT  :M  =  t. 

720  PRINT  I; 

730  FOR  J  =  1  TO  6 

740  PR  I  NT  Tfi8<  K(  J  :•  >  ;F<  I  ,  j  :• ; 

750  NEXT  J 

:'60  PRINT 

770  N  =  H  +  1 

790  NEXT  I 

790  RETURN 
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31  y 

330 
34V3 

330 
;590 

900 

310 


FOREST  DfiTP 

i.JTfiB  2 


REM 
HOME 
PR  =  0 
PRINT  TPe<  l©);" 
PRINT  !  PRINT 
OS   =  CHRS  <4:' 
INPUT  "NPME  OF  MfiP-^ 
"PRINTER  ON'T^ 
=  "V"  THEN  PR  = 


FILE  PFPDER 


FORFST  REPDER 


330 
;340 
350 
960 
370 
380 
990 

1000 
1010 
1020 
1030 
1040 
1050 

1060 

1070 
1080 
1090 

1  1 00 
1110 

U20 

1 130 
1140 
1150 
1160 
1170 
1180 


INPUT 

IF  ZS 

PRINT 

PRINT 

INPUT 

FOR  I 

PRINT 

FOR  J 

INPUT 

NEXT  J 

NEXT  I 

PRINT  0«;"CLOSE 

HOME 

=  1  THEN 
:  PRINT 


";N$ 

-:v.'n;> 


;Z$ 


1 


OS; "OPEN 

0S;"REfiO 

N 

=  1  TO  N 

0«i'*REflD 

=  1  TO  6 

F<I,J> 


;NS;".  FOREST??, 
;NS;". FORESTS, 


L3P" 
R0" 


;N«;". FORESTS,  R";I 


;NS;'*.  FORESTS' 


IF  PR 

PRINT 

PRINT 

PRINT 
M  =  1 

FCfR  I  =  1 

IF  M  »  11 

PRINT  I; 

FOR  J  s  1 
KK  s  K<J> 


NUMBER 


PR#  1 
mi"   FOREST 
OF  FORESTS  ' 


D«TQ": 
';N 


PRINT 


TO  N 
THEh^ 


PRINT  :M  =  1 


TO  6 


*  J 


POKE  36, KK: 
'  >;F<  I  ,J  ); 


PRINT  F.'I 


IF  PR  =  1  THEN 
PRINT  TfiB<  KCJ 
NEXT  J 

PRINT 
M  =  H  +  1 

ne:.<t  I 

PR#  0 

INPUT    "HIT   RETURN  TO   nONTINUF"  ;X-? 

13OTO  50 


GOTO    tl20 
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H.       T2RHAIN    MODEL    PROGRAM 


Itt     REH     UPOPTED  i^l2   P5B  32 
20     REH     DnHSTPC   TERRhIH  MODEL 

30     DIM  :<B<  5 '>,':>B(  5  :',XR<  5  >  .VR<  5  ■',B^<  S,  1 1^  •,  .ev<  5,  tPf  >  ,RX(  5*  IS  > 
I  35     DIM  ZB(  5  >  » ZR<  5  > , HL v 1 0  *  I y  *  1 2  > 
443     0 1 M  eP<  5  >  pPP<  5  "•• , S8<  5  >  > SR<  5  •• . ft V«"^  *  t  >?»  > 
50     DIM  BF<  5 >,RF<  5  >,eT<  5,5 ::',RT<  5,5  J^OBC  5 ■),0R<  5  :• 
80     DIM  SR<5>,ftR<5-' 
70     DIM  B1<5>,R1<5:' 
;i*0  T  =  0:DT  =  30:  MT  =   1500 
90  OS  =     CHRS  <4> 

100  P  =  0.9S:i3  =  0.93 

110  PI  =  3. 141592654: PR  =  0 

120  w  =  2.4:8  =  0.03 

130  BE  =  0.3: RE  =  0.3 

140  8M  =  .5:ftM  =  .5 

150  K,K  =  0 

160  NK  =  MT  /  DT 

165  ZZ  =  1 

170  DIM  B0C  NK  >  »ft0<  NK  :- 

171  HOME 

172  INPUT  "NPME  OF  MAP  FILE'?'  ";MS 

173  PRINT  DSi'MjPEH  ";MS;".MfiP" 

174  PRINT  O^i^REPO  ";MS;".HPP'» 

175  INPUT  L:  INPUT  W 

176  INPUT  l;<:  INPUT  LV 

177  INPUT  HU:  0,0,0:' 
179  FOR  I  =  1  TO  L 
179  FOR  J  =  1  TO  W 
190   INPUT  ML';  I,J,0> 
131   IF  HL'::i,J,0)  =  0  THEN  194 
192  FOR  K  ~   1  TO  HL-:  I  ,J,0;:' 
193   INPUT  ML':.  I,.J,K.> 
184  NEXT  :  NEXT  :  NEXT 

195  PRINT  OS; "CLOSE  "  ,-MS;",  MfiP" 

196  LX  =  LX  -  1:LV  =  LV  -  t 

197  PRINT  OS; "OPEN  " ;HS;". HILLS,  L34" 

198  PRINT  OS;"REflD  "  :MS;". HILLS.  R»^" 

199  INPUT  NH:  DIM  H'rNH,9) 

190  FOR  I  =  1  TO  ^4H 

191  PRINT  OS;"REPO  " ;MS; ". HILLS,  ft"; I 

192  FOR  J  =  1  TO  9 

193  INPUT  H< I,J> 

194  ^€XT  :  NEXT 

195  PRINT  OS; "CLOSE  " ;MS;",HTLLS" 
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I8w  HuhE 

iSS      INPUT  "NPME  OF  FILE  FOR  BLUE  OQTQ''   '*  ;N* 

Zm     PRINT  OS i "OPEN  "jN«i".eLUE" 

210  PRINT  OSi"REflD  ";^«^;",8LUE" 

220   INPUT  M 

230  FOR  I  =  1  TO  H 

240   INPUT  81 

250  FOR  J  =  1  TO  Bl 

260   INPUT  BX<I,J>:  INPUT  ev<  I  ,J> 

270  NE:<T  :  NEXT 

280     PRINT  OS; "CLOSE    '*;N«i'*.8LUE" 

2S0     FOR   I   =   1   TO  M:X8a>  =  BXCI  ,1  >:VB<:  I:-  =  BV<  1,1  vrBPa  >  =   1. 

300     PRINT  81:    FOR  JJ   =   1    TO  81:    PRINT  BXC  I  >J..i  >.8V<  I  ,JJ  >:    NEXT  JJ 

310     PRINT   "SIZE  OF  BLUE  UNIT    %!;:    INPUT  SBC  T  ' 

320  81';l  .:•  =  SB<  I  > 

330  B0<  0  ::■  =  B0<  0  >  +  S8<  I  > 

340  eR<  I  ;«  =  0 

341  •;<  '  XB<  I>:V  =  V6(  I> 

342  X4  =      INT  <X8<I>    •'    i000>   -  l>'i't'4  =      TNT   '  VB*' T  -    .-    lPift0>  -  LV 

343  ijOSUB    1 960 

344  2S<  I  ::■  =  Z 

345  PRINT  "EL£»s.»«TION  OF  BLUE  ";I;'*  13  "  ?Z 

346  FOR  II  =  1  TO  HL';X4,V4,0  >:H2  =  HL<  :<4,V4,  U  ):H«:  H2>0  >  =  i^:  NEXT 
350  NEXT 

360  PRINT 

370   INPUT  "h^HE  OF  FILE  FOR  RED  OpTQ-?  "  ;NS 

380  PRINT  OS i "OPEN  ";NS;'*. RED" 

390  PRINT  0S;"REPO  ";NS?'*.  RED" 

400   INPUT  N 

*10  FOR  I  =  1  TO  N 

420   INPUT  Rl 

■*30  FOR  J  =  1  TO  Rl 

440   INP<JT  RX(  I  ,J  >:  INPUT  RVC  I  ,.J  > 

•J.50  NEXT  :  NEXT 

•*60  PRINT  OS i "CLOSE  ";NS;".RED" 

470     FOR   I    =   1   TO  N:XR<  r  •  =  RX<  T  ,r-.-VP.n  :-   =  RV*:  I .  t  :':RP<  I:-   =    1 

430      PRINT   "SIZE  OF   RED   UNIT    ";!;:    INPUT   SR':' I  • 

490  RIU  >   =  3R<  I  :■ 

500  R0(  0  >  =  R0<  0  >   +  SR(  I  "■' 

510  RR<  IJ  =  5 

511  X   =  XR<  I>:V  =   VR<  I  > 

512  :>4  =      INT  ';XR<i::'   .'•    I000j   -  LX:V4  =      INT  •  VR<  I  -  .-■    1.000)   -  LV 

513  GOSUB    1:360 

514  ZR<  I  .:■  =  Z 

515  PRINT  "ELEUPlTION  OF  RED  'M?"  IS  ";Z 

516  FOR  II  =  1  TO  HL<X4,V4,0>:H2  =  HL-' X4  ,V4, 1 1.  ):H«r  H2  ,0  ::■  =  fi:  NEXT 
520   NEXT  I 
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53y  PRINT  :  INPUT  "00  VOI.I  WQNT  T^^:   opTv-^TPR  "iNT-  ^v.N>  "  ;PS 

535   IF  Hi  =  "V"  THEN  PR  -    1 

54ti  HOM£  :  UT(^  3 

55Hi   IF  PR  =  1  THEN  PR#  I 

5bHr  PRINT 

570  PRINT  "TIHE"i:  POKE  36.?!  PRTNT  "9!.MF"^:  ='OKE  3B^20:  PRINT  "RED" 

5Sw  PR**  (•? 

59w  IF  ZZ  =  1  THEN  itJUSUB  t?6i?»:  GOTQ  Pif^ 

Hij^i^i  GOSUB  1 490 

Hi  10  REH  MP  IN  PROGRfiH 

620  T  =  T  -t-  OT 

t?30      IF    r   <      =  NT  THEN   640 

?35     PRINT   "END  OF  BfiTTLE  DUE  Tfi   TIHF"-    C-fOSUB    1540?    ENO 

640     I3OSUB   790 

650     i5<DSU6   1170 

660  KK  =  KK   +   1 

670  80<KK)  =  0:R0<KK:'   =  0 

690     FOR    I   =    1    TO  M 

690  e0<  KK  >  =  80<  KK  ::■   ■*•  '5B<  I  > 

700     NEXT    I 

7"  10     FOR   J   =    1    TO  N 

720  R0(Kk;>  =  R0c:kK)   +  :5R<J> 

730     NEXT  J 

?40      IF  ZZ  =   1    THEN     i^OSUB   17^0?    t70T0  ^Rf* 

750     I50SU6    1490 

760      IF  80<KK:)   >  BE  *  B0<0>  THEN  770 

765  PRINT   :    PRINT   "ENO   OF  BATTLE  DUE  TO  PTTRITIOW  OF   BLUE" 

766  150SUB   1540:    ENO 

770      IF  R0<KK>    >  RE   *  R0<0>   THEN   790 

775  PRINT   :    PRINT   "ENO  OF  PfiTTLE  OUE  tq  QTTRJTTON  OP  RED" 

776  I50SUB    1540:    ENO 
790     9uT0  620 

790     REM     MOUEHENT  SUBROUTINE 

900     FOR    I   =    I   TO  N 

910  J   =  RP< I >  . 

520  OV  =  RV:  I , J  +   1  :-  -  VR<  T  :• 

'i3Q   OX  =  RXc  I  ^J   +    1  •:•   -  XR«r  I  > 

940      IF  OX   =  0  mNO   OV    >   0  THEN   QN   =  c-T    ••   ?:    r-tr'Tn  990 

950      IF  OX  =  0  QNO  OV   •••   0  THEN  mN   =  3  *  PI    ••   2?    i50T0   990 

;:!60   WN   =      PTH   •:  DV   .•■'   OX  > 

970   IF  OX  <    0  THEN  QN  =  PI  +  hN 

980  01  =  DT  *  fiR<  T  )  *  COS  •:QN':' 

990  02  =  DT  *  RRi.  I  >  *  9 IN  •:' PN ':' 

900  Oh  =  9QR  (Dl  ^  2   •+•  02  ■■•-  2  • 

910  01  =  9QR  (OX  -^  2  +  DV  -^  2-' 

920  IF  01  >  =  OH  THEN  950 

930  01  =  OX: 02  =  DV 

940  Rp^  I  :>  s  J  -i-  1 
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950  VR<  I  >  =      INT  '::VR<I>  +  02  > 
i€&  XR<  I  >  =      INT  '■■  )<R(  I  >  +  01  > 

361  X  s  ;<R<  I  >:V  =  VR<  I> 

362  X4  =      INT  <'APXl.:>  .■■-    11500  >  -  L;<:V4  =      INT  •:  VR<  T  ;•   ••'    1000)  -  LV 
'363     30SUB   1 S60 

m4  ZP.<  I  .:•  =  Z 

i«65  PRINT  "ELEUOTION  OF  RED  '*;!;"  TS  "  ;Z 

:?66  FOR  1 1  =  I  TO  ML-:  X4,V4,0  ):H2  =  WL(  X4,V4,I  I  >:H':' H2,0)  =  0;  NEXT 

370  NEXT 

360  FOR  I  =  1  TO  N 

990  J  s  eP'I) 

1000  Dv  =  8v< i,j  +  1  >  -  ve< I > 

1010  OX  =  BX< I,J  +  1 >  -  XB<  I  > 

1020   IF  DX  =  0  flNO  DV  ■>  0  THEN  PN  =  PT  /  2?  BOTO  1060 

1030  IF  OX  =  0  OHO  DV  <  0  THEN  PN  =  3  *  PT  /  2:  GOTO  1060 

1040  QN  =  ftTN  (DV  .'  0X> 

1050   IF  OX  <  0  THEN  PN  =  PI  +  QN 

1060  01  =  DT  *  BRC  I  >  *  COS  <PN> 

1070  02  =  DT  *  BR<I>  *  SIN  (flN::- 

1080  OH  =  SQR  <01  -^2  +  02  -^  2':' 

1090  01  =  SQR  COX  •'«  2  +  DV  -^  2> 

1100   IF  DI  >  =  DH  THEN  1130 

1110  01  =  OX: 02  =  DV 

1120  BR  I ::-  =  J  -f  t 

1130  VB<:I:j  =  INT  vVB<r:-  +  02  ■> 

1140  XB<  I  :•  =  INT  <X8<  D  +01  ) 

1150  NEXT 

1160  RETURN 

il7'0  REH   ATTRITION  SUBROUTINE 

1130  FOR  I  =  1  TO  M;6F'::i)  =  0!OB(r:'  =  0:  NEXT 

1130  FOR  J  =  1  TO  N:RF<J:)  =  0:OR<J:>  =  0:  NEXT 

1200  FOR  I  =  1  TO  M 

1210  FOR  J  =  1  TO  N 

1220  BT'.'  I  ,J  )  =  0;  RT<  J,  I  :-  =  0 

1 230  R6  =  SQR  ( <  X8<  I  >  -  XR<  J  >  >  -  2  >  •'  VB<  I  >  -  VR<  J  >  >  ■-•   2  > 

1240   IF  RG  :•••  2700  THEN  1310 

U50  RR<  I  >  s  2 

1251      REH     CHECK  LINE   OF  SIGHT 

[NT  ':XB<I>    •'    1000  >  -  LX:Vt    =      TNT   ••  ve<  I  >   .••■    1000;'   -  LV 
INT  '::XR<J':'   .■••    1000:<   -  lX:V2   =      INT   ':VRv.J':'   .•■■    1000>  -  LV 

X2  -  X1:V3   =  V2  -  VI 

II    =  0  TO     wBS  (X3> 

XI    -^11    *     Sj5N  <X3:> 

Jl   =0  TO     MBS  ':V3':' 

VI    +  Jl    *     SGN  '.'V3':' 
i  26 1      PR  I  NT   " CHECK  I NG  GR I D  SQUPRE   "  ; X4 ,- "  ,  "  :  V^ ; "   POR   SLUE    '* ;  I  ;  "   RED 


1252 

XI  = 

1253 

X2  = 

1254 

.^0'  — 

1255 

FOP 

1256 

X4  = 

1257 

FOR 

i253 

V4  = 

164 


1262 
1263 
1264 
1265 
1270 
1230 
i29y 

1300 

1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 

1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 

1500 

1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 

1600 

1610 
1820 
1630 

1640 

1650 
1660 
1670 
1680 
1690 


60SU6   I960 
HEAT    ;    ^€,'<T 
FOR   II    =1    TO 
REM        LOS    IS  i 


8F<I 
RF<.J>  = 
8T<I,J::' 
9J<JA> 

FOR  I 
FOR  J 
IF  RF< 

Dsa 


eF<  I  > 
RF":  J  > 

=  M  * 

=  8  * 
NEXT 
=   1    TO 
==   t    TO 

J>     =    0 

=  DBc  I  ;.' 


N-I:H<  IU0 
OK 
■»•   1 
+   1 

<  1    -  <  P  ■'  ^ 

<  1   -  <  Q  >  .x 


N 

OR  Bf( I  > 
+  BT-::  r ,  J 
>  RT-rj,! 


=  0:    NEXT 


se<i 

SRC  J 


■}  I 


*  SR(  t> 

*  S8<  I  ::■ 


0  THEN 
'  RFv.J-- 
'  8Fa  > 


1370 


Oft<.J>  =  OR<J> 

NEXT   :    NE:<T 

FOR  I  =  1  TO 
SB<T>  =  3Ba> 

IF  S8<I>  <   0 

IF  SSi.  I  ■ 

NEXT 

FOR  J  =  1  TO 
SR<J>  =  3R<j::' 

IF  3R<J>  <   0 

IF  sr<i:j 

NEXT 

RETURN 

REH  PRINTOUT.-'SUMHflRV 

IF  PR  =  1  THEN  PR#  I 

PRINT  T;:  POKE  36,7:  PRINT  80(  KK 

PR#  0 

RETURN 

REM     Si^iS  RESULTS   IN  0  FILE 


M 

-  DBC  I  ■.•- 
THEN  SB<  I  :• 

=  BM  ^  BK  I : 

N 

-  DR<J> 
THEN   SR<J::' 

=  RM  *  Rl< I 


=    0 

•  THEN 


=   0 
•  THEN 


eR<  I  >  =  8 


RR<  I  )  = 


-   \-i 


POKE  36ii20:    PRINT  R0<KK 


PRINT 

IF  «S 
OS  =     I 

INPUT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

FOR  I 
T  =    I    *   3 

PRINT   T: 

NEXT    I 


;    INPUT    "00  vnu  WCINT  P  FILE  OF  THIS  RUN^S  Ol^Tft?  •?/><>    ";PS 

=   "N"   THEN    1750 

HR«  •:4:' 

"NPME  OF  FILE  FOR  FORCE  LE»-€L  DPTQ-?  "  ;Fi 

DS;"OPEN  ";FS;'M::nEFS" 

DS ; " WR I TE  '' ;FS i  " ,  COEFS " 

Q:  PRINT  8:  PRINT  P:  PRINT  ij 

OS  i "  CLOSE  '• ;  FS  i  " .  COEFS  " 

0$ J "OPEN  ";F«J". RESULTS" 

OS ,- "  WR  I  TE  "  ;  FS ; " .  RESUL  TS  " 

NK  -f  1:  PRINT  NK  +  1 

-   0  TO  NK 


PRINT  80^ I 
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17©0  FOR  I  =  y  TO  Mk' 

1?UJ  T  =  I  *  3 

1720  PRINT  T:  PRINT  RiiK  I  > 

1730  NEXT  I 

1740  PRINT  DS;"CLOSE  " ;FS;". RESULTS" 

1750  RETURN 

1760  REH  PRINTOUT-'SUMMfiRV 

1770  HOHE  :  UTfiB  3 

1730  IF  PR  =  I  THEN  PR#  ! 

1790  PRINT  "TIME  IS  ";T 

1300  PRINT  :  PRINT 

1310  PRINT  "BLUE  FORi:ES" 

1320  PRINT  "UNIT";  TRe<  7>;"^<C00RD";  TQe<  t.'5>:"VC0nR0'* ;  Tfi8<  23>j'*SIZE" 

1330  FOR  I  =  I  TO  M 

1340  PRINT    '•    '•;!;   Tpe<    7);Xe<i:);   TPe<    t5::';VB<  I  >;   TP8<    24?;SB<I> 

1:350  NEXT 

1360  PRINT    :    PRINT   "RED  FORCES'* 

1370  PRINT   "UNIT";   TOB(    7>;"XC00R0";   T«B<    l'=i>;"VCOnRO*;   T«B<    23::'; "SIZE" 

1330  FOR    I    =    1.    TO   N 

1330  PRINT   "    ";I;   TUe<    7'-'jXR<r>;   Ti:i8<    15>;VR<I>;   TfiB<    24>;SR<I> 

1300  NEXT 

1310  PRINT  :  PRINT 

1920  PR*  0 

1330  INPUT  "HIT  RETURN  TO  CONTINUE" ;Zf 

1340  RETURN 

1350  REH  ELEUPTION  COHPUTOTION 

I960  CR  =  -  99999999 

1970  Z  =  HL<  0>0ji0.) 

1980  FOR   HI    =    1    TO   HL<X4,V4ii0  :• 

1390  H2  =  HL<  X4,V4, Hr> 

2000  IF  H<H2^0>  =   1    THEN  2150 

2010  :<S  =  :■<  "  H<H2*1>  *   100:VS  =  'J   -  H<H2.2>   *   100 

2020  wl   =     LOi?  ^•H<H2*7'>  .•-  (H(H2,7)  -  50 )  > 

2030  91    =  PI    *  <H<H2^5)>   ^   2 

2040  CI   =  H«;H2,4)  *  PI   .-'   180 

2050  C2  =  H<H2^6>  >^  2 

i060  1^1    =     -  ':P1    *  >:     COS  <  CI  •  •■  -^   2   ->•  91    *  -'    S'lN  ';C1>'   ■*■   2>    '-   C2 

2070  P2  =     -  -:  PI   *  .;■    SIN  ':C1  ••>  -  2  +  81.    ^  <    COS   ''CI  "'"'   ■-  2>   -•■  C2 

2030  P3  =  <2   *     COS  ';C1)   ^     SIN  ''Ct-   *  •:' 81    -  Pt  >  ■-  .-•  C2 

2090  IF  H<H2,7''    >  H(y^2,&')   THEN  CR   =      LOG  <'rH<H2^7:'   -  H«:H2,SV:'    ••■  H':H2^7>::. 

1100  QI    =  PI    *  X3   -^   2   -f-  P2   *  VS  -^  2   +  P3   *  XS   ^  VS 

2110  Hi:H2>0  >   =    1 

2120      IF  QI    <    CR   THEN   2150 

2130   FI    =  H<H2*3::'   +  WH2,7>   *  •:    EXP   <0l>   -    1  ■' 

2140      IF  FI    ■;:•  Z  THEN  Z   =  FT 

2150      NEXT 

2160     RETURN 
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I.   POSTPROCESSOR  PROGRAfl 
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•  Z2  Z*"0  CO  «lUih-<<  •<  X<  :tf  cotoZ  Z  ZOZO  lUt— 
CM  JbjS(/Y_JH>  ♦VJh- ^  CO  COOCQ.  -coo  •  w^  •  ^^K  O  O  O  -^  ►— ^ 
O  aCtUO  '-iUJUYX  <  <  >  ■•V-COt/K  •(%!  aO^  aUiO  UJOUJOUJCsJLU  z*-* 
^0»  —tJUJ*!!*—*  X  •♦-UJO<«JUja:Q._J  -Z^— »H-XCvJXfMX  XO  •-• 
m-J  XXH-X3H-+-JH-  h-t-3XX  OXQl<Z  -  -^Z-JCDOI-  h-  ►-a»-»-  z 
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>>aiujujZu.LLU.a:zuj        at        cou.u_LLLLou.LLU->u-Luu.u_u.a.u-<  "t/i  "cov-ToX        —«»-< 
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a.      COUJ  C3C  X  LL.  U_  U_  O-JO—iO— iO-jOQi-oCUU. 
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a.  I—  o  Q  Q  <H-<i— <»-<i— -"jLuxoa 
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